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Final  Report 

Coherent  Structures  and  Chaos  Control  in 
High-Power  Microwave  and  Charged-Particle  Beam  Devices 
AFOSR  Grant  No.  FA9550-06-1-0269 

This  report  summarizes  research  results  obtained  under  the  auspices  of  Air  Force 
Office  of  Scientific  Research,  Grant  No.  F49620-06-1 -0269  (Chen,  2006).  In  particular, 
we  conducted  vigorous  theoretical  and  computational  investigations  of  coherent 
structures  and  chaos  in  a  wide  range  of  intense  electron  beam  plasmas  relevant  to  the 
development  of  high-power  microwave  and  particle-beam  devices  for  directed  energy 
applications. 

The  following  is  a  brief  summary  of  our  research  accomplishments  in  selected  areas, 
while  detailed  findings  are  described  in  the  preprints  cited  in  this  report. 

1.  High-Power  Magnetron  Research  (Zhou  and  Chen,  2007  and  2008;  Davies,  Zhou 

and  Chen,  2007;  Davies,  Chen  and  Zhou,  2008) 

Under  the  auspices  of  the  grant  (Chen,  2006),  we  have  developed  a  small-signal 
theory  of  a  non-relativistic  magnetron  using  a  planar  model  with  a  thin  electron  cloud 
(Zhou  and  Chen,  2007  and  2008).  The  theory  includes  both  inertial  effects  and 
electromagnetic  effects  in  a  Floquet  expansion.  We  have  derived  an  analytical  dispersion 
relation  of  such  a  planar  magnetron,  and  calculated  the  growth  rate  analytically.  We  have 
shown  that  the  magnetron  instability  involves  the  resonance  between  the  electron  cloud 
and  the  slow  waves  in  the  magnetron  cavities.  We  have  found  good  agreement  between 
the  theory  and  the  self-consistent  two-dimensional  (2D)  particle-in-cell  (PIC)  MAGIC 
simulations.  We  have  predicted  vortex  structures  in  the  equilibrium  relativistic  electron 
tlow  in  magnetrons  (Davies,  Zhou  and  Chen,  2007;  Davies,  Chen  and  Zhou,  2008).  The 
vortex  structures  are  induced  by  the  periodic  corrugations  on  the  magnetron  anode.  In  the 
analysis,  we  have  made  the  guide-center  approximation,  which  is  validate  at  low  electron 
densities.  We  have  validated  the  analysis  using  test-particle  calculations. 

1.1  Small-Signal  Gain  Theory  of  Non-Relativistic  Planar  Magnetrons  (Zhou  and 

Chen,  2007  and  2008) 

The  onset  of  the  unstable  oscillations  in  magnetrons  has  not  been  analytically 
described  to  complete  satisfaction,  although  extensive  particle-in-cell  (PIC)  simulations 
can  make  good  predictions  for  the  instability  characteristics  (see,  for  example,  Chan, 
Chen  and  Davidson,  1993;  Lemke,  Genoni  and  Spencer,  1999).  Previous  analytical 
studies  included  various  models  utilizing  linear  theories.  Earlier  work  focused  on  the 
diocotron  instability  in  the  guiding-center  approximation  (Davidson,  Chan,  Chen,  et  al., 
1991;  Ayres,  Chen,  Stark,  et  al.,  1992),  which  ignores  inertial  effects  in  the  electron  cloud. 
The  recent  work  by  Riyopoulos  on  the  basis  of  a  guiding-center  model  provided  new 
insight  into  the  magnetron  instability  in  the  low-space-charge  limit  (Riyopoulos,  1998).  A 
linear  theory  taking  a  single  rf  mode  in  the  Floquet  expansion  was  developed  to  include 
electromagnetic  effects  (Kaup,  2001,  2004  and  2007).  Despite  these  theoretical  and  PIC 
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simulation  efforts,  quantitative  agreement  between  theory  and  PIC  simulations  has  not 
been  reported  until  our  paper  (Zhou  and  Chen,  2007  and  2008). 


Figure  1.1.1  Schematic  of  a  planar  magnetron  with  a  thin  electron  layer  (from  Zhou  and 
Chen,  2008). 


Assuming  a  multi-stream  thin  electron  cloud  located  at  a  distance  of  h  from  the 
cathode  as  shown  in  Fig.  1.1.1,  we  have  obtained  the  small-signal  gain  equation  or  the 
loaded  dispersion  relation  (Zhou  and  Chen,  2008) 
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.  For  £k  =0,  Eq.  (1.1.1)  is  the  vacuum  dispersion  relation  for  the 


corrugated  stmeture.  For  ek  ^0,  on  the  other  hand,  Eq.  (1.1.1)  permits  calculations  of 

small-signal  gains  in  magnetrons,  which  will  be  discussed  in  Sec.  1.2. 

It  should  be  noted  that  the  simplifying  assumption  of  a  thin  electron  cloud  has 
enabled  us  to  avoid  the  difficulties  in  treating  multiple  poles  in  the  small-signal  gain 
equation  in  general  situation.  We  plan  to  further  study  the  general  small-signal  gain 
equation. 


1.2  Comparison  between  Small-Signal  Gain  Theory  and  Two-Dimensional  Particle- 
in-Cell  Simulations  (Zhou  and  Chen,  2008) 
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To  compare  our  small-signal  gain  theory  with  2D  MAGIC  simulations,  we  have 
restricted  to  the  low-eurrent  regime,  where  the  thm-beam  equilibrium  model  is  a  good 
approximation  to  the  thin  electron  eloud  in  the  2D  MAGIC  simulations. 

As  an  example,  we  have  considered  a  system  with  the  parameters:  L  =  0.478  cm, 
a  =  0.382  em,  b  =  0.478  em,  d  =  4.25  cm,  h  =  0.382  em,  Bz  =  1 80  G,  Ex  =  -5.27  kV/em  , 

crb  -  2. 1 2  x  1 08  cm  2  ,  and  Vb  =  0.098c . 

First,  we  have  computed  the  vaeuum  dispersion  relation  from  Eq.  (1.1.1)  with 
ek  =0,  and  found  good  agreement  with  the  2D  MAGIC  simulation  with  the  absence  of 

the  electron  eloud,  as  shown  in  Fig.  1.2.1. 

Second,  we  have  made  use  of  the  self-consistent  PIC  eode,  2D  MAGIC,  to  simulate 
the  planar  magnetron  system.  Because  2D  MAGIC  can  handle  only  a  few  vanes,  a  3-vane 
slow-wave  corrugated  structure  with  the  same  parameters  as  in  Fig.  1.2.1  is  used  in  the 
simulation.  Penodie  boundary  conditions  are  used  sueh  that  the  2/r/3  mode  is  supported 
by  the  3-vane  structure.  The  uniform  crossed  eleetrie  and  magnetic  fields  are  applied  with 
Bz  -  1 80 G  and  Ex  =-5.27  kV/em.  The  electron  beam  is  initialized  as  a  slab  infinite 
long  in  the  z  direction  and  with  a  width  of  0.05  em  in  the  x  direction.  The  electron  beam 
propagates  with  an  initial  velocity  of  Vb  -  0.098c .  As  the  beam  propagates,  the  instability 

starts  to  build  up  which  is  illustrated  by  the  voltage  cross  the  vane  tip  of  the  slow-wave 
structure.  In  Fig.  1.2.2,  the  amplitude  of  the  oseillating  voltage  filtered  by  a  filter  that 
selects  the  2/r/3  mode  is  plotted.  In  the  early  stage  of  the  instability  (e.g.  for  /  <  15  ns), 
the  oscillation  has  a  very  small  amplitude  and  exhibits  a  relatively  broad  frequency 
speetrum  whieh  is  not  shown  in  Fig.  1.2.2.  Starting  from  t  =  10  ns,  the  2/r/3  mode  grows 
exponentially.  It  saturates  at  about  t  =  26  ns.  For  this  MAGIC  simulation,  the  2/r/3 
mode  is  determined  to  have  a  frequency  of  1.93  GHz  and  an  amplitude  growth  rate  of 
3.12  dB/cm,  shown  as  eirele  and  eross  in  Fig.  1 .2.3,  respectively. 

Finally,  we  have  solved  the  loaded  dispersion  relation  in  Eq.  (1.1.1)  with  £k  ^0  to 

calculate  the  real  frequency  and  the  instability  growth  rate.  Due  to  the  intluenee  of  the 
anode  corrugation,  the  electron  cloud  velocity  prior  to  the  linear  growth  is  observed  to 
vary  sinusoidally  in  the  y  -  direction  with  small  amplitude  around  the  averaged  flow 
velocity  in  the  2D  MAGIC  simulation.  To  model  the  velocity  variation,  we  have  used 
three  electron  cloud  streams  at  the  same  location  x  =  0.0384  cm,  each  with  one  third  of 
the  total  surfaee  eharge  number  Gb\  =  crb2  =  <Jb?>  =  <J&/3  and  a  slightly  different 
velocities,  i.e.,  Vbl  =0.98F),,  Vb2  -  Vb ,  and  Vh3  =  l.02Fb .  The  theoretical  growth  rate  is 
in  good  agreement  with  the  2D  MAGIC  simulation  as  shown  in  Fig.  1 .2.3. 

1.3  Vortex  Structures  in  Relativistic  Magnetrons  (Davies,  Zhou  and  Chen,  2007; 

Davies,  Chen  and  Zhou,  2008) 

The  periodic  corrugations  on  the  anode  have  a  strong  influence  on  the  equilibrium 
flow  in  magnetrons  and  erossed-field  deviees.  In  particular,  we  predieted  the  existence  of 
vortex  structures  induced  by  the  corrugations  on  the  anode  in  non-relativistic  magnetrons 
(Chen,  2004;  Bhatt  and  Chen,  2003  and  2004;  Davies  and  Chen,  2006).  Until  our 
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prediction,  vortex  formation  had  not  been  discussed  in  the  literature  in  the  context  of 
magnetrons  and  crossed  field  devices.  In  fact,  all  conventional  treatments  of  the 


Figure  1.2.1  Vacuum  dispersion  diagram  fora  planar  magnetron  structure  obtained  from 
the  dispersion  relation  in  Eq.  (1.1.1)  with  £k  =  0  (solid  curve)  and  2D  MAGIC 

simulations  (dotted  curve).  The  dashed  curve  is  the  beam  line  /  =  kz/27TVh  .  Here,  the 
parameters  are  L  -  0.478cm,  a  =  0.382  cm,  b  =  0.478  cm,  d  =  4.25cm,  and  Vb  -  0.098c 
(from  Zhou  and  Chen,  2008). 


Figure  1 .2.2  Plot  of  the  amplitude  of  oscillating  voltage  as  a  function  of  time  as  obtained 
from  the  MAGIC  simulation.  Here,  the  parameters  are  L  =  0.478  cm,  a  =  0.382  cm. 
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b  =  0.478cm,  d  =  4.25cm,  h  =  0.382cm,  5.  =  180G,  Ex  =  -5.27  kV/cm, 
<jh  =  2. 12  x  10s  cm  2 ,  and  Vb  =  0.098c  (from  Zhou  and  Chen,  2008). 


Figure  1.2.3  Plot  of  the  temporal  growth  rate  Im(&>)  (solid  curve)  and  the  real  frequency 
Re(&>)  (dashed  curve)  as  a  function  of  wave  number  kz  for  the  lowest  resonant  TM 
mode.  The  solid  and  dashed  curves  are  obtained  from  Eq.  (1.1.1)  using  the  three  electron 
cloud  streams,  whereas  the  circle  and  cross  are  the  real  frequency  and  instability  grow 
rate  from  the  MAGIC  simulation,  respectively.  Here,  the  parameters  are  L  =  0.478  cm, 
a  =  0.382  cm,  b  =  0.478  cm,  d  =  4.25  cm,  //  =  0.382  cm,  5r=180G,  Ex  =  -5.27  kV/cm , 

<rw  =  <rh2  =  <rhi  =  <Jh/3  =  7.07  x  10’  cm  Vb]  =  0.98 Vh ,  Vb2  =  Vb  =  0.098c  ,  and 
Vb$  =  1.02F^  (from  Zhou  and  Chen,  2008). 

equilibrium  electron  flows  in  magnetrons  had  ignored  effects  of  the  anode  corrugations 
(see,  for  example,  Davidson,  Chan,  Chen,  et  ah,  1991). 

By  definition,  an  equilibrium  flow  corresponds  to  a  state  in  magnetrons  and  crossed- 
field  devices  in  the  absence  of  rf  oscillations.  While  it  is  an  ideal  situation  in  an  actual 
device,  it  is  of  critical  importance  in  order  for  us  to  develop  a  better  understanding  of 
magnetrons  and  crossed-field  devices. 

Since  our  report  of  vortex  formation  (Chen,  2004)  at  Magnetron/PIC  Simulation 
Workshop  held  at  AFRL,  May  2-3,  2004,  several  researchers  at  AFRL  (Cartwright,  2007) 
and  elsewhere  (Bosman,  et  ah,  2005)  have  been  exploring  use  of  non-axisymmetric 
cathodes  to  improve  magnetron  performance. 

Under  the  auspices  of  the  grant  (Chen,  2006),  we  have  established  the  theory  of 
vortex  structures  in  equilibrium  electron  flows  in  relativistic  magnetrons  (Davies,  Zhou 
and  Chen,  2007;  Davies,  Chen  and  Zhou,  2008),  as  a  generalization  of  our  earlier  work 
(Davies  and  Chen,  2006).  The  theory  employs  a  planar  geometry  shown  in  Fig.  1.3.1  and 
the  guiding-center  approximation  which  is  valid  when  the  electron  plasma  frequency  is 
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small  compared  to  the  cyclotron  frequency.  Under  the  guiding-center  approximation,  the 
(nonlinear)  equilibrium  equations  are 

V2^(jc,y)  =  4xen(x,y),  (1.3.1) 


y 


Figure  1.3.1  Electron  flow  under  the  influence  of  anode  corrugations  in  a  relativistic 
planar  magnetron.  Here,  the  y  axis  is  the  cathode,  the  curve  x  =  x6(y)  is  the  electron 

layer  envelope,  and  the  dashed  curve  is  a  representative  anode  (Davies,  Zhou  and  Chen, 
2007). 


4  jte 


e.  x  VB:(x,y)  = - /?(.v,y)v(x,>-), 


V(*,>’)  = 


K  x  V<f>(x,y), 


(1.3.2) 

(1.3.3) 


where  —e  is  the  electron  charge,  n{x,y)  is  the  electron  number  density,  V(x,y)  is  the 
electron  fluid  velocity,  and  <f>{x,y )  is  the  electric  potential.  The  electric  field  E(.v,  y)  is 
given  by  E(x,y)=  -Vj(x,y). 

For  a  uniform  electron  layer,  the  solutions  to  the  nonlinear  equilibrium  equations 
(1.3.1 )-( 1.3.3)  are 
</>(x,y) 


-  {E,  (0. y))x  +  2mi0ex 2  +  jr  Fn (eKV  -  e~^x  Jcosfoy) 

n-\ 

V0  +  C0(.r-  D)+  Vpe~^x~D^  cos(/cy) + y  AneK"D  -  e~*"  )cos(ir„y) 

n= 1 

b2(x  y )  =  '  8^[^(-Ty)-  ^(•^(y),>')]+  Bl  (o  <  x  <x J 
W  (xh<x<xa) 


(0<  x<xh) 

(xh<x<xu) 

(1.3.4) 

(1.3.5) 
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where  the  coefficients  are  given  by  analytical  expressions.  One  important  result  is  that  the 
solution  given  in  P2qs.  (1.3.4)  and  (1.3.5)  supports  relativistic  vortex  structures,  as  shown 
in  Figs.  1.3.2  and  1.3.3.  The  other  important  result  is  that  the  constant  magnetic  field 
contours  coincide  with  the  equipotential  contours,  as  seen  in  Eq.  (1.3.5)  and  Fig.  1.3.3. 

As  a  numerical  example,  we  have  considered  a  system  with  the  following  parameters: 

V0  =  5 1 2  k V  ,  Vp  =  49.049  kV  ,  50=20kG  ,  (Ex(0,y))  =  0  ,  k  =  2n  I L  =  n  cmA 

(corresponding  to  £  =  2.0  cm),  and  D  =  0.5cm.  The  geometry,  magnetic  field  5n,  and 
potential  V0  of  this  example  provide  a  rectangular  approximation  to  a  typical  cylindrical 
L-band  relativistic  magnetron  (see,  for  example,  Lemke,  Genoni  and  Spencer,  1999). 
Instead  of  directly  specifying  a  value  of  n0 ,  we  have  chosen  cop/coc  - 1/3  ,  where 

o)p  =  (4mj0e2  hn)1  is  the  nonrelativistic  plasma  frequency  in  the  electron  layer,  and 
coi  -eB^hnc  is  the  nonrelativistic  cyclotron  frequency  in  the  vacuum  region.  From  the 
values  of  B0  and  cop !coc  specified  above,  it  follows  that  n0  =  4.3200  x  1012 . 

Computed  contours  of  constant  potential  are  shown  in  Fig.  1.3.2.  The  cathode  at 
potential  ^  =  0  is  represented  by  the  line  x  =  0.  The  equipotential  contour  at  512  kV 
shown  by  a  dashed  curve  corresponds  to  the  corrugated  conducting  anode.  The  electron- 
layer  boundary  is  represented  by  the  second  dashed  curve  at  93  kV.  The  error  of  ±  1  kV 
in  its  potential  of  (f)(xb(y\y )  =  93  kV  is  due  to  the  neglect  of  higher  order  terms  than  the 
first  in  the  perturbation  calculation. 

Also  shown  in  Fig.  1.3.3  are  the  test  particle  trajectories  which  follow  the 
equipotential  contours.  This  is  a  numerical  validation  of  our  vortex  theory  in  equilibrium 
electron  flows  of  relativistic  magnetrons. 


Figure  1.3.2  Contours  of  constant  potential  as  functions  of  jc  and  v  for  a  system  with 
parameters  V0  =  5  1 2  kV ,  Vp-  49.049  kV  ,  k  =  n  cm  1 ,  (SEX  (0,  y  ))  =  0  ,  B0  =  20  kV  , 
D  =  0.5  cm  ,  and  cop/coc  =  1/3  .  The  corrugated  anode  and  the  electron-layer  boundary 
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(also  contours  of  constant  potential)  are  depicted  by  the  dashed  curves  at  93  kV  and  5 1 2 
kV,  respectively  (from  Davies,  Chen  and  Zhou,  2008). 


Figure  1 .3.3  Detail  of  the  electron-layer  region  for  the  system  of  Fig.  1 .3.2.  Contours  of 
constant  potential  are  shown  by  the  solid  lines  with  the  electron-layer  boundary  shown  by 
the  dashed  curve  at  93  kV.  Because  these  contours  are  also  lines  of  constant  magnetic 
field,  each  contour  is  labeled  with  both  a  value  of  (f)  and  of  B  .  Trajectories  of  test 
particles,  each  launched  from  the  point  marked  with  an  “x”  ,  are  shown  by  circles  (from 
Davies,  Chen  and  Zhou,  2008). 

Our  discovery  of  vortex  structures  in  magnetrons  and  crossed-field  amplifiers 
provides  new  insight  into  a  number  of  areas  in  magnetron  and  CFA  research,  including 

a)  Elimination  of  ion  trapping  and  ion  noise, 

b)  Elimination  of  turbulence  and  noise, 

c)  Improvement  of  magnetron  efficiency  (Cartwright,  2007),  and 

d)  Suppression  of  mode  competition  (Bosman,  et  al,  2005). 

We  plan  further  investigate  them. 

2.  Discovery  of  Adiabatic  Thermal  Beam  Equilibrium  in  a  Periodic  Solenoidal 

Magnetic  Field  (Samokhvalova,  Zhou  and  Chen,  2006,  2007a  and  2007b;  Zhou, 

Samokhvalova  and  Chen,  2006  and  2008;  Samokhvalova,  2008) 

Many  HPM  devices  such  as  klystrons  and  TWTs  employ  a  periodic  solenoidal  (or 
permanent)  magnetic  focusing  field.  Under  the  auspices  of  the  present  grant  (Chen,  2006), 
we  have  discovered  the  thermal  equilibrium  state  of  an  electron  beam  as  it  undergoes 
adiabatic  expansion  and  compression  processes  in  a  periodic  or  axially  varying  solenoidal 
magnetic  field.  In  particular,  we  have  developed  paraxial  kinetic  and  warm-fluid 
equilibrium  theories  for  a  thermal  electron  beam  (Samokhvalova,  Zhou  and  Chen,  2006, 
2007a  and  2007b;  Zhou,  Samokhvalova  and  Chen,  2006  and  2008).  While  our  kinetic  and 
warm-fluid  beam  equilibrium  theories  are  equivalent  and  applicable  to  for  both  the  non- 
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relativistic  and  relativistic  regimes,  to  our  warm-fluid  equilibrium  theory,  our  kinetic 
beam  equilibrium  theory  provides  detailed  information  about  the  electron  distribution  in 
phase  space.  We  have  found  good  agreement  between  our  theories  and  the  recent 
experimental  measurements  of  the  electron  density  profiles  at  University  of  Maryland 
Electron  Ring  (UMER)  (Bernal,  Quinn,  Reiser,  et  al.,  2002;  Zhou,  Samokhvalova  and 
Chen,  2008).  Our  discovery  of  the  adiabatic  thermal  beam  equilibrium  has  filled  a  major 
gap  in  the  understanding  of  periodically  focused  beam  equilibria. 

2.1  Kinetic  Beam  Equilibrium  Theory  (Zhou,  Samokhvalova  and  Chen,  2006  and 
2008;  Samokhvalova,  2008) 

In  our  kinetic  beam  equilibrium  theory  (Zhou,  Samokhvalova  and  Chen,  2006  and 
2008;  Samokhvalova,  2008),  we  have  considered  a  thin,  continuous,  axisymmctric 
(5/50  =  0),  single-species  charged-particle  beam,  propagating  with  constant  axial 
velocity  Ke:  through  an  applied  periodic  solenoidal  magnetic  focusing  field.  The  applied 
periodic  solenoidal  focusing  field  inside  the  beam  can  be  approximated  by 

B“'  (r,  s)  =  -  i  B\  {s)r  cr  +  B,(s)e„  (2.1.1) 

where  5  =  2  is  the  axial  coordinate,  r  =  ^x2  +  y2  is  the  radial  distance  from  the  beam 
axis,  the  prime  denotes  the  derivative  with  respect  to  s ,  and  Bz(s)  is  the  axial  magnetic 
field  which  can  be  either  periodic  along  the  2 -axis  with  periodicity  length  S  or  an 
arbitrary  function  of  .9 .  In  the  paraxial  approximation,  rbrms  «  S  is  assumed,  where  rbrms 
is  the  rms  beam  envelope.  The  transverse  kinetic  energy  of  the  beam  is  assumed  to  be 
small  compared  with  its  axial  kinetic  energy,  i.e.,  |K,( » \VL\ ,  and  the  Budker  parameter 

of  the  beam  is  assumed  to  be  small,  i.e.,  q2Nb/mc 2  «y\f$l,  where  q  and  m  are  the 
particle  charge  and  rest  mass,  respectively,  c  is  the  speed  of  light  in  vacuum, 

X 

Nb  -2 n^drrnb{r,s)  is  the  number  of  particles  per  unit  axial  length,  and  yb  is  the 

0 

relativistic  mass  factor,  which,  to  leading  order,  is  yb  =  const  =  (l  -  Pb ) '  with 

Pb=Vb/c  =  VJc . 

Because  the  beam  is  axisymmetric,  the  canonical  angular  momentum  P0  is  a  constant 
of  motion,  i.e., 

dP 

^L  =  0.  (2.1.2) 

ds 

After  performing  a  two-step  canonical  transformation,  we  have  also  found  that  the  sealed 
transverse  Hamiltonian  for  the  single-particle  motion 

E  =  w2(s)H±(x,y,Px,Py,s)  (2.1.3) 

is  an  approximate  invariant  (for  detailed  analyses  and  definitions  of  H  ,  .v, 

/^,etc.,  see  Zhou,  Samokvalova  and  Chen,  2008).  We  have  chosen  the  beam  equilibrium 
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distribution  in  the  form  similar  to  the  Maxwell-Boltzmann  distribution  (Zhou, 
Samokhvalova  and  Chen,  2008;  Samokhvalova,  2008) 

fb  =  Cexv[-p{E-cobPg)\  (2.1.4) 

where  C,  /?  and  coh  are  constants.  C  is  an  integration  constant,  /?  is  related  to  the  beam 
emittance,  and  cob  is  the  rotation  frequency  relative  to  the  Larmor  frame.  The  distribution 
function  fb  defined  in  Eq.  (2.1.4)  is  a  Vlasov  equilibrium,  i.e.,  dfb/ds  =  0.  The  kinetic 

beam  equilibrium  theory  makes  the  following  two  important  predictions: 

1 .  The  thermal  beam  emittance  is  a  constant. 


S,u  =  ■ 


py 


(x2)({»x~kY)=- 


_  ypAiYLM  - 


2/2 


2yhmp;c 


const , 


(2.1.5) 


2.  The  equation  of  state  is  adiabatic,  i.e.,  T± (s)rbrms  (s)  =  const . 

The  rms  envelope  equation  (2.2.4)  is  (Zhou,  Samokhvalova  and  Chen,  2008) 


-  Yl  ^  -  %T  [Q/,  (-0  +  UKm,  (s)  ~  -  ~  ~  3 

ds  Phc  2  rhms(s)  rhrms{s ) 


4c2 


(2.1.6) 


where  cohCj Phc /2rbrms{s)- flc(.v)/2  and  Q.  (s)-qB_(s)//hmc.  The  beam  density 

profile  is 


—  4^C^( 


>h(r’s)=  ™2(s)\\fdPxdry  = 


r 2 

brms 


th_ 

,(s) 


exp)  - 


A"+  4  si 


2  ^L(-0. 


<  ylkBTL{s) 


(2.1.7) 

where  the  scalar  potential  for  the  self-electric  field  is  determined  by  the  Poisson  equation 

1 

r  dr 


dr 


=  -4nqnh{r,s). 


(2.1.8) 


It  should  be  pointed  out  that  because  the  derivation  of  the  theory  does  not  assume 
specific  magnetic  profile  as  defined  in  Eq.  (2.1.1)  it  is  valid  not  only  for  the  periodic 
solenoid  magnetic  field  but  also  for  an  arbitrary  varying  solenoid  magnetic  field. 
Therefore,  our  results  apply  for  the  periodic  focusing  channel  as  well  as  for  the  matching 
section  between  the  electron  source  and  the  periodic  focusing  channel,  which  is  important 
in  the  design  of  electron  beams  for  HPM  applications  and  beam  experiments  such  as 
University  of  Maryland  Electron  Ring  (Bernal,  Quinn,  Reiser,  et  al.,  2002). 


2.2  Warm-Fluid  Beam  Equilibrium  Theory  (Samokhvalova,  Zhou  and  Chen,  2006, 
2007a  and  2007b;  Samokhvalova,  2008) 


We  have  recovered  the  results  of  the  macroscopic  qualities  in  the  kinetic  beam 
equilibrium  theory  [i.e.,  Eqs.  (2.1.5)-(2.1.8)]  by  solving  the  following  adiabatic  warm- 
fluid  equilibrium  equations  (Samokhvalova,  Zhou  and  Chen,  2007b;  Samokhvalova, 
2008) 


nbV  ■V{ybm\)=nbq 


-Vfel/  +  —  x(B"'  +B 

C 


seif 


-V-P(x),  (2.2.1) 


v-(»*v)  =  o, 


(2.2.2) 
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S72fe'f  (r,s)=  -4nqnb{r,s), 
pl{r,s)  =  nb{r,s)kBTL(s), 
*i($)iL(s)  =  const . 


(2.2.3) 

(2.2.4) 

(2.2.5) 


The  detailed  analysis  is  available  in  our  paper  (Samokhvalova,  Zhou  and  Chen,  2007b) 
and  Samokhvalova’s  doctoral  thesis  (Samokhvalova,  2008). 

2.3  Comparison  between  Theory  and  Experiment  (Zhou,  Samokhvalova  and  Chen, 

2008;  Samokhvalova,  2008) 

Using  our  adiabatic  thermal  beam  equilibrium  theories,  we  have  replicated  the  beam 
density  profiles  at  different  axial  distances  in  good  agreement  with  the  experimental 
measurements  conducted  on  the  University  of  Maryland  Electron  Ring.  Our  equilibrium 
theory  is  applicable  to  this  experiment  from  the  anode  aperture  to  a  distance  prior  to  the 
wave  breaking  initiated  by  high  order  density  distribution  fluctuations  induced  by  a 
pressure  force  at  the  anode  aperture.  Wave  breaking  occurs  at  about  one  quarter  of 
plasma  wavelength,  which  is  about  30  cm  in  this  example.  Our  equilibrium  theories  do 
not  explain  the  density  distribution  distortion  in  the  present  form,  but  it  is  possible  to 
develop  a  perturbation  theory  based  on  the  equilibrium  in  the  future. 

By  solving  Eqs.  (2.16)-(2.1.8),  we  have  calculated  the  beam  transverse  density 
profiles  of  the  UMER  5  keV,  6.5  mA  electron  beam  at  three  axial  distances:  *s*  =  6.4  cm, 

1 1.2  cm,  and  17.2  cm,  as  shown  in  solid  curves  in  Fig.  2.3.1.  The  dashed  curves  are  the 
equivalent  Kapchinskij-Vladirmirskij  (KV)  beam  density  profiles  (Kapchinskij  and 
Vladirmirskij,  1959).  Compared  with  the  experimental  measurements  (dotted  curves),  the 
calculated  beam  density  profiles  are  in  good  agreement.  As  the  beam  radius  increases,  the 
beam  density  profile  approaches  to  the  KV  (uniform)  beam  density  profile,  because  the 
beam  temperature  must  decrease  in  order  to  keep  T±  {s)r^rms  (5)  at  a  constant.  In  this 

adiabatic  process,  the  Debye  length  XD  =  yJ/lkBTL (s)/ 4 7icfnh (0, 5)  =  0.54  mm  is  constant. 

3.  Design  of  High-Brightness  Circular  Electron  Beams  (Bemis,  Bhatt,  Chen  and 

Zhou,  2007a  and  2007b) 

An  experimental  demonstration  of  the  thermal  beam  equilibrium  over  long 
propagation  distances  requires  a  high-brightness  circular  electron  beam  which  is  well 
matched  into  a  periodic  solenoidal  magnetic  focusing  channel.  To  gain  experience  in  the 
design  of  high-brightness  electron  beams,  we  have  developed,  under  the  auspices  of  the 
present  grant,  a  method  for  the  design  of  a  high-brightness  non-relativistic  circular  beam 
system  including  a  charged-particle  emitting  diode,  a  diode  aperture,  a  circular  beam 
tunnel,  and  a  focusing  magnetic  field  that  matches  the  beam  from  the  emitter  to  the  beam 
tunnel.  The  applied  magnetic  field  has  been  determined  by  balancing  the  forces 
throughout  the  gun  and  transport  sections  of  the  beam  system.  The  method  has  been 
validated  by  three-dimensional  simulations. 

While  the  detailed  method  is  described  in  our  paper  (Bemis,  Bhatt,  Chen  and  Zhou, 
2007b)  and  our  US  patent  application  (Chen,  Bemis,  Bhatt  and  Zhou,  2007),  Figs.  3. 1-3.4 
summarize  the  results  for  a  high-perveance  electron  beam  at  a  voltage  of  2.3  kV.  This 
system  is  scalable  to  moderately  high  voltages. 
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Figure  3.1  shows  the  final  results  of  an  OMNITRAK  simulation  for  a  circular 
electron  beam  which  is  emitted  from  an  emitter  (cathode)  with  a  radius  of  1.52  mm,  a 
current  of  0.1 1  A,  a  cathode-to-anode  distance  of  4.1 1  mm  at  radius  r  =  1.52  mm.  The 
diode  voltage  is  2300  V. 

Figure  3.2  shows  plots  of  Er  versus  z  at  r  =  1.0  mm  and  5.(0,z)  versus  z  from  the 
OMNITRAK  simulation  after  two  iterations.  The  radial  electric  field  vanishes  at  the 


Figure  2.3.1  Normalized  beam  transverse  density  profiles  of  a  5  keV,  6.5  mA 
(4s~ rms  =  30  mm-mrad)  electron  beam  at  three  axial  distances:  s  =  6.4  cm,  1 1.2  cm,  and 
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17.2  cm.  The  solid  curves  are  from  theory,  the  dotted  curves  are  the  experimental 
measurements,  and  the  dashed  lines  are  the  equivalent  KV  beam  density  distributions. 
The  densities  are  normalized  to  the  equivalent  KV  beam  density  (from  Zhou, 
Samokhvalova  and  Chen,  2008). 


0  z  (mm)  10.0 


Figure  3.1  OMN1TRAK  simulation  of  the  dynamics  of  an  electron  beam  emitted  from  a 
flat  circular  cathode  with  a  radius  of  1 .52  mm  and  a  current  of  0.1  1  amperes  in  an 
optimized  magnetic  field.  Here,  the  cathode-to-anode  distance  is  4.1 1  mm  at 
r  =  1 .52  mm,  the  circular  anode  aperture  has  a  radius  of  1 .8  mm,  and  the  beam  tunnel  has 
a  radius  of  6.0  mm.  The  diode  voltage  is  2300  V  (from  Bemis,  Bhatt,  Chen  and  Zhou, 
2007b). 
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Figure  3.2  Plots  of  Er  versus  2  at  r  =  1  mm  (solid  curve)  and  £z(0,z)  versus  z  (dashed 
curve)  from  the  OMNITRAK  simulations  after  two  iterations  (from  Bemis,  Bhatt,  Chen 
and  Zhou,  2007b). 
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Figure  3.3  Plot  of  the  electron  distribution  at  z  =  8  mm  in  the  phase  plane  (jt, y)  (from 
Bemis,  Bhatt,  Chen  and  Zhou,  2007b). 

emitter  (i.e.,  at  z  =  0)  ,  achieves  a  maximum  magnitude  at  z  =  4.11  mm,  and  then 
approaches  to  a  constant  value  well  inside  the  beam  tunnel.  The  applied  axial  magnetic 
field  at  the  emitter  vanishes,  increases  to  640  G  at  the  aperture,  and  then  falls  to  about 
160  G  well  inside  the  beam  tunnel. 

Figure  3.3  shows  the  electron  distribution  in  the  phase  plane  (x,y)  at  z  -  8  mm.  The 
beam  distribution  maintains  transverse  uniformity.  Indeed,  the  normalized  fourth  moment 

if  which  is  equal  to  4/3  for  a  transversely  uniform  density  distribution,  remains 

4/3  within  ±1%.  The  outer  beam  radius  remains  the  same  as  the  emitter  radius  within  1 
to  3%.  The  beam  is  very  bright,  and  its  normalized  rms  emittance  is  predicted  to  be  0.33 
mm-mrad,  which  is  1.15  times  the  intrinsic  normalized  rms  emittance  of  0.33  mm-mrad. 

Our  design  method  is  limited  to  a  nonrelativistic  beam  matching  into  a  uniform 
magnetic  field  We  plan  to  extend  our  design  method  to  include  thermal  effects  and  beam 
matching  into  periodic  solenoidal  (permanent)  magnetic  focusing  field.  We  will  also 
extend  our  design  method  to  relativistic  electron  beams  in  which  both  relativistic  effects 
and  self-magnetic  field  effects  are  important. 

4.  Development  of  Elliptic  Beam  Theory  for  High-Power  Microwave  Device 

Applications  (Bhatt,  2006;  Zhou,  2006;  Chen  and  Zhou,  2007;  Zhou,  Bhatt  and 
Chen,  2008) 
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Existing  HPM  devices  employ  either  a  pencil  electron  beam  or  an  annular  electron 
beam.  One  of  the  main  objectives  in  high-power  microwave  (HPM)  research  has  been  to 
develop  innovative  science  and  technology  which  could  lead  to  improve  the  efficiency, 
output  power  level,  and  other  performance  (such  as  weight  reduction)  of  high-power 
microwave  (HPM)  devices.  A  promising  approach  is  to  use  a  large-aspect-ratio  elliptic 
electron  beam  rather  than  the  conventional  pencil  or  annular  electron  beam  as  the  energy 
source.  Compared  with  a  conventional  HPM  device,  a  HPM  device  powered  by  an 
elliptic  electron  beam  has  the  following  attractive  features: 

a)  Power  scales  as  instead  of  /“2  (where  /  is  the  frequency),  thus  more 
power; 

b)  Low  effective  beam  perveance,  thus  higher  efficiency; 

c)  Use  of  permanent  magnets  for  beam  focusing,  thus  lower  energy  consumption; 
and 

d)  Low  cathode  loading,  thus  longer  device  lifetime. 

In  order  to  experimentally  demonstrate  high-power  elliptic  electron  beams  which  can  be 
employed  in  HPM  and  vacuum  electron  devices,  we  must  advance  our  understanding  of 
the  generation,  compression,  focusing  and  transport  of  elliptic  electron  beams. 

Building  upon  our  recent  work  on  elliptic  electron  beam  formation  (Bhatt  and  Chen, 
2005;  Bhatt,  Bemis  and  Chen,  2005  and  2006)  and  periodically  twisting  elliptic  electron 
beams  (Zhou,  Bhatt  and  Chen,  2006),  we  have  developed  a  kinetic  theory  of  periodically 
twisting  elliptic  beams  and  both  cold-fluid  and  kinetic  theories  of  non-twisting  elliptic 
electron  beams  under  the  auspices  of  the  grant  (Chen,  2006). 

4.1  Kinetic  Equilibrium  Theory  of  Periodically  Twisting  Elliptic  Electron  Beams 

(Zhou,  2006;  Zhou  and  Chen,  2006a  and  2006b) 

We  have  developed  our  kinetic  equilibrium  theory  to  examine  effects  of  beam 
temperature  on  periodically  twisting  elliptic  electron  beam.  In  the  kinetic  equilibrium 
theory  (Zhou,  2006;  Zhou  and  Chen,  2006a  and  2006b),  we  have  derived  a  constants  of 
motion  analogous  to  the  Courant-Snyder  invariant  (Courant  and  Snyder,  1958).  We  have 
constructed  a  Vlasov  beam  equilibrium  distribution  of  the  Kapchinskij-Vladimirskij  form 
using  the  two  constants  of  motion.  We  have  obtained  the  generalized  envelope  equations 
which  include  beam  temperature  effects.  In  the  cold-fluid  limit,  the  generalized  envelope 
equations  recover  those  in  the  cold-fluid  beam  equilibrium  theory  (Zhou,  2006;  Zhou, 
Bhatt  and  Chen,  2006).  Detailed  results  are  available  in  our  paper  (Zhou  and  Chen, 
2006b)  and  Zhou’s  doctoral  thesis  (Zhou,  2006). 

To  illustrate  the  effects  of  beam  temperature,  we  have  considered  a  relativistic  elliptic 
beam  with  1^  =  198.5  keV  ,  current  Ih=  85.5  A,  aspect  ratio  a/b  =  5,  and  non- 

axisymmetric  periodic  permanent  magnet  focusing  with  B0  =  2.4  kG  ,  S  =  2.2  cm,  and 
A0v  k0x  =1.52.  Such  a  relativistic  elliptic  beam  could  be  used  in  a  10  MW  L-Band 

ribbon-beam  klystron  (RBK)  for  the  International  Linear  Collider  (ILC). 

As  shown  in  Fig.  4.1.1,  the  solid  curves  represent  the  beam  semi-axis  envelopes  and 
twist  angle  with  zero  temperature  which  is  corresponding  to  a  cold  beam,  while  the 
dotted  curves  represent  the  beam  envelopes  and  twist  angle  with  2.5  keV  on-axis 
temperature.  It  is  evident  in  Fig.  4.1.1  that  the  temperature  effects  on  the  beam  envelopes 
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and  twist  angle  are  negligibly  small.  Since  an  actual  relativistic  elliptic  beam  in  a  well 
designed  system  will  have  a  temperature  which  will  be  considerably  less  than  2.5  keV, 
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Figure  4.1.1  Plots  of  (a)  envelopes  a(s)  and  /)(s )  and  (b)  twist  angle  tf(s)  versus  the 
axial  distance  5  for  the  relativistic  periodically  twisting  elliptic  electron  beam.  The  solid 
curves  are  the  generalized  envelope  solution  for  a  zero-temperature  beam,  whereas  the 
dotted  curves  are  for  a  2.5  keV  on-axis  temperature  beam  (from  Zhou  and  Chen,  2006b). 

the  results  in  Fig.  4.1.1  implies  that  the  temperature  effect  on  the  beam  envelopes  and 
twist  angle  is  expected  to  be  negligibly  small  (Zhou  and  Chen,  2006b). 

4.2  Cold-Fluid  Equilibrium  Theory  of  Non-Twisting  Elliptic  Electron  Beams  (Bhatt, 
2006;  Zhou,  2006;  Chen  and  Zhou,  2007;  Zhou,  Bhatt  and  Chen,  2008) 

For  HPM  and  vacuum  electron  device  applications,  non-twisting  elliptic  electron 
beams  are  desirable.  Under  auspices  of  the  present  grant  (Chen,  2006),  we  have 
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developed  a  cold-fluid  equilibrium  theory  of  a  non-twisting  elliptic  electron  beam  (Bhatt, 
2006;  Chen  and  Zhou,  2007;  Zhou,  Bhatt  and  Chen,  2008),  which  is  a  generalization  of 
our  cold-fluid  equilibrium  theory  of  a  periodically  twisting  elliptic  electron  beam  (Zhou 
2006;  Zhou,  Bhatt  and  Chen,  2006). 

We  have  used  the  combination  of  a  periodic  non-axisymmetric  magnetic  field  and  a 
quadrupole  magnetic  field  to  focus  a  nearly  straight  large-aspect-ratio  elliptic  beam, 
whose  twist  angle  vanishes  approximately.  The  (nonlinear)  cold-fluid  equilibrium 
equations  are: 


Ac— n,+VA.(»,V1)=  0, 
os 

=-4mjn„. 


(4.2.1) 

(4.2.2) 

.  (4.2.3) 

For  the  beam  dimensions  small  relative  to  the  characteristic  scale  of  magnetic  variations, 
'•e->  (koxx)2/6«  • and  [hyyf /6  «  1 ,  the  combined  magnetic  field  can  be  described  to 
the  lowest  order  in  the  transverse  dimensions  as 
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S  is  the  axial  periodicity  length,  and 


We  seek  solutions  to  Eqs.  (2.4.2.  l)-(  2.4.2. 3)  of  the  form 


«a(xx,js)= 


Nh 


-0 


1- 
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mi(s)b(s)  a2  (.9)  b2(s) h 

vi(xi^)=  +  [/q (,v).v +  «vCv)xKce,  .  (4.2.6) 

In  Eqs.  (4.2.5)  and  (4.2.6),  x±  =  xe~  is  a  transverse  displacement  in  the  twisted 

coordinate  system;  is  the  twist  angle  of  the  ellipse;  ®(x)=  1  if  jc  >  0  and  ®(x)~  0 
if  .y<0;  and  the  functions  a(s),  fe(s),  //v(s),  /iv(s),  ax{s),  #v(s)  and  #(s)  are  to  be 

detennined  self-consistently  [see  Eqs.  (4.2.8)-(4.2.12)].  The  self-electric  and  self- 
magnetic  fields  are  well  known  for  an  elliptical  beam  with  density  distribution  specified 
in  Eq.  (4.2.5),  i.e., 

2qNh 


P=Pb-'A‘  = 
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-2  \ 


\a  b  ) 


(4.2.7) 


Using  the  expressions  in  Eqs.  (4.2.4)-(4.2.7),  we  have  shown  that  both  the 
equilibrium  continuity  equation  (4.2.1)  and  force  equation  (4.2.3)  are  satisfied  if  the 
dynamical  variables  0(5),  b(s),  /jx(s)=  a  'da/ ds ,  /jv(s)  =  b~'db/ ds ,  «,(■?),  ax  (•O  and 

#(.9)  obey  the  (cold-fluid)  generalized  beam  envelope  equations  (Chen  and  Zhou,  2007; 
Zhou,  Bhatt  and  Chen,  2008) 
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Equations  (4.2.8)-(4.2.12)  support  non-twisting  elliptic  electron  beam  solutions  in  the 
non-relativistic  regime  which  is  of  interest  to  vacuum  electron  device  applications,  as 
well  as  in  the  relativistic  regimes  which  is  of  interest  to  HPM  applications. 


4.3  Kinetic  Equilibrium  Theory  of  Non-Twisting  Elliptic  Electron  Beams  (Bhatt, 

2006;  Chen  and  Zhou,  2007;  Zhou,  Bhatt  and  Chen,  2008) 

We  have  developed  our  kinetic  equilibrium  theory  to  examine  effects  of  beam 
temperature  on  non-twisting  elliptic  electron  beam.  In  the  kinetic  equilibrium  theory 
(Bhatt,  2006;  Zhou,  Bhatt  and  Chen,  2008),  we  have  derived  a  constants  of  motion 
analogous  to  the  Courant-Snyder  invariant  (Courant  and  Snyder,  1958).  We  have 
constructed  a  Vlasov  beam  equilibrium  distribution  of  the  Kapchinskij-Vladimirskij  form 
using  the  constant  of  motion.  We  have  obtained  the  (kinetic)  generalized  envelope 
equations  which  include  beam  temperature  effects.  In  the  cold-fluid  limit,  the  (kinetic) 
generalized  envelope  equations  recover  the  cold-fluid  ones  (Chen  and  Zhou,  2007;  Zhou, 
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Bhatt  and  Chen,  2008).  Detailed  results  are  available  in  our  paper  (Zhou,  Bhatt  and  Chen, 
2008)  and  Bhatt’s  doctoral  thesis  (Bhatt,  2006). 

4.4  Application  and  Simulation  Validation  of  Cold-Fluid  and  Kinetic  Equilibrium 
Theories  (Bhatt,  2006;  Chen  and  Zhou,  2007;  Zhou,  Bhatt  and  Chen,  2008) 

As  an  example,  we  have  considered  a  relativistic  elliptic  beam  that  can  be  used  in  a 
10  MW  L-Band  ribbon-beam  klystron  (RBK)  for  the  International  Linear  Collider  (1LC). 
The  beam  has  a  current  of  I  b  =  1 1 1 .1  A,  a  voltage  of  Vb  =  120  kV  and  an  aspect  ratio  of 

20:1  .  Solving  the  (cold-fluid)  generalized  envelope  equations  (4.2.8)-(4.2.12),  the 
hybrid  magnetic  fields  are  determined  to  be  the  form  of  Eq.  (4.2.4)  with 
B:(s)  =  -2000sin(/:0.v)  G  ,  B'q  =  80.8  G/cm  ,  S  =  2.2  cm  ,  and  k0v/k0x  =  20  .  In  Fig. 

4.4.1,  the  solid  curves  are  the  beam  semi-axes  a(s)  and  calculated  from  the  (cold- 
fluid)  generalized  envelope  equations,  whereas  dotted  curves  are  from  the  self-consistent 
PIC  PFB2D  simulation.  The  twist  angle  vanishes,  i.e.,  #(s)  =  0  ,  in  this  example. 

To  study  the  temperature  effects  in  the  111.1  A  and  120  kV  elliptic  electron  beam,  we 
have  solved  the  (kinetic)  generalized  envelope  equations  with  nonzero  initial  thermal 

emittances,  i.e.,  x  =  kBTo2(s  =  0)/ myh  fibc~  and  s}h  v  -  kBTb^(s  -  0)y  myb  pbc2  .  As 

shown  in  Fig.  4.4.2,  the  elliptic  beam  envelopes  are  calculated  for  three  different 
temperature  choices:  0  eV,  50  eV  and  100  eV.  Compared  with  the  cold  beam  envelopes 
shown  as  solid  curves  in  Fig.  4.4.1,  the  warm  beam  envelopes  are  found  to  increase 
slightly  as  the  beam  temperature  increases,  while  the  aspect  ratio  of  the  beam  decreases 
from  20:1  to  1 6.4: 1  as  the  beam  temperature  increases  from  0  eV  to  100  eV. 
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Figure  4.4.1  Plots  of  the  beam  envelopes  a(s)  and  6(s)  versus  the  axial  distance  s  for 
the  111  1  A,  120  kV  relativistic  elliptic  electron  beam  with  zero  temperature.  The  solid 
curves  are  the  (cold-fluid)  generalized  envelope  solution,  whereas  the  dotted  curves  are 
from  the  PFB2D  simulation  (from  Zhou,  Bhatt  and  Chen,  2008). 
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Figure  4.4.2  Plots  of  the  beam  envelopes  (a)  and  (b)  b(s)  versus  the  axial  distance 
s  for  the  111.1  A,  120  kV  relativistic  elliptic  electron  beam  for  three  different 
temperature  choices:  0  eV,  50  eV  and  100  eV  (from  Zhou,  Bhatt  and  Chen,  2008). 

For  a  practical  HPM  device,  we  have  assumed  that  the  elliptic  beam  is  generated  from 
an  electron  gun  with  an  intrinsic  temperature  of  0.1  eV  and  a  current  density  of  1.5  A/cm. 
The  elliptic  beam  has  to  be  compressed  by  a  factor  of  471.5  in  area  to  achieve  a  current 
density  of  707.3  A/cm'  in  the  focusing  channel.  During  the  compression,  the  temperature 
increases  by  a  factor  of  471.5  to  47.2  eV.  Therefore,  in  our  calculations,  a  temperature  of 
50  eV  is  a  reasonable  assumption.  As  shown  in  Fig.  4.4.3,  the  envelopes  of  the  elliptic 
beam  with  a  temperature  of  50  eV  are  obtained  by  solving  the  (kinetic)  generalized 
envelope  equations  (solid  curves)  and  by  the  PFB2D  PIC  simulations  (dotted  curves). 
Both  results  showed  a  slight  increase  in  the  two  envelope  dimensions  and  the  aspect  ratio 
of  the  elliptic  beam  decreases  to  17.8,  compared  with  the  cold  elliptic  beam. 
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Figure  4.4.3  Plots  of  the  beam  envelopes  and  />(s)  versus  the  axial  distance  5  for 

the  1 1 1 . 1  A,  120  kV  relativistic  elliptic  electron  beam  for  a  temeprature  of  50  eV.  The 
solid  curves  are  the  (kinetic)  generalized  envelope  solution,  whereas  the  dotted  curves  are 
from  the  PFB2D  simulation  (from  Zhou,  Bhatt  and  Chen,  2008). 
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Kinetic  equilibrium  of  a  periodically  twisted  ellipse-shaped  charged-particle  beam 
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A  Vlasov  equilibrium  of  the  Kapchinskij-Vladimirskij  form  is  obtained  for  a  periodically  twisted 
ellipse-shaped  charged-particle  beam  in  a  nonaxi symmetric  periodic  magnetic  focusing  field  The  single¬ 
particle  Hamiltonian  dynamics  is  analyzed  self-consistently.  A  constant  of  motion  analogous  to  the 
Courant-Snyder  invariant  is  found  The  equilibrium  distribution  function  is  constructed  The  statistical 
properties  of  the  beam  equilibrium  are  studied.  In  the  zero-temperature  limit,  the  generalized  envelope 
equations  derived  from  the  kinetic  equilibrium  theory  recover  the  generalized  envelope  equations  obtained 
in  the  cold-fluid  equilibrium  theory.  Examples  of  periodically  twisted  elliptic  beam  equilibria  are 
presented,  and  potential  applications  are  discussed.  For  ribbon-beam  amplifier  and  ribbon-beam  klystron 
applications,  the  kinetic  equilibrium  theory  predicts  that  the  effect  of  beam  temperature  on  the  beam 


envelopes  is  negligibly  small. 
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l.  INTRODUCTION 

A  fundamental  understanding  of  the  kinetic  equilibrium 
and  stability  properties  of  high-intensity  electron  and  ion 
beams  in  periodic  focusing  fields  is  desirable  in  the  design 
and  operation  of  particle  accelerators  [1-14],  such  as  stor¬ 
age  rings  and  rf  and  induction  linacs,  as  well  as  vacuum 
electron  devices,  such  as  klystrons  and  traveling- wave 
tubes  with  periodic  permanent  magnet  (PPM)  focusing. 
There  are  two  well-known  equilibria  for  periodically 
focused  intense  beams,  including  the  Kapchinskij- 
Vladmirskij  (KV)  equilibrium  [7-9]  in  an  alternating- 
gradient  quadrupolc  magnetic  focusing  field  and  the  peri¬ 
odically  focused  rigid-rotor  Vlasov  equilibrium  [10,1 1]  in 
a  periodic  solenoid  magnetic  focusing  field.  In  general,  for 
linear  focusing  forces,  self-consistent  beam  distributions 
can  be  formally  constructed  using  a  matrix  formulation 
[12.13]. 

It  was  shown  formally  in  Ref.  [12]  that  self-consistent 
beam  distributions  can  be  obtained  that  allow  elliptical 
space-charge  beams  of  arbitrary  aspect  ratio  and  with 
arbitrary  rotation  angle  of  the  ellipse  as  long  as  the  field 
is  linear.  However,  obtaining  concrete  equilibria  with  non¬ 
upright  ellipses  is  nontrivial  The  previous  explicitly 
known  Vlasov  equilibria  of  KV  form  [7,9-1  1]  for  high- 
mtensity,  space-chargc-dominatcd  charged-particle  beams 
propagating  in  the  alternating-gradient  quadrupole  mag¬ 
netic  focusing  field  and  the  periodic  solenoid  magnetic 
focusing  field  charged-particle  beams  are  circular  on  aver¬ 
age  that  is,  the  averages  of  the  beam  envelopes  in  different 
transverse  directions  over  one  period  are  the  same. 

There  is  considerable  interest  in  the  research  and  devel¬ 
opment  of  high-intensity  charged-particle  beams  with  a 
large  aspect  ratio  transverse  to  the  direction  of  propagation. 
First,  large-aspect-ratio  elliptic  beams  (or  ribbon  beams) 
can  transport  larger  amounts  of  beam  currents  at  reduced 


PACS  numbers:  29.27  -a.  52.59.Sa,  41.75.-i,  52.25.Dg 

intrinsic  space-charge  forces  and  energies.  Second,  they 
couple  efficiently  to  planar  or  rectangular  rf  structures.  The 
combination  of  the  space-charge  reduction  and  efficient 
coupling  allows  efficient  rf  generation  in  vacuum  elec¬ 
tronic  devices,  and  efficient  acceleration  in  particle  accel¬ 
erators.  Third,  elliptic  beams  provide  an  additional 
adjustable  parameter  (c.g.,  the  aspect  ratio)  which  may 
be  useful  for  better  matching  a  beam  into  a  periodic 
focusing  channel  [14]. 

One  important  application  of  ribbon  beams  is  in  the 
development  of  advanced  radiation  devices  such  as 
ribbon-beam  amplifiers  (RBAs)  and  ribbon-beam  klystrons 
(RBKs),  which  have  advantages  over  the  corresponding 
conventional  (circular-beam)  devices  in  terms  of  efficiency 
and  operational  parameters.  Other  applications  include  the 
development  of  advanced  accelerators  capable  of  generat¬ 
ing  nonconventional  beams,  e.g.,  a  planar  radio-frequency 
(rf)  hnac  producing  ribbonlike  bunches  of  charged 
particles. 

Although  sheet  beams  have  been  discussed  in  the  litera¬ 
ture  for  four  decades,  the  Vlasov  equilibrium  of  a  high- 
intensity,  space-charge-dominated  beam  with  a  large-as- 
pcct-ratio  elliptic  cross  section  has  not  been  discovered 
until  this  paper.  Sturrock  [15]  first  suggested  use  of  a 
periodic  magnetic  focusing  consisting  of  an  array  of 
planar-wiggler  magnets  for  rectilinear  beams.  Zhang  et 
al  [16]  had  some  modest  success  in  the  experimental 
demonstration  of  the  transport  of  a  low-intensity  (10  A, 
500  kV)  sheet  beam  in  a  planar-wiggler  magnetic  field,  and 
observed  considerable  beam  loss.  Researchers  made  use  of 
the  multiple-timc-scalc  analysis  and  the  paraxial  approx¬ 
imations  to  obtain  the  smooth-beam  approximation  of 
high-intensity  sheet  beam  equilibria  [17,18].  Recently, 
Russell  et  al  demonstrated  the  transformation  of  a  circular 
beam  into  a  sheet  beam  using  asymmetric  lenses  [  19].  The 
authors  discovered  a  cold-fluid  equilibrium  for  a  high- 
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intensity  periodically  twisted  elliptic  beam  in  a  nonasym- 
metric  periodic  magnetic  focusing  field  [6]. 

In  this  paper,  it  is  shown  that  there  exists  a  Vlasov 
equilibrium  for  a  periodically  twisted  large-aspect-ratio 
intense  charged-particle  beam  with  uniform  density  in 
the  transverse  plane  propagating  through  a  nonaxisymmet- 
ric  periodic  magnetic  focusing  field.  The  single-particle 
Hamiltonian  of  such  a  periodically  twisted  large-aspect- 
ratio  elliptic  beam  is  investigated.  The  constant  of  motion 
analogous  to  the  Courant-Snyder  invariant  [20]  is  found. 
The  equilibrium  beam  distribution  is  constructed.  The 
beam  envelope  equations  and  flow  velocity  equations  are 
derived.  In  the  zero-temperature  limit,  they  are  consistent 
with  the  generalized  envelope  equations  derived  from  cold- 
fluid  equilibrium  theory  [6].  Statistical  properties  and  pos¬ 
sible  applications  of  the  present  beam  equilibrium  are 
discussed. 


y 


II.  VLASOV  EQUILIBRIUM  THEORY 


We  consider  an  ellipse-shaped,  continuous,  intense 
charged-particle  beam  propagating  with  constant  axial 
velocity  /3hcez  through  an  applied  nonaxisymmetric  peri¬ 
odic  magnetic  focusing  field.  The  applied  nonaxisymmet¬ 
ric  periodic  magnetic  focusing  field  inside  the  thin  beam 
can  be  approximated  by 


r  k 2  k2 

Bext  =  —  BQs\n(kQs)e.  4  B0cos(fc0.s)  -~xev  -I- 

L  *o  *  o 


0) 


where  s  =  z  is  the  axial  coordinate,  /c0  =  27 r/S,  k%x  + 
/cqv  =  and  S  is  the  axial  periodicity  length.  The  3D 
magnetic  field  in  Exp  (4)  is  fully  specified  by  the  three 
parameters  B{h  S ,  and  k0x/k 0v.  The  associated  magnetic 
vector  potential  can  be  expressed  as 


—  B0  sin(fc0s)r  —  +  -jrxey 

L  *o  *o 


(2) 


To  determine  the  sclf-electric  and  self-magnetic  fields  of 
the  beam  sclf-consistently  in  the  present  paraxial  approxi¬ 
mation,  we  assume  that  the  density  profile  of  the  beam  is 
uniform  inside  the  beam  boundary,  i.c.. 


nh(x,  y) 


Nh 

t ra(s)b(s) 


b2(s)f 


(3) 


In  Eq.  (3),  x  —  acos[0(.v)]  4  y  sin[0(jr)]  and  y  = 
— xsin[0(,v)]  4  ycos[0(.v)]  represent  the  twisted  coordinate 
as  illustrated  in  Fig.  1 ;  0(s)  is  the  twist  angle  of  the  ellipse; 
0(a)  =  I  if  x  >  0  and  0(a)  =  0  if  x  <  0.  The  density  of 
the  elliptic  beam  with  semimajor  axis  a(jr)  =  a(s  T  S) 
and  semiminor  axis  b(s)  =  b{s  4-  S)  is  uniform  in  the 
beam  interior  (x2/cr  T  yr/b 2  ^  1).  The  semimajor 
and  semiminor  axes  have  the  same  periodicity  as  the 


applied  magnetic  field  with  S  =  2i r/k^.  Nb  = 
2tt  j°°  dxdynh(x,  y ,  s)  =  const  is  the  number  of  particles 
per  unit  axial  length  In  the  paraxial  approximation,  the 
Budker  parameter  of  the  beam  is  assumed  to  be  small,  i.c., 
q2Nb/mc 2  «  yb ,  and  the  transverse  kinetic  energy  of  a 
beam  particle  is  assumed  to  be  small  compared  with  its 
axial  kinetic  energy.  Here,  c  is  the  speed  of  light  in  vacuo, 
yh  =  (1  —  is  the  relativistic  mass  factor,  and  q  and 

m  arc  the  particle  charge  and  rest  mass,  respectively. 

From  the  equilibrium  Maxwell  equations,  we  tind  that 
the  sclf-electric  and  self-magnetic  fields,  Ef  and  B\  are 
given  by 


E  s{x,  y,  s) 


2qNb 


n(.v)  4  b(s)  Ln(.s') 


(4) 


2  qNh 

tf(jr)  4-  b{s) 


4 


(5) 


in  the  beam  interior  (a 2 /a2  4  yr/b 2  ^  1).  It  is  convenient 
to  express  the  self-fields  in  terms  of  the  scalar  and  vector 
potentials  defined  for  xr/a2  4  y2/b~  <  1  by 


9>Ha,  y,  s)  =  ft hlAt( a,  y,  s) 


2gNb 

a(s)  4  b(s) 


_ 

A** 


_ 
V“ 


_a(s)  b(s)_ 


(6) 


where  A* (a,  v,  s)  =  A!(a,  y,  s)er,  E*(x,  v,  .v)  = 
—  (Cfd/dx  4  Cjd/dy)<t>x,  and  Bv  =  (— c+d/dx  4 
exd/dy)AL 

In  the  paraxial  approximation,  the  transverse  motion  for 
an  individual  particle  in  the  combined  self-fields  and  ap¬ 
plied  magnetic  field  is  described  by  the  normalized  per¬ 
pendicular  Hamiltonian  H L  =  H ±/ybPbmc2, 
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where  (jc,  Pv)  and  (y,  Py)  arc  canonical  conjugate  pairs, 

=  Q Bzti)/2yb/3bntc 2,  K  =  2  q1Nb/yl02bmc2  is 
the  self-field  perveance,  and  the  normalized  transverse 
canonical  momentum  Px  —  ( Px ,  Py)  is  related  to  the  trans¬ 
verse  mechanical  momentum  px  by  Px  =  ( yh/3hmc)~l  X 
(Pl  +  q\el'/c). 

It  is  convenient  to  transform  the  Hamiltonian  in  the 
Cartesian  canonical  coordinate  (jc,  y,  Px,  Py)  to  a  new  ca¬ 
nonical  coordinate  (jcj,  y\,  PrI,  Py\),  so  that  the  new 
Hamiltonian  assumes  a  simpler  form  from  which  the  in¬ 
variants  of  motion  are  easily  identified.  The  transformation 

ot  the  Hamiltonian  from  the  Cartesian  canonical  coordi- 
_ I 


nate  (jc,  y,  Px ,  Py)  to  the  canonical  coordinate 

(jcj ,  >  | ,  Px\,  Pyi)  involves  two  steps. 

We  first  transform  the  Cartesian  canonical  coordinate 
(jc,  y,  Px,  Py)  to  a  twisted  canonical  coordinate  (.*,  y,  P,,  Pv) 
using  the  generating  function 

F2(x,  y;  Px,  Py,  s)  =  Px[xcos6{s)  4*  y  sintf(.v)] 

+  Pv[— jcsin0($)  4-  y  cos0(a*)].  (8) 

It  follows  from  Eq.  (8)  that 


Px  —  Px  cos#(.v)  —  Py  sin0(.v),  (9a) 

Py  =  Pxs'\n0(s)  4  Pvcos0(,v),  (9b) 

Jc  =  *cos0(.y)  4-  ysin0(.y),  (9c) 

y  =  —  jcsin0(j»)  4-  ycosd(s).  (9d) 


The  Hamiltonian  in  the  twisted  canonical  coordinate  is 
then  expressed  as 


H±(x,  y,  Px,  Pyt  s)  =  H±(x,  y,  Px,  Py,  s)  4-  dF2/ds 

1 


_ j-2  _ f.  2 

Pxcos6(s)  —  PySinO(s)  4-  Jkz(s)-^x  sin0(y)  4-  J kz(s)  y  cos 0(s) 

*0  "  N) 


4-  ■ 


1 


PY  sin 


I - kr  i - k~ 

in0(j)  4-  Pvcos0(j»)  —  ^Kz(s)-^xcosO(s)  4-  yj kz(s)  y  ^inO(s) 


kL 


K  (x1  y\  -  ^dd(s) 

\-  +  T)  +  (pxy-P>V- 


a  +  b\a  b  J  ds 

The  equations  of  motion  associated  with  the  Hamiltonian  in  Eq.  (10)  are 

•*'  =  =  Px  +  C(s)x  4-  [0'(jr)  -  <*,(*)]>. 

o  r  x 

y'  =  =  Py  -  C{s)y  ~  [0'(jt)  -  a,(i)]i, 

<>P> 

_ dH± _ j  j&cosHg(j)]  +  <.-Qysin2[6l(5)]  _  2 K 

dx  {  /:q  4-  b(s)] 

!'4  _  i.4 


*0v  “  ko.x\  sin[20(.y)]  _ 

~ — y> 


(  \(K 0y  K0.x\ 


2  K 


_  BHi  _  f  ^o.tsiir[g(j)]  +  ^.cos^flU)] _ 

dy  [K:  S  Id  fe(.v)[a(.v)  +  /;(.?)] 


x  -  C(s)Px  +  [0'(s)  -  afv(j)]Pv 


v  +  C(s)Py  -  [0'(.v)  -  ar(.v)]P, 


K.(s) 


koy  ~  kox\  sin[20(.r)]  . 


JC, 


*S  J  2 

where  prime  denotes  derivative  with  respect  to  s. 


(10) 


(Ha) 

(lib) 


(He) 


(lid) 


i - k2  -  k2 

C{s)  =  ^[i<Js)~- - — 


Ik2 


sin[20(.r)], 


(12) 
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«,(•'■)  =  -yjKz(s) 
and 


^ysin2[0(.?)]  +  -yeos2[0(.v)] 
*5  *o 


(13) 


r)  F 

Px\’Px\.s)  =  Wi(.v.  y,  Pt.  Pv,  j)  +  —=■ 

as 


p2*'  +  ZiL  +  VT 


*’?($)  wl(s)  u;(.v) 


aY(s)  = 


7#cz(o 

^eos2[0(s)]  +  ^sin2[0(s)] 

■  (14) 

1 

4- 

V 

*o  *5  J 

WyM. 

The  functions  ar  and  ay  are  related  to  the  variables  ax  and 
ay  in  the  eold-fluid  equilibrium  theory  [6].  Indeed,  by 
adding  Eqs.  (10)  and  (11)  and  subtracting  Eq.  (10)  from 
Eq.  ( 1 1 )  in  [6],  and  carrying  out  the  integrations  on  result¬ 
ing  equations  with  the  initial  condition  av(0)  =  a'y(O)  = 
0,  it  is  readily  shown  that  the  functions  ax  and  ay  in 
Eqs.  (13)  and  (14)  correspond  to  a  particular  solution  to 
Eqs.  (10)  and  (11)  in  [6], 

As  a  second  step,  we  apply  another  transformation  from 
the  twisted  canonical  variables  (Jc,  y,  Px,  Py)  to  the  canoni¬ 
cal  variables  (A|,  y]t  Px Py\)  using  the  generating  function 


F2(x,y,Pxl,P,t,s)  = 


■<(s) 

.Wx(s) 


-  c(.v) 


jf2 


1 

+  — 

2 


I'vW 


+ 


+  C(s) 
»VtU)  »»\.(j)’ 


XPr 


■  + 


(15) 


where  h\(j)  =  w  v(.r  +  S )  and  wv(s)  =  wv(.y  +  S)  are  pc- 
riodic  functions  solving  the  differential  equations 


<(*) 

M^j) 


-  C'(s)  -  C2(s)  -  a v  - 


2  K 


a($)[a($)  +  b(s)] 


U’4t(i)  ’ 


(16) 


"\(.v) 


+  C'(.r)  -  C2(*)  -  ax  - 


2  K 


b(s)[a(s)  +  b(s)] 


It  follows  from  Eq.  (15)  that 
rw'x(s) 


p*  = 

Py  = 

A  |  = 

Vi  = 


W.x(x) 

>y(-0 

H'v(.t) 

A 


-  c{s) 

+  C(s) 


Jf  + 

y  + 


vl 


w\($)  ’ 

Pyl 

w,(s)’ 


wx(s)’ 

y 

w’v(.r) ' 


(17) 

(18a) 

(18b) 

(18c) 

(18d) 


The  Hamiltonian  in  the  canonical  coordinate 
(r i ,  V|,  Px[,  Px])  is  then  expressed  as 


+  ^ ,  _  )  (19) 
as 


where  we  have  introduced  and  demanded 


d(f(s)  _  WyW 


ds 


d0(s) 


ds 


"',(■?) 

Wy(l) 


d0(s) 


c Is 


(20) 


Following  Eq.  (20),  it  can  be  shown  that  the  twisted  angle 
0(.y)  has  to  satisfy  the  differential  equation 


dd(s)  _  Wt(s)«v(j)  -  W’yfcW-v) 
w2(.r)  -  u’2.(.v) 


ds 


(21) 


The  motion  described  by  the  simplified  Hamiltonian  in  the 
new  canonical  coordinate  (x\,  V|,  Px\,  Py\)  in  Eq.  (19)  is 
described  by  the  equations 


dH 


ii 


Px 


9P,i 

dH  u 


,  d<p(s) 

+  - : - Vi 


\v^(i)  ds 
Py  I  dif(s) 


dP 


Kx  =  - 


Pit 


dH 


H'vM 


ds 


-Vi. 


U 


dH  |i 


w-(s) 

y  i 


■  + 


d(f(s) 

ds 


(22a) 

(22b) 

Pyl. 

(22c) 

Px, • 

(22d) 

(23) 

wj(.f)  ds 
From  Eq.  (22),  it  is  readily  shown  that 

E  -jrf +.vf +  PJ,  +  #*, 


is  an  exact  single-particle  constant  of  the  motion  for  the 
Hamiltonian  in  Eq.  (19). 

We  consider  the  following  trial  choice  of  the  Vlasov 
equilibrium  distribution  function: 


fh(Xt’yuPx\’Px\’S)=-T- 

TTEj 


h  S(x]+y2l+P2xl+P;,-sT). 

(24) 


where  dfb/ds  =  0,  er  =  const  >  0  is  an  effective  emit- 
tance,  and  5(a)  is  the  Dirac  8  function.  As  will  be  shown  in 
Sec.  Ill,  the  density  profile  of  the  beam  described  by  the 
distribution  function  fb  is  consistent  with  the  uniform- 
density  profile  within  an  ellipse,  which  is  the  key  require¬ 
ment  for  the  quantity  E  =  a*  4-  y]  +  P;,  +  Pjj  to  be  a 
constant  of  motion.  Therefore,  the  distribution  function 
defined  in  Eq.  (24)  is  indeed  a  Vlasov  equilibrium,  i.c.. 
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Hi  =  0.  (25) 

ds 


111.  STATISTICAL  PROPERTIES 


The  distribution  function  described  in  Eq.  (24)  has  the 
following  statistical  properties.  First,  the  distribution  is 
consistent  with  the  assumed  density  profile  in  Eq.  (3),  i.e.. 

nb(.l  y,S)=J—  ff  fdPxldPyl 

=  Nh/{i TeTwxwyX  if  x\/eT  +  y\/eT  <  1. 
0,  otherwise. 

(26) 


The  beam  has  the  uniform-density  profile  given  in  Eq.  (3), 
provided  that  a  —  and  b  =  yfsrwx  are  satisfied. 

Under  these  self-consistent  conditions,  Eqs.  (16),  (17),  and 
(21)  can  be  expressed  as 


d2a  f  b2(a2x  —  2axay)  +  a2  a2. 


d-a  _  f 
~d? 


2  i  2 

a  —  b 


l<2  _  k2 

+  — x— — —  cos(A’o*v)  sin(20)  -2^K^aysm{k^s) 


k0 


2  K 


a  4-  b  a 3 

d2b  f  a2(a 2  —  2axay)  +  b2a2 


(27) 


d~b  fc 
ds2  + 


1  i 1 

a  —  b~ 


k2  -  k 2 

*0.t  K0v 


cos(/cos)sin(20)  +  2jK~^axsin(k0s) 


2  K  e2 


a  +  b  b5 


d6  croty  -  b2  a  x 

ds  cr  -  b2 


(28) 

(29) 


Equations  (27)— (29)  are  written  in  a  form  similar  to  the 
generalized  envelope  equations  in  the  cold-fluid  equilib¬ 
rium  theory  [6].  They  are  identical  to  the  generalized 
envelope  equations  of  a(s),  b(s),  and  6(s)  in  the  cold-fluid 
equilibrium  theory,  except  that  the  emittanee  terms  appear¬ 
ing  on  the  right-hand  side  of  Eqs.  (27)  and  (28)  are  zero  in 
the  cold-fluid  equilibrium  theory.  Therefore,  they  are  more 
general  than  the  eold-fluid  equilibrium  theory. 

Second,  in  the  normalized  units,  the  average  (macro¬ 
scopic  flow)  transverse  velocity  of  the  beam  equilibrium 
described  by  Eq.  (24)  is  given  in  the  tw  isted  coordinates  by 


Vi  =  (  -  f  Vi  fdPxldPtX 

\7 TETWXWJ  WXWY  J 


(30) 


The  flow  velocity  in  Eq.  (30)  is  identical  to  the  flow 
velocity  derived  by  the  cold-fluid  equilibrium  theory  [6] 
provided  that  the  relations  fix  =  a* /a  =  w'/w\  and  /xy  = 
bf/b  =  Wy/\vy  are  satisfied. 

As  a  third  statistical  property,  the  beam  equilibrium 
described  by  Eq.  (24)  has  the  effective  transverse  tempera¬ 
ture  profile  (in  dimensional  units) 


Ti(x,  y,  s) 


Vi  )2fdPxldPyl 


2 


(31) 


As  the  fourth  property,  the  rms  emittances  of  the  beam  in 
the  x  and  the  y  directions  are 


F-Ihx  =  ■  (32a) 

Pcc  v  4 

e,h9  =  =  <32b) 


Finally,  the  Vlasov  elliptic  beam  equilibrium  has  two 
limiting  eases  which  are  well  known  It  recovers  the  famil¬ 
iar  periodic  (circular)  rigid-rotor  Vlasov  equilibrium  [10] 
by  setting  the  major-axis  equal  to  the  minor  axis  of  the 
beam  ellipse.  It  also  recovers  the  familiar  constant-radius, 
uniform-density  rigid-rotor  Valsov  equilibrium  [9]  by  tak¬ 
ing  the  limit  of  a  uniform  magnetic  field  with  B:  =  B  = 
const. 


IV.  EXAMPLES 

We  illustrate  examples  of  periodically  twisted  Vlasov 
elliptic  beam  equilibria  in  a  periodic  non  ax  i  symmetric 
magnetic  foeusing  field  and  the  temperature  effects  with 
numerical  calculations.  A  numerical  module  in  the  PFB2D 
code  [6,21]  has  been  developed  to  solve  the  generalized 
envelope  equations  (13),  (14),  and  (27)-(29),  which  deter¬ 
mines  the  rotational  flow  velocity,  the  outer  equilibrium 
major  axis  </($)  and  minor  axis  b(s)  of  the  beam  ellipse,  and 
the  twisted  angle  0(s). 

In  particular,  we  consider  a  nonrelativistie  elliptic  beam 
with  voltage  Vb  =  2.29  keV,  current  1b  =  0. 1 1  A,  aspect 
ratio  a/b  =  6,  and  nonaxisymmetric  periodic  permanent 
magnet  focusing  with  B0  =  337.5  G,  5  =  1.912  cm,  and 
koy/k0x  =  1.6,  which  is  corresponding  to  a  beam  design  for 
a  high-efficiency  200  W  RBA  under  development  at 
Massachusetts  Institute  of  Technology  (MIT)  and  beam 
power  technology  for  wireless  communication  For  such 
a  system  the  matched  solution  of  the  generalized  envelope 
equations  (13),  (14),  and  (27)-(29)  is  calculated  numeri¬ 
cally  as  shown  in  Figs.  2  and  3  for  several  maximum  (on- 
axis)  temperature  choices  and  fixed  parameters:  k()x  = 
1.90  cm-1,  k0y  =  3.03  em"1,  =  1.04  citT1,  and 

K  =  1.52  X  10-2.  The  solutions  to  the  generalized  enve- 
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FIG.  2.  Plots  of  (a)  envelopes  c/(.v)  and  b(s)  and  (b)  twist  angle 
tf{s)  versus  the  axial  distance  s  for  the  nonrelativistic  twisted 
elliptic  beam  The  solid  curves  are  the  generalized  envelope 
solution  for  a  zero-temperature  beam,  whereas  the  dotted  curves 
are  for  a  1  eV  on-axis  temperature  beam 

lope  equations  ( 1 3),  ( 1 4),  and  (27)-(29),  displayed  as  solid 
and  dotted  curves  in  Figs.  2  and  3,  show  that  the  semiaxes 
of  the  elliptical  beam  remain  almost  constant  with  small- 
amplitude  oscillations,  that  the  orientation  of  the  ellipse 
twists  periodically  with  an  amplitude  of  ten  degrees. 

The  solid  lines  in  Figs.  2  and  3  represent  the  beam 
envelopes  and  twisted  angle  with  zero  temperature  which 
is  corresponding  to  a  cold  beam,  while  the  dotted  curves 
represent  the  beam  envelopes  and  twisted  angles  with  1  eV 
on-axis  temperature  in  Fig.  2  and  10  eV  on-axis  tempera¬ 
ture  in  Fig.  3,  respectively.  The  aspect  ratio  of  the  beam 
reduces  from  6  to  4  as  the  on-axis  temperature  of  the  beam 
increases  from  0  to  10  eV,  i.e.,  the  elliptic  beam  becomes 
more  circular.  However,  the  twisted  angle  is  almost  un¬ 
changed  as  the  on-axis  temperature  increases  from  0  to 
10  eV.  For  the  elliptic  beam  designed  for  the  200  W  ribbon- 
beam  amplifier,  the  temperature  of  the  beam  is  estimated  to 


FIG.  3.  Plots  of  (a)  envelopes  <7(5)  and  b(s)  and  (b)  twist  angle 
0(.y)  versus  the  axial  distance  s  for  the  nonrelativistic  twisted 
elliptic  beam  The  solid  curves  are  the  generalized  envelope 
solution  for  a  zero-temperature  beam,  whereas  the  dotted  curves 
are  for  a  10  eV  on-axis  temperature  beam. 

be  0.1  eV  from  simulations  [22,23].  In  such  a  case,  the 
temperature  effect  is  negligible. 

To  further  illustrate  the  effects  of  beam  temperature,  we 
consider  a  relativistic  elliptic  beam  with  Vh  =  198.5  keV, 
current  Ib  =  85.5  A,  aspect  ratio  a/b  =  5,  and  nonaxi- 
symmetric  periodic  permanent  magnet  focusing  with  = 
2.4  kG,  S  =  2.2  cm,  and  &oyAo*  =  152.  [Such  a  relativ¬ 
istic  elliptic  beam  could  be  used  in  a  10  MW  /.-band  RBK 
for  the  International  Linear  Collider  (1LC).]  For  such  a 
system  the  matched  solution  of  the  generalized  envelope 
equations  (13),  (14),  and  (27)— (29)  is  calculated  numeri¬ 
cally  with  the  corresponding  parameters:  k{)x  = 
I  57  cm-1,  =  0.732  cm-1,  and  K  =  1.13  X  lO-2. 

As  shown  in  Fig.  4,  the  solid  lines  represent  the  beam 
semiaxis  envelopes  and  twist  angle  with  zero  temperature 
which  is  corresponding  to  a  cold  beam,  while  the  dotted 
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s  (cm) 


s  (cm) 

FIG.  4.  Plots  of  (a)  envelopes  a(s)  and  b(s)  and  (b)  twist  angle 
0(s)  versus  the  axial  distance  s  for  the  relativistic  twisted  elliptic 
beam.  The  solid  curves  are  the  generalized  envelope  solution  for 
a  zero-temperature  beam,  whereas  the  dotted  curves  are  for  a 
2.5  keV  on-axis  temperature  beam. 

curves  represent  the  beam  envelopes  and  twist  angle  with 
2.5  keV  on-axis  temperature.  It  is  evident  in  Fig.  4  that  the 
temperature  effects  on  the  beam  envelopes  and  twist  angle 
are  negligibly  small.  Since  an  actual  relati vistie  elliptic 
beam  in  a  well  designed  system  will  have  a  temperature 
whieh  will  be  considerably  less  than  2.5  keV,  the  results  in 
Fig.  4  imply  that  the  temperature  effect  on  the  beam 
envelopes  and  twist  angle  is  expected  to  be  negligibly 
small. 

V.  CONCLUSIONS 

The  single-particle  Hamiltonian  of  a  periodically 
twisted  large-aspect-ratio  elliptic  beam  in  a  nonaxisym- 
metric  periodic  magnetic  foeusing  field  has  been  investi¬ 
gated.  A  constant  of  motion  analogous  to  the  Courant- 
Snyder  invariant  has  been  found  such  that  the  self- 


consistent  beam  equilibrium  can  be  constructed  as  a  func¬ 
tion  of  the  constant  of  motion.  The  beam  envelope  equa¬ 
tions  and  flow  velocity  equations  have  been  derived.  They 
are  consistent  with  the  generalized  envelope  equations 
derived  from  the  cold-fluid  equilibrium  theory  [6J  when 
the  temperature  is  taken  to  be  zero.  Statistical  properties  of 
the  present  Vlasov  elliptic  beam  equilibrium  have  been 
studied.  For  current  applications  of  interest,  namely,  the 
RBA  and  RBK,  the  temperature  effeets  have  been  found  to 
be  negligibly  small  on  periodically  twisted  large-aspect- 
ratio  elliptic  beams. 

ACKNOWLEDGMENTS 

This  research  was  supported  by  the  U.S.  Department  of 
Energy,  Office  of  High-Energy  Physics,  Grant  No.  DE- 
FG02-95ER40919,  Office  of  Fusion  Energy  Science, 
Grant  No.  DE-FG02-01ER54662,  and  in  part  by  Air 
Force  Office  of  Scientific  Research,  Grant  No.  FA9550- 
06-1-0269. 


[1]  R.C.  Davidson  and  C.  Chen,  Part.  Accel.  59.  175  (1998); 
Nucl.  Instrum  Methods  Phys.  Res.,  Sect  A  415.  370 
(1998). 

[2]  R.C,  Davidson,  P.  Stoltz,  and  C.  Chen,  Phys.  Plasmas  4, 
3710  (1997). 

[3]  R.  Pakter  and  C.  Chen,  Phys.  Rev.  E  62,  2789  (2000). 

[4]  C.  Chen  and  R.  Pakter,  Phys.  Plasmas  7,  2203  (2000). 

[5]  S.  Bernal  et  al ,  Proceedings  of  the  2005  Particle 
Accelerator  Conference,  2005,  p.  892. 

[6]  J.  Zhou,  R.  Bhatt,  and  C.  Chen,  Phys.  Rev  ST  Accel. 
Beams  9,  034401  (2006). 

[7]  I  M.  Kapchinskij  and  V.  V  Vladimirskij,  in  Proceedings 
of  the  International  Conference  on  High  Energy 
Accelerators  (CERN,  Geneva,  1959),  p.  274. 

[8]  S.M.  Lund  and  B  Bukh,  Phys.  Rev.  ST  Accel  Beams  7, 
024801  (2004). 

[9]  R.C.  Davidson,  Physics  of  Nonneutral  Plasma  (Addison- 
Wesley,  Reading,  MA,  1990). 

[10]  C.  Chen,  R.  Pakter,  and  R.C.  Davidson,  Phys.  Rev.  Lett 
79,  225  (1997). 

[11]  R.C.  Davidson  and  H.  Qin,  Physics  of  Intense  Charged 
Particle  Beams  in  High  Energy  Accelerators  (World 
Scientific,  Singapore,  2001),  p.  242. 

[12]  F.  J  Sacherer,  Ph.D.  thesis,  Lawrence  Radiation 
Laboratory,  University  of  California,  Berkeley,  1968, 
pp.  15-19. 

[13]  V.  Danilov.  S  Cousineau,  S  Henderson,  and  J  Holmes, 
Phys  Rev.  ST  Accel.  Beams  6,  094202  (2003). 

[14]  C.  Chen,  R  Pakter,  and  R.C.  Davidson,  Nucl.  Instrum. 
Methods  Phys.  Res.,  Sect.  A  464,  518  (2001). 

[15]  P.  A.  Sturrock,  J  Electron.  Control  7,  162  (1959). 

[16]  Z.  X.  Zhang  et  al .,  IEEE  Trans.  Plasma  Sci.  21,  760 
(1993). 

[17]  J.H.  Booske,  A  H.  Kumbasar,  and  M.  A.  Basten,  Phys. 
Rev.  Lett.  71,  3979  (1993). 


104201-7 


JING  ZHOU  AND  CHIPING  CHEN 


Phys.  Rev.  ST  Accel.  Beams  9,  104201  (2006) 


[18]  M.A.  Basten  and  J  H  Booske,  J.  Appl  Phys.  85,  6313 
(1999). 

[191  S.  J.  Russell  ei  «/.,  Phys.  Rev.  ST  Accel.  Beams  8,  080401 
(2005). 

[20]  E.  D.  Courant  and  H.S  Snyder,  Ann.  Phys.  (N.Y.)  3,  1 
(1958). 


[21]  J.  Zhou,  PhD.  thesis,  Massachusetts  Institute  of 
Technology,  2006. 

[22]  R  Bhatt  and  C.  Chen,  Phys.  Rev.  ST  Accel.  Beams  8, 
014201  (2005). 

[23]  R.  Bhatt,  T.  Bemis,  and  C.  Chen,  IEEE  Trans.  Plasma  Sci. 
34,  187  (2006). 


104201-8 


PHYSICS  OF  PLASMAS  14.  103102  (2007) 


Warm-fluid  equilibrium  theory  of  a  thermal  charged-particle  beam 
in  a  periodic  solenoidal  focusing  field 

Ksenia  R.  Samokhvalova,3'  Jing  Zhou,  and  Chiping  Chen 

Plasma  Science  atul  Fusion  Center,  Massachusetts  Institute  of  Technology. 

Cambridge,  Massachusetts  02139 .  USA 

(Received  9  May  2007:  accepted  10  August  2007;  published  online  5  October  2007) 

A  warm-fluid  equilibrium  theory  is  presented  which  describes  a  new  thermal  equilibrium  of  a 
rotating  charged-particle  beam  in  a  periodic  solenoidal  focusing  field.  Warm-fluid  equations  arc 
solved  in  the  paraxial  approximation.  The  rms  beam  envelope,  the  density  and  flow  velocity  profiles, 
and  the  self-consistent  Poisson  equations  arc  derived.  Density  profiles  arc  calculated  numerically  for 
high-intensity  and  low-intensity  beams.  Temperature  effects  in  such  beams  arc  investigated.  Radial 


confinement  of  the  beam  is  discussed.  ©  2007  Ameri c 
[DOl:  10.1063/1.2779281] 

I.  INTRODUCTION 

Many  charged -particle  beam  experiments  and  applica¬ 
tions,  such  as  particle  accelerators,  spallation  neutron 
sources,  high-power  microwave  sources  and  high-energy 
colliders,  use  high-intensity  beams  of  charged  particles.  For 
such  systems,  beams  of  high  quality  (i.c..  low  emittancc, 
high  current,  small  energy  spread,  and  low  beam  loss)  are 
required.  Exploration  of  equilibrium  states  of  charged- 
particle  beams  and  their  stability  properties  is  critical  to  the 
advancement  of  basic  particle  accelerator  physics,  as  well  as 
the  design,  construction,  commissioning,  and  operation  ol 
high-brightness  particle  beams  and  accelerator  systems. 

Of  particular  concern  arc  emittancc  growth  and  beam 
losses  which  arc  related  to  the  evolution  of  particle  beams  in 
their  noncquilibrium  states.  To  minimize  emittancc  growth 
and  control  beam  losses,  it  is  critical  to  find  equilibrium  dis¬ 
tributions  of  high-brightness  charged-particle  beams  in  accel¬ 
erators  and  beam  transport  systems. 

Several  kinetic  equilibria  have  been  discovered  for  peri¬ 
odically  focused  intense  charged-particle  beams,  for  ex¬ 
ample,  the  Kapchinsky-Vladimirsky  (KV)  self-consistent 
two-dimensional  Vlasov  equilibrium  for  an  intense  charged- 
particle  beam  in  an  alternating-gradient  quadrupolc  magnetic 
focusing  field  and  the  periodically  focused  rigid-rotor  Vla¬ 
sov  equilibrium  in  a  periodic  solenoidal  magnetic  focusing 
held  Both  equilibria  use  a  /^-function  phase-space  distri¬ 
bution,  which  is  unphysical.  Although  there  is  a  well-known 
rigid-rotor  thermal  equilibrium  in  a  uniform  magnetic  field, 
periodically  focused  thermal  beam  equilibrium  has  not  been 
reported  until  our  present  work,  which  includes  both  a  ki¬ 
netic  treatment  presented  elsewhere6'  and  a  warm-fluid 
treatment  presented  in  this  paper. 

In  this  paper,  we  present  a  warm-fluid  equilibrium  theory 
of  a  new  thermal  charged-particle  beam  in  a  periodic  solc- 
noidal  focusing  field.  Solving  the  warm-fluid  equations  in  the 
paraxial  approximation,  the  beam  density  and  flow  velocity 
are  obtained.  The  self-consistent  rms  beam  envelope  equa¬ 
tion  is  derived.  The  self-consistent  Poisson  equation,  govern- 
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ing  the  beam  density  and  potential  distributions,  is  also  de¬ 
rived.  For  such  thermal  beam  equilibria,  temperature  effects 
are  found  to  play  an  important  role.  Due  to  temperature  ef¬ 
fects,  the  beam  density  profile  has  a  smooth  edge,  which  is  a 
more  realistic  representation  of  the  beam  density  than  the 
uniform  density  profile  in  previous  theories  (sec.  lor  ex¬ 
ample,  Refs.  1,  2,  and  8).  Finally,  we  discuss  the  radial  con¬ 
finement  of  the  beam. 

The  organization  of  this  paper  is  as  follows.  In  Sec.  11, 
the  basic  assumptions  in  the  present  model  are  presented. 
Warm-fluid  equilibrium  equations  are  used  to  derive  expres¬ 
sions  for  the  flow  velocity  profile  and  beam  density  distribu¬ 
tion,  a  root-mean-squared  (rms)  beam  envelope  equation, 
and  a  self-consistent  Poisson  equation.  In  Sec  111,  a  numeri¬ 
cal  technique  for  computing  the  warm-fluid  beam  equilibria 
is  discussed.  Several  examples  of  the  thermal  beam  equilib¬ 
rium  arc  presented.  The  radial  confinement  of  the  beam  is 
demonstrated.  In  See.  TV,  the  conclusions  arc  presented. 

II.  WARM-FLUID  BEAM  EQUILIBRIUM 

We  consider  a  thin,  continuous,  axisymmetrie  (0100 
=  0),  single-species  charged-particle  beam,  propagating  with 
constant  axial  velocity  Vzc.  through  an  applied  periodic  so¬ 
lenoidal  magnetic  focusing  field.  The  applied  periodic  sole¬ 
noidal  focusing  field  inside  the  beam  can  be  approximated  by 

BCI1(r,.v)  =  -  \B'Z  ( s)rc,  +  ( 1 ) 

where  s=z  is  the  axial  coordinate,  r=\x2+y 2  is  the  radial 
distance  from  the  beam  axis,  prime  denotes  the  derivative 
with  respect  to  .y,  and  B:(s)=Bz(s+S)  is  the  axial  magnetic 
field,  which  is  periodic  along  the  c- axis  with  periodicity 
length  S. 

In  the  paraxial  approximation.  rbrms<S  is  assumed, 
where  rbrrns  is  the  rms  beam  envelope.  The  transverse  kinetic 
energy  of  the  beam  is  assumed  to  be  small  compared  to  its 
axial  kinetic  energy,  i.c.,  |  V-l  > \V x  |.  The  Budker  parameter 
of  the  beam  is  assumed  to  be  small,  i.e.,  q2Nbhnc2  <  yb, 
where  q  and  m  arc  the  particle  charge  and  rest  mass,  respec¬ 
tively,  c  is  the  speed  of  light  in  vacuo ,  Nb=27rf^dnnh(r.s)  is 
the  number  of  particles  per  unit  axial  length,  and  yb  is  the 
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relativistic  mass  factor,  which,  to  leading  order,  is 
y/;=con.st=(l with  (ih=Vhlc  =  Vzlc . 

It  is  convenient  to  express  the  self-electric  and  self- 
magnetic  fields,  produced  by  the  space  charge  and  axial  cur¬ 
rent  of  the  beam,  in  terms  of  the  scalar  and  vector  potentials, 
i.e.,  E sclf(r..y)=-V<^hClf(r,.v)  and  Bse,1(r,.v)  =  V X  A5Clf(r^v).  In 
the  paraxial  approximation,  the  self-field  potentials  0self  and 
A*’11  arc  related  by  the  familiar  expression  A%cX^A^ez 
=  /^</>sdf(/\.v)e-.  Consequently,  the  self-magnetic  field  is 

In  the  paraxial  approximation,  the  warm-fluid  cquilib- 
rium  (r9/<9r=0)  equations  arc 


nh\  •  V(y^iV)  =  nhq 


V 

■V<//d,  +  -  x  (B“l  +  BwU) 
c 


-Vp, 


(2) 


V-KV)  =  0,  (3) 

V2^eVf.v)  =  -4m/;j/,(r,.5)l  (4) 

p{i\s)  =  nh(rts)kaT ,  (.v),  (5) 

7’l(-v)rbnns(jr)=  c°nst.  (6) 


In  Hqs.  (2M6),  p(ry.s)  is  the  thermal  pressure,  7\(.v)  is 
the  transverse  beam  temperature  which  remains  constant 
across  the  cross  section  of  the  beam,  and  /‘brmsW  the 
root-mean-square  (mis)  radius  of  the  beam,  defined  by 
1 2 77 Jyd rr^nb( r , .v) .  Equation  (6)  states  that  the 
beam  motion  is  adiabatic.  Note  that  for  the  axisymmctric 
beam  in  the  paraxial  approximation,  we  can  approximate 
V2  =  (l /r)(f)/ f)r)rrl/ dr  to  leading  order  in  the  Poisson  equa¬ 
tion  (4). 

We  seek  a  solution  for  the  equilibrium  beam  profile  of 
the  form 


nb(r,s)  =  /(i)  exp 


2Asri(.v) 


|r"nro(A)  n^)-[n;,(.o  +  n,(,)j 

1  ^brms(A) 


Y,knT  i  (a) 


(10) 


where  /(.v)  is  the  peak  density  at  the  center  of  the  beam, 
which  is  an  arbitrary  function  of  s  to  be  determined  later  [sec 
Eq.  (18)].  Using  the  density  profile  given  in  Eq.  (10).  wc 
obtain  a  useful  expression  for  the  rnts  beam  radius,  i.e.. 


'bmja  =  2 


kBTi  (f) 

,j2  2  + 


| 

1  ^bnm(A  1 


q 1N_b 

2  vlmPbC*. 

U/,(.v)[lUv)  +  U,(.v)] 


-i 


(11) 


where  wc  have  assumed  that  the  beam  density  is  infinitely 
small  at  r=  oc. 

Since  k„T  (s)=(myb/2)((\  -V.)2),-=« ryh({i\x-\\)2)v. 
(A'2)[  =t  Jrim(.v)/2,  wc  can  express  the  rm.s  thermal  emittancc 
of  the  beam  as 


4=(/V)‘V>r<(y.v-^)2>r  = 


kj£_  i(a)4„,,(a) 

2m  yhpf,c2 


(12) 


where  the  statistical  average  of  \  is  defined  in  the  usual 
manner  by  (x)\'=N^1  $  xfi>d*dy  dp  xdps  with  j),  being  the  par¬ 
ticle  distribution  function  corresponding  to  the  warm-fluid 
beam  equilibrium.'5  Combining  Eqs.  (11)  and  (12)  yields  the 
following  rms  beam  envelope  equation 


n»(,v)[nfc(.o  +  nf(f)] 

">  ^bnnsW/ 

fir 


K 

2^*brnw(1') 


(13) 


Vr(r,s)  =  r 


'hrmsfc) 

^brins(^) 


(7) 


V0(r.s)  =  rilh(s),  (8) 

which  corresponds  to  a  beam  undergoing  rotation  with  the 
angular  frequency  0/,(.y)  to  be  determined  sclf-consistcntly 
later  [see  Eq.  (20)]. 

The  radial  component  of  the  momentum  equation  (2)  can 
be  rewritten  as 


—  ln[#i/,(/*.  a)] 
dr 


kBr  (s) 


£cirj™^l  _  n;,(.v) .  [njt,)  +  iic.(.v)] 


'  brms 


(a) 


<><TV,a) 


ylkBTJs) 


dr 


(9) 


where  use  of  Eqs.  (5),  (7 1,  and  (8)  have  been  made  and 
ilc{s)-qB7{s)hncyh  is  the  relativistic  cyclotron  frequency. 
Equation  (9)  can  be  integrated  to  give  the  density  profile 


where  K  =  2Nhq2 1  ypn^c2  is  the  self-field  perveanee. 

Substituting  Eq.  (13)  into  Eq.  (10)  wc  obtain  the  simpli¬ 
fied  expression  for  the  equilibrium  beam  density 


nh(r,s)  =/(.v)exp 


1  — 

K  +  4cfh  ’ 

q<t>x'Hr. a)  \ 

1  44 

. 2  'i,tJA) . 

ylk.T  (,)  J 

(14) 


where  the  scalar  potential  for  the  self-clcctric  field  satisfies 
the  Poisson  equation 


1  d  d  tr 

- r—  <A5tV.v) 

r  dr  dr 


=  -  4my/(.v)cxpj  -  ^3 


ihL 


K 


jgth 

'br, „,(•'') 


r,,s  I  I 

YhkRT  (a)  J 

(15) 


Note  that  when  tf.(.y)  =  const,  the  beam  density  in  Eq.  (14) 
recovers  the  well-known  thermal  rigid-rotor  equilibrium  in  a 
uniform  magnetic  field. 

The  density  profile  in  the  form  of  Eq.  (14)  and  the  ve¬ 
locity  profiles  (7)  and  (8)  have  to  satisfy  the  continuity  Eq. 
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(3).  Substituting  Eqs.  (7),  (8),  anti  (14)  into  Eq.  (3),  and 
integrating  over  the  cross  section  of  the  beam  yields 


d_ 

ds 


'  hrm> 


df(s) 


2  + 

^hrm.s  d-S 


K  d 
4e^/.v 


'  hrms 


U) 


'  hrms 


,(.v) 


<P*x\r%s) 

c}Nh 


=  0. 


(16) 


Note  that  Eq.  (16)  is  equivalent  to  the  conservation  of  the 
total  number  of  particles  per  unit  axial  length,  i.e.. 


dNh 

—  =0, 
ds 


or  Nh  =  const. 


(17) 


Setting  the  sum  of  the  first  two  terms  in  Eq.  (16)  to  zero 

gives 


/(■v)  - 


C 

*  ) 


(18) 


where  C  is  a  constant  of  integration. 

We  solve  the  Poisson  equation  (15)  to  determine  the 
electric  self-field  potential,  with  f{s)  satisfying  Eq.  (18), 
where  the  electric  self-field  potential  on  axis  </>sdt(r=0,.v)  is 
determined  hy  setting  the  sum  of  the  third  and  fourth  terms 
in  Eq.  (16)  to  zero.  A  numerical  scheme  for  determining 
<£se,f(r=0,.v)  will  be  described  in  Sec.  III.  The  electric  self- 
held  potential  energy  on  axis  r/<£sclt(0 ,,v)  is  very  small  com¬ 
pared  to  the  beam  transverse  thermal  energy,  which  will  be 
demonstrated  in  Sec.  III. 

To  gain  further  insight  into  the  azimuthal  motion  of  the 
beam,  we  make  use  of  Eqs.  (7)  and  (8)  to  re-express  the 
azimuthal  component  of  the  momentum  equation  (2),  i.e., 


nh(r,x) 


'bnm(') 

r - 

^brms(^) 


a  a 

pbc~  +  pbt~ 
dr  ds 


ilh(s)r  +  ’-i\.{s)r2 


=  0. 


(19) 


Consistent  with  Eq.  (8),  we  lind  the  solution  to  Eq.  (19)  as 


)  + 


'bnwW  ' 


(20) 


where  ioh  and  rw,  are  constants.  In  Eq.  (20),  the  term 
(°br]iQ^ rbnm^  represents  the  azimuthal  beam  rotation  relative 
to  the  Larmor  frame,  which  rotates  at  the  frequency 
-nc(.v)/2  relative  to  the  laboratory  frame. 

Substituting  Eq.  (20)  into  Eq.  (13),  wc  obtain  the  follow¬ 
ing  alternative  form  of  the  rms  hcam  envelope  equation: 


''hrms)')  + 


K-(.v) 


2  4 


*  brrnsC') 


K  _ 

"^brms(^) 


(21) 


where  \  Kz(s)  =  qB:(s)l2yhph/nc2  is  the  focusing  parameter 
In  the  limit  e.h=0,  Eq.  (21)  recovers  the  previous  envelope 
equation  for  the  eold-fluid  beam  equilibrium.* 

Note  that  the  term  proportional  to  a rlrms(s)  in  Eq 
(21 )  plays  the  role  of  an  effective  cmittance  contribution  to 
the  envelope  equation  associated  with  the  average  azimuthal 


beam  rotation  relative  to  the  Larmor  frame,  and  that  the  rms 
beam  envelope  equation  (21)  agrees  with  the  well-known 
rms  envelope  equation,  (  with  the  interpretation  of  the  total 
cmittance 


et~ 


16^  +  4 


u,brbO 


(22) 


111.  EXAMPLES  OF  WARM-FLUID 
BEAM  EQUILIBRIA 


In  this  section,  wc  present  a  numerical  technique  Re¬ 
computing  the  warm-fluid  beam  equilibria.  Wc  calculate  the 
beam  density  by  solving  the  self-consistent  Poisson  equation 
and  present  several  examples  of  warm-lluid  beam  equilibria. 
Wc  show  that  thermal  beam  equilibria  exist  for  a  wide  range 
of  parameters  and  discuss  the  radial  confinement  of  the 
beam. 

To  determine  the  warm-fluid  beam  equilibrium  numeri¬ 
cally,  we  obtain  the  matched  rms  beam  envelope  b>  solving 
the  rms  beam  envelope  equation  (21 )  with  periodic  boundary 
conditions.  Then,  we  use  the  matched  rms  beam  envelope  in 
the  calculation  of  the  beam  density  and  potential  at  any  given 
s  from  Eqs.  (14)  and  (15). 

We  calculate  the  scalar  potential  for  the  sclf-clec iric  field 
using  the  Poisson  equation  (15).  We  rewrite  the  Poisson 
equation  (15)  as 


]_£ 

r  dr 


dr 


477Z/C 

'hrn  «(•*) 


exp 


q4>xW(r  =  0,5) 
%kBT  (.v) 


f  jL 

44 

<y  A  </>(>.. \)  | 

1  44 

_  2  4m»(i)  _ 

l  (■')  j 

(23) 


where  A<Mr.s)  =  <ft,e,r(r,s)-<£scl,(r=0..v)  and  use  has  been 
made  of  Eq.  (18).  We  solve  Eq.  (23)  with  the  boundary  con¬ 
ditions 


A<M0,s)  =  0 


and 


f?A  </>(  r,  s ) 
dr 


=  0. 


0 


(24) 


We  integrate  Eq.  (23)  from  r=0  to  a  few  /bnm,  paying  special 
attention  to  the  singularity  at  r= 0.  To  avoid  the  singularity, 
wc  analytically  integrate  Eq.  (23)  with  boundary  conditions 
(24)  near  the  z-axis  from  r=  0  to  r=Ar  (with  Ar<^rbrms), 
treating  the  beam  density  as  a  constant  Then,  we  approxi¬ 
mate  A <(>(r,s)  by  the  scalar  potential  of  the  space -charge- 
dominated  beam  with  SK/Eth>  1  as 


A  </>(',  v)  ==  - 


2rtnrn(s)r  ’ 


for  r  =S  v2rbrlns. 


(25) 


Using  this  potential  wc  numerically  integrate  Eq.  (23)  out¬ 
wards  from  r=Ar. 

For  the  purposes  of  numerical  calculations,  it  is  useful  to 
rewrite  Eq.  (17)  as 
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FIG.  1.  Normalized  beam  envelope  profiles  for  S\V5Ul=a0 
+«j  cos(27 ts/S),  o{)=a j  =  1.14.  a  warm- fluid  (solid  curve)  beam  equi¬ 

librium  with  A*=l(),  and  a  cold-fluid  (dashed  curve)  beam  equilibrium  with 

c. 


1 


f 

Jq 


-  1  -  C exp  - 


2imhrdr 

r/0“V  =  O,.v) 


Yhk[)T  (a-) 


r  {  z-2 

~  K  +  *4  " 

(/A <£(/-, .v)  I 

CXP|  “  4  2 

■M)  [ 

.  2  4, !»(•»')  - 

■fonTj.  (.*)  j 

\rdr 


=  0. 


(26) 


In  our  numerical  calculations,  an  iterative  procedure  is  ap¬ 
plied  to  solve  Eq.  (26),  and  A  is  less  than  10“4. 

Figure  1  shows  the  rms  envelope  profiles  for  S\k-(s) 
=a0+d\  cos(27ta75),  </0=£/j=  1.14.  0^=0,  a  warm-fluid 
(solid  curve)  heam  equilibrium  with  the  “scaled1’  normalized 
pcrvcancc  K  =  KS/4?.lh=  10,  and  a  cold-fluid  (dashed  curve) 
beam  equilibrium  with  K=*.  The  rms  beam  radius  rj^s(-v) 
for  the  cold  fluid  beam  equilibrium  is  determined  from  Eq. 
(21)  with  the  right-hand  side  equal  to  zero.  In  Fig.  1,  the 
effects  of  the  Unite  temperature  enlarge  the  rms  beam  enve¬ 
lope  by  1  %. 

In  Fig.  2  we  plot  the  on-axis  electric  self-field  potential 
energy  relative  to  the  heam  transverse  thermal  energy, 
r/0sell(O,.v)/  %kBT  (.y),  as  a  function  of  s/S  for  £=0.1,  1,  and 
10.  The  rest  of  the  parameters  are  the  same  as  in  Fig.  1.  The 
integration  constant  C  is  chosen  such  that  </>self(0,574)=0. 
The  electric  self-field  potential  on  axis  is  indeed  small, 
which  is  consistent  with  the  paraxial  approximation. 

Figure  3  shows  the  density  profiles  for  the  warm-fluid 
(solid  curve)  and  cold-fluid  (dashed  curve)  heam  equilibria 
corresponding  to  the  examples  discussed  in  Fig.  1 .  The 
warm-lluid  beam  density  is  nearly  uniform  up  to  the  beam 
edge  where  it  falls  rapidly  within  a  few  Debye  lengths.  Here, 
the  Debye  length  is  defined  as 


FIG.  2.  Ploi  of  the  on-axis  electric  self-field  poicniial  energy  relative  to  the 
beam  transverse  thermal  energy  as  a  function  of  s/S  for  A  =0.1,  I,  and  10. 
The  olher  system  parameters  are  the  same  as  in  Fig.  1. 


I  yjTtTTT 

\  4m/2/j/j(0,.y) 


(27) 


For  the  warm-fluid  beam  equilibrium,  rhrms ^  1 5.4\7>.  The 
density  of  the  cold  [T  (.v)^0]  beam  is  (see,  for  example. 
Ref.  8) 


n 


cold 

b 


for  /•  =s  \  2 


u/-rkl 

N  briTV 


L(-v). 


(28) 


The  effect  of  the  beam  temperature  on  beam  density  dis¬ 
tribution  is  illustrated  in  Fig.  4.  As  we  increase  the  beam 
temperature  and  keep  other  system  parameters  the  same,  K 
decreases,  and  the  density  profile  makes  the  transition  from  a 
step-function  profile  (for  T  =0)  to  a  bell-shape  profile,  as 
shown  in  Fig.  4. 

There  is  a  wide  range  of  parameters  for  which  the  warm 
fluid  beam  equilibrium  exists  in  a  periodic  focusing  channel. 
For  practical  purposes,  it  is  useful  to  determine  the  radial 
confinement  in  an  average  sense.  In  Fig.  5,  we  plot  the  nor- 


FIG.  3.  Plot  of  the  relative  beam  density  vs.  rfku  for  a  warm-lluid  beam 
equilibrium  (solid  curve)  and  a  cold-fluid  beam  equilibrium  (dashed  curve) 
at  5=0  for  the  same  parameters  as  in  Fig.  1.  Here,  1 5A\it  for  the 

warm-fluid  beam  equilibrium. 
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FICi.  4.  Ploi  of  lhe  relative  density  profiles  at  .9=0  at  several  temperatures: 
K=^  (cold).  10.  5.  and  2.  The  other  system  parameters  are  kept  the  same  as 
in  Fig.  I. 

malized  angular  frequency  of  beam  rotation  in  the  Larmor 
frame,  (S/a,jfyc)(fi£>($)+fit.(jr)/2),  as  a  function  of  the  ef¬ 
fective  self- field  parameter  (se)  =  S2((^2pb(s))/2^o^f3lc  for 
£=0.1,  0.2,  1,  and  10.  The  beam  propagates  in  a  periodic 
solenoidal  focusing  field  with  S\Kz(s)=aQ+aiCOs(2irs/S), 
where  (Iq=(1\  =  1.14.  The  beam  current  is  kept  the  same  while 
the  rms  thermal  emittancc  eth  of  the  beam  decreases  Here, 
o)ph(s)  —  (47rq2nb(Q%s)f  yhm)U2  is  the  plasma  frequency,  rrc, 
=  JqUq2(5)Ja  is  the  vacuum  phase  advance  over  one  axial 
period  5,  the  amplitude  function  h;q(.v)  satisfies  the  following 
equation  (see,  for  example.  Ref.  1 1): 

H’n(.v)  +  kz(s)w0(s)  =  ,  (29) 

W()(.V) 

and  (f(s))= S~ 1  f'lf(s)ds  denotes  the  average  of  the  function 
/(.s)  over  one  axial  period  of  the  system. 

While  Fig.  5  is  computed  for  the  specific  periodic  sole¬ 
noidal  focusing  field  with  S\jKz(s)=aQ+a\  cos(27 ts/S), 


FIG.  5.  Plot  of  lhe  normalized  angular  frequency  of  beam  rolalion  in  the 
Larmor  frame  as  a  function  of  lhe  effective  self-field  parameter  for  normal¬ 
ized  perveances  A'=0.l,  0.2,  1,  and  10. 


FIG.  6.  Plot  of  lhe  critical  effective  self-field  parameter  (a,)  as  a  function  of 
K  =  ACS74c,h.  The  shaded  region  gives  lhe  parameter  space  lor  radial  beam 
confinement. 


where  a()=a]  =  1.14,  vve  observe  no  change  in  Fig.  4  if  wc 
vary  the  values  of  a()  and  provided  that  the  vacuum  phase 
advance  crv  of  the  magnetic  field  docs  not  change.  For  at 
=0,  Fig.  5  recovers  the  thermal  beam  equilibrium  in  a  uni¬ 
form  magnetic  focusing  field  (sec  Ref.  5). 

As  shown  in  Fig.  5,  each  curve  at  a  particular  value  of  A 
has  two  branches.  For  any  value  of  the  effective  self-field 
parameter  (se)  below  a  critical  value,  a  confined  beam  can 
rotate  at  two  angular  frequencies,  either  positive  or  negative 
relative  to  the  Larmor  frame.  For  each  value  of  A',  the  maxi¬ 
mum  (critical)  value  of  the  effective  self-field  parameter  for 
a  confined  beam  is  reached  when  the  beam  docs  not  rotate 
relative  to  the  Larmor  frame.  In  Fig.  6,  the  critical  effective 
self- field  parameter  (se)  is  plotted  as  a  function  of  K 
=  KSf 4elh.  The  parameter  space  for  radial  beam  confinement 
is  indicated  by  the  shaded  region  in  Fig.  6. 

IV.  CONCLUSIONS 

We  presented  a  warm-fluid  equilibrium  beam  theory  of  a 
new  thermal  charged-particle  beam  propagating  through  a 
periodic  solenoidal  focusing  field.  We  solved  the  warm  fluid 
beam  equations  in  the  paraxial  approximation.  We  derived 
the  rms  beam  envelope  equation  and  solved  it  numerically. 
Wc  also  derived  the  self-consistent  Poisson  equation,  govern¬ 
ing  the  beam  density  and  potential  distributions.  We  com¬ 
puted  the  density  profiles  numerically  for  high-intensity  and 
low- intensity  beams.  We  investigated  the  temperature  effects 
in  such  beams.  We  found  that  the  thermal  beam  equilibrium 
has  a  density  profile  with  a  smooth  edge  and  a  uniform  tem¬ 
perature  profile  across  the  beam  cross  section.  Finally,  wc 
discussed  the  radial  confinement  of  the  beam. 
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A  method  is  presented  for  the  design  of  a  high-brightness  nonrclativistie  circular  beam  system 
including  a  charged-particle  emitting  diode,  a  diode  aperture,  a  circular  beam  tunnel,  and  a  focusing 
magnetic  field  that  matches  the  beam  from  the  emitter  to  the  beam  tunnel.  The  applied  magnetic 
field  is  determined  by  balancing  the  forces  throughout  the  gun  and  transport  sections  of  the  beam 
system.  The  method  is  validated  by  three-dimensional  simulations.  €)  2007  American  Institute  of 
Physics .  [ DOl:  10.1063/1.2815938] 


High-brightness,  space-charge-dominated  charged- 
particle  beams  arc  of  great  interest  because  of  their  applica¬ 
tions  in  particle  accelerators,  x-ray  sources,  vacuum  electron 
devices,  and  material  processing  such  as  ion  implantation. 
When  the  beam  brightness  increases,  the  beam  becomes 
space-charge  dominated.  In  the  spacc-charge-dominated  re¬ 
gime,  the  beam  equilibrium  is  characterized  by  a  beam  core 
with  a  transversely  uniform  density  distribution  and  a  sharp 
edge  where  the  beam  density  falls  rapidly  to  zero  in  a  few 
Debye  lengths.  For  particle  accelerators,  high-brightness, 
space-charge-dominated  charged-particle  beams  provide 
high  beam  intensities.  For  medical  accelerators  and  x-ray 
sources,  they  provide  higher  and  more  precise  radiation  dos¬ 
age.  For  ion  implantation,  they  improve  deposition  unifor¬ 
mity  and  speed.  For  vacuum  electron  devices,  they  permit 
high-efficiency,  low-noise  operation  with  depressed 
collectors. 

An  essential  component  of  ehargcd-partiele  beam  sys¬ 
tems  is  the  beam  generation  and  acceleration  diode,  consist¬ 
ing  of  a  charged-particle  emitter  and  a  gap  across  which  one 
or  more  electrostatic  potential  differences  arc  maintained. 
Conventionally,  Pierce-type  diodes  with  or  without  a  grid 
arc  employed  to  produce  the  Child-Langmuir  emission."* 
Compression  is  often  used  in  Pierce-type  diodes  in  order  to 
generate  the  desired  beam  radius.  Scrapers  arc  often  used 
also  to  chop  off  the  nonuniform  beam  edges.4  The  grid,  com¬ 
pression,  and  scrapers  introduce  a  mismatch  into  the  beam 
systems  and  degrade  beam  brightness.  An  important  aspect 
of  charged-particle  beam  system  design  is  the  transition  from 
the  diode  to  the  beam  focusing  tunnel.  While  the  rigid-rotor 
equilibrium  is  well  known  for  a  uniform  solenoidal  focusing 
field,  a  perfect  matching  of  a  circular  beam  from  a  Piercc- 
type  diode  into  the  rigid-rotor  equilibrium  has  not  been  re¬ 
ported  until  this  letter. 

In  this  letter,  a  method  is  presented  for  the  generation, 
acceleration,  focusing,  and  collection  of  a  high-brightness, 
spacc-charge-dominated  circular  charged-particle  beam.  As 
illustrated  in  a  cross  sectional  view  shown  in  Fig.  1,  the 
beam  system  comprises1* 

(1)  a  fiat  circular  emitter  which  emits  charged  panicles, 
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(2)  a  diode  with  one  electrode  at  the  emitter  and  at  least  one 
additional  electrode  which  accelerates  the  charged  par¬ 
ticles, 

(3)  a  beam  tunnel  which  is  connected  electrically  to  at  least 
one  of  the  additional  electrodes, 

(4)  an  applied  axisymmctric  magnetic  field  for  charged- 
panicle  beam  focusing,  and 

(5)  a  depressed  collector  w  hich  collects  the  charged-particle 
beam. 

The  emitter  consists  of  a  fiat,  circular  disk.  A  circular 
chargcd-partielc  beam  is  emitted  from  the  emitter  with  a  uni 
form  density.  The  current  emission  is  space-charge  limited, 
obey  ing  the  Child-Langmuir  law'.  The  electrodes  and  applied 
axisymmctric  magnetic  field  arc  designed  to  preserve  the 
beam  cross  section  in  the  accelerating  section.  The  method 
for  designing  the  required  electrodes  and  applied  axisymmet- 
ric  magnetic  field  is  described  as  follows. 

As  a  first  step,  the  beam  dynamics  is  modeled  with  a 
three-dimensional  OMNITRAK  simulation  with  no  applied 
magnetic  field.  This  provides  the  electric  field  database 
which  is  used  to  compute  iteratively  the  applied  magnetic 
field  required  to  preserve  the  cross  section  of  the  eharged- 
particlc  beam  at  the  accelerating  and  transport  sections. 

As  a  second  step,  an  estimate  of  the  required  applied 
magnetic  field  is  obtained  by  balancing  all  of  the  radial 
forces  on  the  charged  particles  on  a  line  whose  radius  corre¬ 
sponds  to  the  root-mean-square  (rms)  radius  of  the  emitter. 
The  line  starts  at  the  cathode  disk  and  continues  through  the 
anode  aperture. 


FIG  L  Cross  sectional  view  of  a  high-brightness  circular  charged- panicle 
beam  system 
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For  a  thin  charged-particle  beam,  all  forces  are  propor¬ 
tional  to  the  radial  displacement,  and  an  expression  for  the 
required  magnetic  field  can  be  derived  to  achieve  the  radial 
force  balance  at  any  radius  in  the  beam  core.  An  applied 
axisymmetric  magnetic  field  is  expressed  in  terms  of  the  vec¬ 
tor  potential  in  a  cylindrical  coordinate  system  as 


1  x  .  r)A(} 

I  l  =  c.--(M„)-er— 

r  r)r  flz 

In  the  thin  beam  approximation, 

„  I  DB,(0,z)  . 

1$  =  B.(  0.r.)e,  - - — -re,. 

2  Jz 

Inspection  of  Eqs.  (1 )  and  (2)  gives 
A^r,z)=r-B.(  O.z). 


(1) 

(2) 

(3) 


The  conservation  of  canonical  angular  momentum  y  ields  the 
expression  for  the  azimuthal  velocity 

t >dr,z)  =  -  ~[B,( O.z)  -  B.(0,0)]r.  (4) 

2m 


where  use  has  been  made  of  i^(r,0)=0  at  the  emitter,  and  q 
and  m  arc  the  particle  charge  and  rest  mass,  respectively.  The 
radial  force  balance  equation  is 


< lEr(r,z )  ifor.z)  q 

—  -=- - -  —  VfHr,z)B.(r.z) 

m  r  m 


(5) 


Substituting  Eq.  (4)  into  Eq.  (5)  yields 


qErir,z) 
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b:{  o,z). 


(6) 


where  use  has  been  made  of  the  approximation  ZT(r,c) 
=Bz(0,z)  and  the  boundary  condition  ZT(0,0)  =  0.  Equation 
lb)  produces  a  relationship  between  Erlr  and  #z(0,z). 

By  iterating  the  second  step  described  above,  better  es¬ 
timates  of  the  required  magnetic  field  arc  obtained.  Typically, 
results  converge  after  two  or  three  iterations.  As  an  illustra¬ 
tive  example.  Fig.  2  shows  the  final  results  of  OMNITRAK 
simulations  with  a  grid  with  A.r=  Av=0.05  mm  and 
Az=().l  mm  for  a  circular  electron  beam  which  is  emitted 
from  an  emitter  (cathode)  with  a  radius  of  1  52  mm.  a  cur¬ 
rent  of  0.1  1  A,  a  cathodc-to-anodc  distance  of  4.1  1  mm  at 
radius  r-  1.52  mm.  The  diode  voltage  is  2300  V.  In  Fig.  2, 
the  electrodes  at  the  cathode  and  the  anode  arc  equipotential 
surfaces  which  are  analytically  computed  to  yield  a  Child- 
Langmuir  flow  in  the  absence  of  the  anode  aperture.  A 
circular  aperture  at  the  anode  is  chosen,  and  it  has  a  radius  of 
1.8  mm.  A  larger  circular  beam  tunnel  with  a  radius  of 
6.0  mm  is  connected  to  the  anode.  The  total  axial  length  of 
the  system  is  10  mm. 

Figure  3  shows  plots  of  Er  vs  z  at  r=  1 .0  mm  and  /T(0,z) 
vs  ^  from  the  OMNITRAK  simulations  after  two  iterations.  The 
radial  electric  field  vanishes  al  the  emitter  (i.e.,  at  z=0), 
achieves  a  maximum  magnitude  at  z= 4.1 1  mm.  and  then  ap¬ 
proaches  to  a  constant  value  well  inside  the  beam  tunnel.  The 
applied  axial  magnetic  field  at  the  emitter  vanishes,  increases 
to  640  G  al  the  aperture,  and  then  falls  to  about  160  G  well 
inside  the  beam  tunnel.  In  the  simulation,  Eq.  (6)  is  exactly 
satisfied  for  all  the  electrons  on  the  rms  beam  radius,  and 
approximately  satisfied  for  the  electrons  not  at  the  rms  beam 


0  z  (mm)  10  0 


FIG.  2.  (Color  online)  omnitrak  simulation  of  the  dynamics  of  an  electron 
beam  emitted  from  a  flal  circular  calhode  with  a  radius  of  1.52  mm  and  a 
currenl  of  0.1 1  A  in  an  optimized  magnetic  field.  Here,  the  cathodc-to-anodc 
distance  is  4.1 1  mm  at  /  =  1.52  mm.  the  circular  anode  aperture  has  a  radius 
of  1 .8  min,  and  the  heam  tunnel  has  a  radius  of  6.0  mm.  The  diode  collage 
is  2300  V. 


radius.  The  difference  in  the  electric  and  applied  magnetic 
forces  in  Eq.  (6)  results  in  less  than  2 %  deviations  of  the 
electron  radii  from  their  initial  values.  Figure  4  shows  the 
electron  distribution  in  the  phase  plane  (v,y)  at  .=8  mm. 
The  beam  distribution  maintains  transverse  uniformity.  In 
deed,  the  normalized  fourth  moment  (r*)l{r2)2  which  is 
equal  to  4/3  for  a  transversely  uniform  density  distribution, 
remains  4/3  within  ±1%.  The  outer  beam  radius 
remains  the  same  as  the  emitter  radius  within  l%-3%. 
At  z= 8  mm,  the  normalized  angular  rotation  velocity 
is  (ob={x,y)/(y2)= -0.0244  rad/mm,  and  the  unnormalized 
rms  thermal  emitiances  in  the  x-  and  y-direelions 
arc  E{rh  =  \ !(x2)((x'  -  (oby)2)=  1 .96  mm  mrad  and  kJJ1 
=  \(y2)((y'  +  ov02)=  1  96,  respectively,  which  correspond  to 
a  normalized  rms  thermal  omittance  of  #4h  =  —  /V-f 

=0.186  mm  mrad.  Adding  this  (numerical)  emittance  to  the 
intrinsic  normalized  rms  emittance  of  0.33  nun  mrad  for  a 
thermionic  cathode  of  radius  1,52  mm  predicts  a  normalized 
thermal  emittance  of  e®c,ua1^  0.38  mm  mrad  for  the  actual 
electron  beam. 

Because  the  beam  is  in  laminar  flow,  a  depressed  collec¬ 
tor  is  designed  using  the  same  geometry  as  the  charged- 
particle  emitting  diode,  where  the  circular  emitting  disk  is 


FIG.  3.  Plots  of  Er  vs  z  al  r-  !  mm  (solid  curve)  and  B  (0.:)  vs  :  (dashed 
curve)  from  ihe  omnitrak  simulations  afler  two  iterations. 
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FIG.  4.  Plot  of  the  electron  distribution  at  c=K  mm  in  (he  phase  plane  (,v,.v). 

the  beam  collecting  surface,  and  diode  voltage  is  slightly 
lower  (i.e.,  a  fraction  of  a  percent  lower)  than  the  diode 
voltage  but  has  a  negative  bias.  Such  a  depressed  collector 
yields  a  collection  efficiency  of  nearly  100%. 

To  summarize,  a  high-brightness,  spacc-chargc- 
dominated  circular  chargcd-particlc  beam  system  was  de¬ 


scribed.  The  method  was  presented  for  the  system,  including 
the  beam  generation,  acceleration,  focusing,  and  collection 
processes.  The  method  was  validated  by  three-dimensional 
simulations. 
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An  adiabatic  equilibrium  theory  is  presented  for  an  intense,  axisymmetric  charged-particle  beam 
propagating  through  a  periodic  solenoidal  focusing  field.  The  thermal  beam  distribution  function  is 
constructed  based  on  the  approximate  and  exact  invariants  of  motion,  i.e..  a  scaled  transverse 
Hamiltonian  and  the  angular  momentum.  The  approximation  of  the  scaled  transverse  Hamiltonian 
as  an  invariant  of  motion  is  validated  analytically  for  highly  ennttancc-dominated  beams  and  highly 
space-charge-dominated  beams,  and  numerically  tested  to  be  valid  for  cases  in  between  with 
moderate  vacuum  phase  advances  {crv  <90°).  The  beam  root-mcan-square  (rms)  envelope  equation 
is  derived,  and  the  self-consistent  nonuniform  density  profile  is  determined.  Other  statistical 
properties  such  as  How  velocity,  temperature,  total  emittance  and  mis  thermal  emittance,  equation 
of  state,  and  Debye  length  arc  discussed.  Numerical  examples  are  presented,  illustrating  the  effects 
of  the  beam  pervcancc,  emittance,  and  rotation  on  the  heam  envelope  and  density  distribution.  Good 
agreement  is  found  hetween  theory  and  a  reeent  high- intensity  beam  experiment  performed  at  the 
University  of  Maryland  Electron  Ring  [S.  Bernal,  B.  Quinn,  M.  Reiser,  and  P.  G.  O’Shea,  Phys.  Rev. 

ST  Accel.  Beams  5,  064202  (2002)].  ©  2008  American  Institute  of  Physics. 
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I.  INTRODUCTION 

A  fundamental  understanding  of  the  kinetic  equilibrium 
and  stability  properties  of  high-intensity  electron  and  ion 
beams  in  periodic  focusing  fields  plays  a  central  role  in  high 
energy  density  physics  research,  in  the  design  and  operation 
of  particle  accelerators,  such  as  storage  rings,  rf  and  induc¬ 
tion  linacs,  and  high-energy  colliders,  as  well  as  in  the  design 
and  operation  of  vacuum  electron  devices,  such  as  klystrons 
and  traveling-wave  tubes  with  periodic  permanent  magnet 
(PPM)  focusing.  For  such  systems,  beams  of  high  quality 
(i.e.,  low  emittance,  high  current,  small  energy  spread,  and 
low  beam  loss)  are  required.  Exploration  of  equilibrium 
states  of  charged-particlc  heams  and  their  stability  properties 
is  critical  to  the  advancement  of  basic  particle  accelerator 
physics. 

Several  kinetic  equilibria  have  been  discovered  for  peri¬ 
odically  focused  intense  charged-particle  beams.  Well-known 
equilibria  for  periodically  focused  intense  beams  include  the 
Kapchinskij-Vladmirskij  (KV)  equilibrium  in  an 
alternating-gradient  quadrupole  magnetic  focusing  field  and 
the  periodically  focused  rigid-rotor  Vlasov  equilibrium4  in  a 
periodic  solenoidal  magnetic  focusing  field.  Both  beam 
equilibria  *  have  a  singular  (^function)  distribution  in  the 
four-dimensional  phase  space.  Such  a  <$-function  distribution 
gives  a  uniform  density  profile  across  the  beam  in  the  trans¬ 
verse  directions,  and  a  transverse  temperature  profile  that 
peaks  on  axis  and  decreases  quadrat ically  to  zero  on  the  edge 
of  the  heam.  Because  of  the  singularity  in  the  distribution 
functions,  both  equilibria  are  not  likely  to  occur  in  real 
physical  systems  and  cannot  provide  realistic  models  for  the¬ 
oretical  and  experimental  studies  and  simulations  except  for 
the  zero- temperature  limit.  For  example,  the  KV  equilibrium 
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model  cannot  be  used  to  explain  the  beam  tails  in  the  radial 
distributions  observed  in  recent  high-intensity  beam 
experiments/ 

In  general,  beams  tend  to  relax  to  a  thermal  equilibrium 
in  the  transverse  plane  such  that  the  temperature  across  the 
transverse  beam  cross  section  is  uniform.  A  theoretical  un¬ 
derstanding  of  thermal  equilibrium  and  stable  transport  is 
desirable.  Kinetic  and  warm-lluid  theories  of  a  thermal  equi¬ 
librium  in  a  uniform  magnetic  focusing  field  have  been  stud¬ 
ied  in  Ref.  3.  A  formal  multiple  scale  analysis  (a  third-order 
averaging  technique)  has  been  applied  to  obtain  an  approxi¬ 
mate  periodically  focused  thermal  equilibrium  in  periodic 
solenoidal  and  periodic  quadrupole  magnetic  fields/  Such  an 
averaging  procedure  is  expected  to  be  valid  for  sufficiently 
small  vacuum  phase  advances.  The  concept  of  adiabatic 
breathing  of  high-intensity  charged-particlc  heam  was  pro¬ 
posed  in  Ref.  7.  Recently,  a  warm-lluid  equilibrium  theory  of 
a  periodically  focused  intense  charge-particle  beam*  has 
been  developed  under  the  assumption  of  the  adiabatic  equa¬ 
tion  of  state.  In  this  paper,  a  kinetic  thermal  equilibrium  in  a 
periodic  magnetic  focusing  field  is  developed  to  provide 
more  insight  into  the  equilibrium  properties  of  a  periodically 
focused  thermal  beam. 

We  present  a  kinetic  theory  describing  an  adiahatic  ther¬ 
mal  equilibrium  of  an  intense  charged-particlc  beam  propa¬ 
gating  through  a  periodic  solenoidal  magnetic  focusing  field. 
For  continuous  beams  with  long  pulses,  the  longitudinal  en¬ 
ergy  spread  is  small  such  that  the  longitudinal  motion  can  be 
treated  as  “cold”  and  decoupled  from  the  transverse  motion, 
which  is  kept  nonrelalivistic.  The  beam  pulsates  in  transverse 
directions  adiabatically  like  an  ideal  gas  in  an  adiabatic  pro¬ 
cess,  in  which  the  invariant  is  the  product  of  the  transverse 
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temperature  and  an  effective  beam  area.  It  differs  from  the 
usual  thermal  equilibrium  in  which  the  temperature  is  kept 
constant  In  the  present  work,  the  Hamiltonian  for  single- 
particle  motion  is  analyzed  to  nnd  the  approximate  and  exact 
invariants  of  motion,  i.e.,  a  scaled  transverse  Hamiltonian 
(nonlinear  space  charge  included)  and  the  angular  momen¬ 
tum,  from  which  the  beam  equilibrium  distribution  is  con¬ 
structed  The  approximation  of  the  sealed  transverse  Hamil¬ 
tonian  as  an  invariant  of  motion  is  validated  analytically  for 
highly  cmittance-dominatcd  beams  and  highly  spaec-ehargc- 
dominated  beams,  and  numerically  tested  to  be  valid  for 
cases  in  between  with  moderate  vacuum  phase  advances 
(rrL.  <90°).  The  beam  envelope  and  cmittanecs  arc  deter¬ 
mined  sclf-consistcntly  with  the  beam  equilibrium  distribu¬ 
tion.  Because  the  distribution  function  has  a  Maxwell- 
Boltzmann  form,  it  solves  not  only  the  Vlasov  equation  but 
also  the  Fokkcr-Planek  equation  It  is  expected  to  be  stable  in 
a  similar  manner  as  the  beam  thermal  equilibrium  in  a 
smooth- focusing  approximation.9 

The  paper  is  organized  as  follows.  In  See.  II,  the  theo¬ 
retical  model  is  introduced.  Exact  and  approximate  constants 
of  motion  are  found  for  the  single-particle  Hamiltonian  in  the 
paraxial  approximation,  and  the  equilibrium  distribution  is 
constructed.  In  Sec.  Ill  the  statistical  properties  of  the  beam 
equilibrium,  such  as  the  beam  envelope  equation,  emittanccs, 
and  beam  temperature,  arc  discussed.  In  Sec.  IV,  the  numeri¬ 
cal  calculations  of  the  beam  density  and  potential  arc  pre¬ 
sented  In  Sec.  V,  the  comparison  between  our  theory  and  a 
recent  experiment  from  University  of  Maryland  Electron 
Ring  (UMER)  is  discussed.  Finally,  the  paper  is  concluded 
in  Sec.  VI 

II.  BEAM  EQUILIBRIUM  DISTRIBUTION 

We  consider  a  continuous,  intense  charged-particle  beam 
propagating  with  constant  axial  velocity  /3bce.  through  an 
applied  periodic  solenoidal  magnetic  focusing  field.  The  pe¬ 
riodic  solcnoidal  magnetic  focusing  field  is  described  by 


=  {m2c4+c2P:)1 2.  Because  v!  yhl^b<  1,  the  longitudinal  par¬ 
ticle  motion  can  be  decoupled  from  the  transverse  particle 
motion,  and  the  total  Hamiltonian  for  single-particle  motion 
is  approximated  by  H~ybmc2  +  HL.  Here,  the  longitudinal 
Hamiltonian  H^=ybmc2  is  a  constant  and  the  normalized 
transverse  Hamiltonian  H  =//,  / y^nfiyc1  is  expressed  as 


H x (.v.v. Px, Py,s)  =  |{[PV  +  \  *.(*) vf  +  [Pv  -  \k .(ihf } 


+  <r'", 


(2) 


where  \j  k\(s)  =qB:(s)l2ybm/3bc 2  is  a  measure  of  the 
strength  of  periodic  solcnoidal  magnetic  field,  P 
=  P Llyhmfihc.  ^eU=q<ffclilylml3ici=K^fl2(iNh.  and  K 
=  2q2Nh/  Yjjnfiic1  is  ihc  beam  pervcanec.  The  scalar 
and  vector  potentials  for  the  self-electric  and  self-magnetic 
fields  satisfy  V2  <p*lf=-4irqnb(x,y  .jr),  V2  Ase,f 
=-4Tr/3bcqiib(x,  \\s)ez,  and  A sdt=A!d  e.=/^<£sd,e-.  The  as¬ 
sociated  equations  of  motion  with  the  Hamiltonian  Eq.  (2) 
are 


cfx 

(is1 


- dy 

2  \*z(s)-f 

ds 


d 


\  K.js) 
ds 


v  + 


dx 


=  0, 


(3) 


d2v  - dx  d\ K~(s)  d<j)**^ 

— t  +  2\  k-(.v) —  + - - — x  + - =  0. 

ds2  ds  ds  dy 


(4) 


In  order  to  simplify  the  transverse  Hamiltonian 
H j  (x,\\PxyPVys),  we  perform  a  two-step  canonical  transfor¬ 
mation  The  first  step  is  to  transform  from  the  Cartesian  co¬ 
ordinates  into  the  Larmor  frame  which  rotates  with  one-half 
of  the  cyclotron  frequency  relative  to  the  laboratory  frame. 
The  second  step  is  a  Courant-Snyder  type  of  transformation. 
The  first  transformation  uses  the  second  type  of  the  generat¬ 
ing  function 


F2(a-.)> ;  Px,  Py,x)  =  [x  cos  c p{s)  -  y  sin  <p(.v)]P, 


,  I  cW,(s) 

BCM  =  «,(.v)e:  -  -  ~  (ac,  +  vev) ,  ( I ) 

2  ds 

where  s=z  is  the  axial  coordinate,  Bz(s)  =  B,(s+S)  is  the  axial 
magnetic  field,  and  S  is  the  fundamental  periodicity  length  of 
the  focusing  field.  The  condition  S>rhrms  is  assumed,  where 
r brtns  is  the  characteristic  radius  of  the  rms  beam  envelope. 

The  single-particle  Hamiltonian  can  be  written  as  H 
=  [wrr4  +  (ci>—<7A)2]I/2+f/0sel*.  where  the  canonical  momen¬ 
tum  P  is  related  to  the  mechanical  momentum  p  by  P=^ 
+  r/A/c,  A=Aexl+ Asdf  is  the  vector  potential  for  the  total 
magnetic  field,  Asdl  is  the  vector  potential  for  the  self- 
magnetic  field,  Aexl(A,  v,A*)=fi2U)(-yev+Aev)/2  is  the  vector 
potential  for  the  applied  magnetic  field,  0sdt  is  the  scalar 
potential  for  sclf-clectrie  Held,  in  and  q  arc  particle  rest  mass 
and  charge,  and  c  is  the  speed  of  light  in  vacuum. 

In  the  paraxial  approximation,  we  assume  vlyjfi2h<  1, 
where  r^q2Nbhnc2  is  the  Budker  parameter  of  the  beam, 
Nb=S*nh{.x,y)dxdy= const  is  the  number  of  particles  per  unit 
axial  length,  and  yh=(  I -$)~Xl2  is  the  relativistic 
mass  factor.  The  axial  energy  is  approximately  ybmc2 


+  [x  sin  <p(.y)  +  y  eos  tp(s)]Py% 
where  q>(s)=p{)\Kz(s)ds.  The  transformation  is 


v  = 


.v  cos  ip(s )  -  v  sin 


dPx 

db  i 


— 11  =x  sin  <p(s)  +  v  cos  <p(s), 
Wy 


(5) 

(6) 

(7) 


Px  eos  ip(s)  +  Py  sin  ^(.v). 


(8) 


dFi 

PY  =  — -  =  ~  Px  sin  <p(s)  +  PY  cos  <p(s).  (9) 

dy 

The  transverse  Hamiltonian  after  the  11  rst  transformation  is 
expressed  as 
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-  ~  dF-> 

H  Lv,v,A\,/Vj)  =  //  (a ,v,Pf,/V)  +  — = 

ds 

=  “[^7  +  +  ^-(v)Cr +.V2)]  + 

(10) 

Note  that  (fil  dxl+(Pl  dy2)^x\x,y  ,s)  =  (&!  dxl 

+  r>!/rVv2)0selt(r,  v,.5).  The  equation  of  motion  associated  with 
the  transverse  Hamiltonian  in  Eq.  (10)  is 


tPx  _ 
— 7  +  KZ(S)X  + 
ds~ 


dx 


=  0, 


(ID 


dy  _  d<f> 
->  +  k7(s)\  + 
ds~  '  d\ 


self 


=  0. 


(12) 


The  second  canonical  transformation  uses  the  second 
type  of  the  generating  function 


F2(xJ;P„Pv.s)  = 


w(s) 


1  dw(s) 

P,  + - x 

A  2  ds 


I  dw(s) 

'V,  —y 


vv(.y)  L  2  ds 

where  vv(.v)  satisfies  the  differential  equation 

K 


d~\v(s) 

— Ti  +  ^(.v)wlv)- 
ds~ 


2/-L«(v) 


U'(.v)  = 


»■(•?) 


(13) 


(14) 


and  t’brms(s)  is  the  rms  beam  radius.  It  will  be  shown  in  See. 
Ill  that  the  function  u’(.v)  is  related  to  the  rms  beam  radius 
[see  Eq.  (29)],  The  transformation  is 


dx_dH__ 
ds  jpx 


(20) 


dy  dH  PY 

ds  ap,  w'2(y)  ’ 


dPx  dH±  x_  K  , 

ds  dx  w2(s)  Px  2r^nu(.v)'  '  ^ 


(22) 


dP,_  9H^_  l  _  _  K  , 

ds  ay  w*(s)  ay  2rllnJs)" 


(23) 


In  order  to  construct  a  beam  equilibrium  distribution,  we 
need  to  hnd  constants  of  motion  of  the  system.  Two  con¬ 
stants  of  motion  can  be  found  using  the  transverse  Hamil¬ 
tonian  in  Eq.  (19).  It  is  readily  shown  that  the  angular  mo¬ 
mentum  P0  is  a  constant  of  motion,  i.c.. 


dPe 

ds 


•  =  0. 


(24) 


In  deriving  Eq.  (24),  we  have  used  Eqs.  (20)— (23),  and  the 
axial  symmetry  property  of  the  self- held  potential,  i.c.,  </>scH 
is  only  a  function  of  r=  \  x 2+y2  and  s. 

We  also  hnd  that  the  scaled  transverse  Hamiltonian  for 
single-particle  motion, 


£=  w2(s)H  l  (a,v,A\,Pv,.v), 


(25) 


sf2 

X 

<y\ 

U’(.V)' 

#2 

V 

*Py 

H’(.v)  ’ 

#2 

I 

dx 

H'(.V) 

(IF  2 

I 

dy 

\v(s) 

Jw(s) 

+  X - 

ds 


5  -dw(s) 
Py+>— — 
ds 


(15) 


(16) 


(17) 


(18) 


Using  Eqs.  ( 1 53— ( 1 8).  the  transverse  Hamiltonian  is  trans¬ 
formed  into 


H  (V. v,/\.  Prs)  =  -T- [p;  +  P;  +  x2  +  V2]  +  d>xlf(x,y,s) 
2h"(a) 


+  ,  J—.  n2(.v)(A2  +  y2).  (19) 

4'^ms(y) 


The  equations  of  motion  associated  with  the  Hamiltonian  in 
Eq.  (19)  are 


is  an  approximate  invariant  The  transverse  Hamiltonian  is  a 
highly  oscillating  function.  We  use  the  periodic  I  unction 
w2(s)  to  scale  the  transverse  Hamiltonian  and  to  eliminate 
the  oscillations  such  that  the  scaled  transverse  Hamiltonian  is 
an  approximate  invariant  with  small  residual  oscillations.  As 
will  be  discussed  in  Sec.  IV,  the  small  residual  oscillations 
are  numerically  estimated  to  be  a  few  percent.  Using  Eqs. 
(20) — (.23),  the  derivative  of  the  sealed  transverse  Hamil¬ 
tonian  can  be  evaluated,  using 


clE 

ds 


d_ 

ds 


+  f  +  P2  +  P2)  +  w2(.v)  <r'lf(r„y) 


*»■*(»)  _,1 
4 rWm(s)  J 

4  H2(.v)<r'f(F..r)  + 


Kw\s) 

4rhrms('y)r  . 


(26) 


It  is  readily  shown  that  dElds  is  approximately  zero  in 
two  limiting  cases:  (a)  a  highly  spaec-chargcd-doininatcd 
beam  with  SK/f 1  and  (b)  a  highly  emittanee-doininated 
beam  with  SK/e^—* 0,  where  slh  is  the  thermal  beam  cmit- 
tance  defined  later  in  Eq.  (35).  For  a  highly  space-charge- 
dominated  beam  with  SK/Elh>  1, 
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&&7.S)  =  -  rfr-/’  ^  '•  <  v2^n«  (27) 

4rbrms(*v) 

and  dElds=Q.  For  a  highly  ernittancc-dominatcd  beam  with 
SK/e^  —  O,  ^ nih  ~ 0  and  dE/ds=0.  For  eases  in  which  the 
space-charge  effect  is  comparable  to  the  emittancc  effect,  we 
will  demonstrate  dElds^O  numerically  in  See.  IV. 

We  choose  the  beam  equilibrium  distribution  in  the  form 
similar  to  the  Maxwell-Boltzmann  distribution,  i  c., 

fh  =  Ccw[-(i(E-<ohPe)l  (28) 

where  C,  /3,  and  u)b  are  constants.  C  is  an  integration  con¬ 
stant,  is  related  to  the  beam  emittancc,  and  toh  is  the  rota¬ 
tion  frequency  relative  to  the  Larmor  frame.  Note  that  u)h 
=0  for  Brillouin  flow  and  a)h^  0  for  general  Hows  in  which 
there  is  magnetic  flux  on  the  emitter.  Since  P 0  and  E  are  the 
constants  of  motion,  the  distribution  function  fh  defined  in 
Eq.  (28)  is  indeed  a  Vlasov  equilibrium,  i.c.,  dfb/ds= 0. 


III.  STATISTICAL  PROPERTIES 


In  this  section,  we  will  discuss  the  statistical  properties 
of  the  thermal  equilibrium  developed  in  Sec.  II,  including  the 
rms  beam  radius,  rms  emittancc  and  thermal  emittancc,  flow 
velocity,  beam  temperature,  and  beam  density  profile. 

The  distribution  function  described  in  Eq.  (28)  has  the 
following  statistical  properties.  First,  the  rms  beam  radius 

^brm,(')  -  <*2  +y2)=^h '  JJ  (r  +y1)fbdxdydPxdPy  can  be 

evaluated  to  yield 

'br.mW  =  ,  ~  2  w2(jt)  ==  yw2(j),  (29) 

p(  1  -  (0b)  2 


where  we  have  introduced  the  concept  of  the  total  emittancc 
t:f =4/r‘(l -<uy“)_1.  Substituting  Eq.  (29)  into  Eq.  (14),  we 
arrive  at  the  rms  envelope  equation. 


^  ^brim  (  v  ^ 

2  *AA)^brms~~  ^  _ 


dsl 


'  brm> 


4  r 


(30) 


brms 


Second,  the  rms  beam  emittancc  of  the  beam  equilibrium 
described  in  Eq.  (28)  is  given  in  the  Larmor  frame  by 

e*nu  =  \  <-V2)<JC,2>  -  (xx’)2  =  k,74  =  const.  (3  1 ) 


Similarly,  ezmx$=ET/4=consl.  Note  that  Eq.  (30)  agrees  with 
the  well-known  rms  envelope  equation  in  Ref.  10  with  the 
interpretation  of  the  total  omittance  in  Eq.  (29).  As  a  third 
statistical  property,  in  dimensional  units,  the  average  (mac¬ 
roscopic  flow)  transverse  velocity  of  the  beam  equilibrium  is 
given  in  the  Larmor  frame  by 


v  (r.s)  =  [«/(.v)]-1  J  J  V  JhclPKdPv 


'brms^ 

^bnnsO'") 


ftbcrcr + 


2rinJ*) 


-  V  k-(.v) 


ftnc,  (32) 


As  the  fourth  property,  the  beam  equilibrium  described  by 
Eq.  (28)  has  the  transverse  temperature  profile  (in  dimen¬ 
sional  units) 


kBnU)  =  /  (yi  -  V  ■  ) 

(1  -  o^myb^E} 

S'-Lsb)  ' 

where  kp  is  the  Bolt/.mann  constant.  Note  from  Eq.  (33)  that 
the  product  'P1(s)r^rnit(s)  is  a  conserved  quantity  (d/ds= 0) 
as  the  beam  pulsates  transversely;  that  is,  the  equation  of 
stale  is 

7*1  W'iiwW®  const.  (34) 

Since  2 7rr^.ms  is  a  measure  of  the  effective  area  of  the  beam, 
Eq.  (34)  is  analogous  to  the  equation  of  state  for  a  two- 

t  o 

dimensional  adiabatic  plasma.  As  the  lifth  property,  the 
thermal  beam  emittancc  in  the  Larmor  frame  is 


«*■ 


-T;(x2)<(t\-  Vt)2)  = 

Pif 


^rT  Cv)rhrms('v) 

2  Vb'nfoc1 


=  const . 


(35) 


It  follows  from  Eqs.  (31)  and  (35)  that  r:)=  I6e?h 
+  16 where  the  term  16e>£«^m  corresponds  to  the 
contribution  from  the  average  azimuthal  motion  in  ihe  Lar¬ 
mor  frame  to  the  total  emittancc.  The  rms  envelope  equation 
(30)  can  also  be  expressed  as 


f  brms  (-y) 

ds2 


-5r[n»w + awiuj)  - 
fir 


K 


bomb') 


4«n 


.3 

'  brms 


b)’ 


(36) 


where  i'lh(s)  =  atheTfihc/2rfirmi{s)-il,  .(.v)/2  and  11,. (.v) 

=qB.{s)/  y hmc.  Finally,  the  beam  density  profile  is 


nh(r,s)  =  w1_2(.v) 


fdPJPy 


47rC4 

4n»(*) 


K 

~  + 
2 


4«n 


'  bnns 


b)  J 


r 


<1 


(37) 


where  the  scalar  potential  for  the  scif-cicctric  held  is  deter¬ 
mined  by  the  Poisson  equation 


L  — 

r  dr 


r— ~  j  =-4minb(r,s). 


(38) 


It  is  worth  noting  that  in  the  paraxial  approximation,  the  total 
number  of  particles  per  unit  length  is  kept  constant,  i.e.,  Nb 
— 5^nb(r,s)2Tirdr.  This  requires  the  on-axis  self-electric  po¬ 
tential  <£selt(r=(),.v)  to  vary  as  a  function  of  the  longitudinal 
distance,  which  will  be  determined  numerically  in 
Sec.  IV. 

It  is  readily  shown  that  the  thermal  equilibrium  distribu¬ 
tion  in  Eq  (37)  recovers  the  well-known  thermal  rigid-rotor 
equilibrium  in  a  constant  magnetic  focusing  field  ^  by  set¬ 
ting  d2rbrins! ds2  =  0  in  Eq  (36). 

It  is  worth  pointing  out  that  because  the  derivation  of  the 
theory  docs  not  assume  specitic  magnetic  proiile  as  defined 
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TIG.  1.  Plots  of  the  normalized  axial  magnetic  lield  k\(j)  (dashed  curve) 
and  beam  mis  envelope  rhrmsf  \4cth5  (solid  curve)  versus  the  axial  propaga¬ 
tion  distance  s/S  for  a  periodically  focused  adiabatic  thermal  beam  equilib¬ 
rium  in  an  applied  magnetic  field  described  by  the  periodic  step-function 
lattice  in  Fq.  (39 1.  Here,  the  choice  of  system  parameters  corresponds  to 
Sv#r.„=2.l2,  77=0.4  (£7U=80°),  SK/4t:lh=  10.0,  and  <a*=0. 


FIG.  2.  Plot  of  the  relative  beam  density  versus  r/\D  at  .9=0  for  the  same 
beam  and  focusing  lield  as  in  Fig.  I.  Here.  rbnm=l6. \\f)  and  the  beam 
densities  are  normalized  to  the  peak  density. 


in  Eq.  (1),  it  is  valid  not  only  for  the  periodic  solenoid  mag¬ 
netic  held  hut  also  for  arbitrary  varying  solenoid  magnetic 
field.  Therefore,  our  results  apply  for  the  periodic  focusing 
channel  as  well  as  for  the  matching  section  between  the 
source  and  the  periodic  focusing  channel.  We  will  discuss 
numerical  examples  in  a  periodic  foeusing  channel  in  Sec. 
IV  and  compare  the  theoretical  results  with  the  UMER  ex¬ 
perimental  measurements  in  a  short  matching  solenoid  chan¬ 
nel  in  See.  V. 


IV.  NUMERICAL  CALCULATIONS 


In  this  section,  we  illustrate  examples  of  adiabatic  ther¬ 
mal  beam  equilibria  in  a  periodic  solcnoidal  focusing  field 
and  the  temperature  and  beam  rotation  effects  with  numerical 
calculations.  We  also  demonstrate  numerically  that  dE/ds 
=  0,  as  promised  in  Sec.  II.  A  numerical  module  has  been 
added  to  the  PFB2D  code1  to  solve  the  rms  envelope  equa¬ 
tion,  (36)  which  determines  the  rms  beam  radius  given  the 
periodic  solcnoidal  magnetic  field  and  beam  perveanee,  and 
Eqs.  (37)  and  (38),  which  determine  the  beam  density  and 
scalar  potential  for  the  self-clectric  field. 

We  consider  a  thermal  beam  focused  by  a  periodical 
solcnoidal  foeusing  magnetic  held  defined  by  the  ideal  peri¬ 
odic  step  function  kz(s)  =  kz(s+S)  with 


|  \  kz 0  =  const,  -  i}/2  <  s/S  <  77/2, 
[0,  ?//2  <  s/S  <  1  -  77/2, 


(39) 


where  77  is  the  filling  factor  of  the  solenoid  magnetic  field.  In 
Fig.  1,  the  profile  of  the  normalized  axial  magnetic  field 
S\  a\(.v)  is  plotted  as  a  dashed  curve,  and  the  normalized  rms 
beam  envelope  PbTm^  =  rbrms/  v4elhS  *or  the  thermal  beam  is 
plotted  as  a  solid  curve.  The  system  parameters  arc 
S \'Kj(0)  =2. 1 2,  77=0.4,  £=5A74eth=10,  and  cob=0.  The 
vacuum  and  spaee-eharge-depressed  phase  advances  of  the 
particle  betatron  oscillations  over  one  lattice  period  arc 
evaluated  to  be  <r()  =  e  TJ  yds  /  2  rbrmjl0 =78.9°  and  a 


-£TJflds/2r^Tms=  10.7°,  respectively.  Here,  rbrms0  is  the  rms 
equilibrium  beam  envelope  when  K- 0. 

In  Fig.  2,  the  beam  density  relative  to  the  peak  density 
n{r,s)/ n(0,s)  is  plotted  as  a  function  of  the  radius  relative  to 
the  Debye  length  at  ,y=0  for  the  same  beam  as  in  Eig.  1. 
Here,  the  Debye  length  is  defined  as  X/} 
=  \  7^A:aT1(0)/47ri72/i/;(0,0).  The  density  has  a  Hat  top  near 
the  center  of  the  beam  and  drops  to  zero  within  a  few  Debye 
lengths  near  the  edge  of  the  beam. 

In  Fig.  3,  thermal  beam  density  profiles  arc  plotted  lor 
£=0.1,  1.0,  3.0,  and  10.0  with  the  focusing  field  in  Eq.  (39), 
xz(0)  =  2. 12,  and  no  beam  rotation  in  the  Larmor  frame 
(i.e.,  o^=0).  Here,  the  beam  density  is  normalized  to  the 
peak  density  n0  of  the  beam  with  SA74«lh=10.  The  beam 
density  becomes  a  flat  profile  near  the  beam  axis  as  the  nor¬ 
malized  perveanee  K  =  SK/4Elh  increases,  i.c.,  as  the  beam 
current  increases  or  the  temperature  decreases. 


FIG.  3.  Plot  of  the  relative  beam  density  versus  rf  \4eUl5  for  several  beams 
with  5A74elh=0.l,  1.0,  3.0.  and  10.0,  and  other  system  parameters  the  same 
as  in  Fig.  1.  Here,  the  beam  densities  arc  normalized  to  the  peak  density  of 
the  beam  with  SK/4nlh=  10.0. 
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FIG.  4,  Plol  of  the  on-axis  sell-electric  potential  relative  to  the  beam  trans¬ 
verse  thermal  energy  as  a  function  of  the  propagation  distance  s/S  for  sev¬ 
eral  beams  with  SKI 4e,h=0,l,  1.0,  3.0  and  10.0,  and  other  system  param¬ 
eters  the  same  as  in  Fig.  I. 


The  on-axis  self-clectric  potential  is  determined  numeri¬ 
cally,  requiring  the  total  number  of  partieles  per  unit  length 
to  be  constant.  For  the  detailed  numerical  method,  please 
refer  to  Ref.  <S.  In  Fig.  4,  the  computed  on-axis  self-clcctric 
potential  energy  relative  to  the  beam  transverse  thermal  en¬ 
ergy,  q(/?clf(Q ,s)/ yf}kBT  (5),  is  plotted  as  a  function  of  s/S 
for  5A’/4eth=0.1,  1.0,  3.0,  and  10.0,  and  other  system  param¬ 
eters  the  same  as  in  Fig.  1.  The  variation  of  the  on-axis 
self-clcctric  potential,  i.e.,  the  axial  electric  field,  is  indeed 
small. 

To  illustrate  the  inllucncc  of  the  beam  rotation  rate  in  the 
Larmor  frame  on  the  periodically  focused  thermal  beam 
equilibrium,  we  plot  the  relative  beam  density  profiles  for 
three  choices  of  the  rotation  parameter:  o>/,  =  0,  0.9,  and  0.99 
in  Fig.  5.  The  system  parameters  are  S\/V(0)=2.12  and  K 
=  10.  As  the  beam  rotation  increases,  the  beam  radius  in- 


F1G.  5.  Plot  of  the  relative  beam  denMly  versus  r!  \4elh5  for  several  beams 
with  SK74clh  =  10,0,  oj/,=0,  0.9.  and  0,99,  and  other  system  parameters  the 
same  as  in  Fig.  1.  Here,  the  beam  densities  are  normalized  to  the  peak 
density  of  the  beam  with  a^=0. 


FIG.  6.  Plot  of  the  quantity  (SI2ir(E))cIE/(is  versus  s/S  for  the  four  radial 
displacements  of  the  beam  with  system  parameters  S\ k^-2.  1 2,  77=  0.4 
(170=80°),  SK/4*lh=I.O,  and  tah= 0. 

creases,  and  the  peak  density  at  the  beam  axis  decreases. 
Flowever.  the  Debye  length  is  intact  as  the  beam  rotation  rate 
varies. 

Finally,  we  demonstrate  the  approximate  invariant  of  the 
sealed  transverse  Hamiltonian  as  detined  in  Fq.  (25)  for  the 
cases  SK/g7 — >1.  Instead  of  showing  dE/ds**  0  for  each 
single  particle,  which  requires  very  intensive  numerical  cal¬ 
culations,  we  demonstrate,  by  numerical  calculations,  that 
the  scaled  transverse  Hamiltonian  E  is  slowly  varying  at  a 
few  radial  displacements.  In  Fig.  6,  the  quantity 
(Sl27r(E))dEl ds  is  plotted  as  a  function_of  .s  for  various 
radial  displacements  r=().5rbrms,  rbrms,  \2;brms.  and  2rbrms 
with  other  system  parameters  (70=80o,  K=  1,  and  oj/;=0. 
I  lerc,  (E)  =  n~h 1 JJ  EfbdPxdPy =\(r2  +  £T/2)  +  vr2(.v)  </>sel,(r .  s) 
+  [ATv,4(5)/4rj;rms(.v)]r2  is  the  scaled  transverse  Hamiltonian 
averaged  over  the  particles  located  at  the  same  radial  dis¬ 
placement.  Indeed,  a  maximum  value  of  \(S/27r(E))dElds\ 
=0.06,  which  is  achieved  at  s/S= 0.2  and  0.8,  assures  that 
dE/ds  =  0  in  the  paraxial  approximation. 

V.  COMPARISON  BETWEEN  THEORY 
AND  EXPERIMENT 

In  this  section,  we  compare  the  present  adiabatic  thermal 
beam  equilibrium  theory  with  experimental  measurements. 
The  system  is  a  5  keV  electron  beam  focused  by  a  short 
solenoid  magnetic  in  one  of  the  experiments  of  UMER  In 
Ref.  5,  the  electron  beam  was  generated  by  a  gridded  gun 
and  exited  the  gun  through  an  anode  aperture  at  5=0.  Bell¬ 
shaped  beam  density  profiles  were  imaged  by  a  fluorescent 
screen  while  the  detailed  velocity  space  distribution  was  not 
accessible.  The  bell-shaped  beam  density  protile  and  the 
change  of  the  beam  density  shape  as  the  beam  propagates 
have  not  been  well  understood  theoretically  using  previous 
equilibrium  theories,  such  as  the  KV  beam  equilibrium 

Using  our  adiabatic  thermal  beam  equilibrium  theory,  wc 
replicate  the  beam  density  profiles  at  different  axial  distances 
in  good  agreement  with  the  experimental  measurements.  Our 
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PIG.  7.  Plots  of  calculated  heam  radius  \'2 rhtim  (solid  curves)  for  three 
5  kcV  electron  beams  with  currents  (omittances):  43  mA 
=  71  mm-mrad),  6.5  mA  (4eITT11<=3[)  mm-mrad),  and  1.9  mA  (4einTls 
=  20  mm-mrad).  The  dotted  curves  are  the  availahle  experimental  measure¬ 
ment  for  two  beams:  43  mA  (4eiVlT1.  =  7l  mm-mrad)  and  1.9  mA  {4cItit,s 
=  20  mm  mrad).  The  on-axis  magnetic  field  is  shown  as  a  dashed  curve. 


equilibrium  theory  is  applicable  to  this  experiment  from  the 
anode  aperture  to  a  distance  prior  to  wave  breaking  initiated 
by  high-order  density  distribution  fluctuations  induced  by  a 
pressure  force  at  the  anode  aperture.  Wave  breaking  4  occurs 
at  about  one-quarter  of  plasma  wavelength,  which  is  about 
30  cm  in  this  example.  Our  equilibrium  theory  cannot  ex¬ 
plain  the  density  distribution  distortion  in  the  present  form, 
but  it  will  be  possible  to  develop  a  perturbation  theory  based 
on  the  equilibrium  in  the  future. 

The  calculated  rms  beam  radii  by  solving  Eq.  (30)  are 
shown  to  agree  with  the  available  experimental  rms  beam 
radius_ measurements.  In  Fig.  7,  the  calculated  beam  radii 
R  —  \  2/*brms  are  plotted  as  solid  curves  by  solving  Eq.  (30) 
for  three  5  keV  electron  beams  with  currents  (emit- 
tanees).  43  mA  (4*-:A-Tms=7 1  mm-mrad),  6.5  mA  (4e?rrm 
=30  mm-mrad),  and  1.9  mA  20  mm-mrad).  The 

three  beams  are  focused  by  a  short  solenoid  magnet  whose 
on-axis  magnetic  field  is  shown  as  a  dashed  curve.  The  cal¬ 
culated  beam  radii  for  the  two  beams  with  currents  43  and 
1 .9  tnA  agree  with  the  previous  experiment  measurements 
(dotted  curves)  and  calculations  in  Ref.  5.  as  expected.  The 
calculated  beam  radius  for  the  6.5  mA  beam  will  be  used  for 
the  following  density  ealeulalions. 

By  solving  Eqs.  (37)  and  (38),  we  calculate  the  beam 
transverse  density  profiles  of  ihe  UMER  5  keV,  6.5  mA 
electron  beam  at  three  axial  distances:  .v=6.4,  11.2,  and 
17.2  cm,  as  shown  by  the  solid  curves  in  Fig.  8.  The  dashed 
curves  are  the  equivalent  KV  beam  density  profiles.5"  Com¬ 
pared  with  the  experimental  measurements  (dotted  curves)," 
the  calculated  beam  density  profiles  are  in  good  agreement. 
As  the  beam  radius  increases,  the  beam  density  profile  ap¬ 
proaches  the  KV  (uniform)  beam  density  profile,  because  the 
beam  temperature  must  decrease  in  order  to  keep 
T  1(s)rlrms{s)  at  a  constant.  In  this  adiabatic  process,  the  De¬ 
bye  length  \p—  \  %kBT  (,v)/4m/2/^(0,.v)=0.54  mm  is  con¬ 
stant. 


FIG.  8.  Normalized  transverse  beam  densily  profiles  of  a  5  keV,  6.5  mA 
(4cfrms=30  mm-mrad)  electron  beam  at  three  axial  distances:  =  6.4.  1 1.2, 
and  17.2  cm.  The  solid  curves  are  from  theory,  lhe  doited  curves  are  the 
expenmental  measurements,  and  the  dashed  lines  are  lhe  equivalent  KV 
beam  density  dislribulions.  The  densilies  are  normalized  to  the  equivalent 
KV  beam  density  at  ,r=6.4  cm. 

VI.  CONCLUSION 

In  conclusion,  an  adiabatic  thermal  equilibrium  was  dis¬ 
covered  for  an  intense,  axisymmetric  charged-particle  beam 
propagating  through  a  periodic  solenoidal  focusing  field.  The 


Downloaded  08  Feb  2008  to  198.125  179  144.  Redistribution  subject  to  AIR  license  or  copyright:  see  http://pop.aip  org/pop/copyrighi.jsp 


023102-8  Zhou,  Samokhvalova,  and  Chen 


Phys.  Plasmas  15,  023102  (2008 ) 


thermal  beam  distribution  function  was  constructed.  The 
beam  rms  envelope  equation  was  derived,  and  the  self- 
consistent  nonuniform  density  profile  was  calculated.  Other 
statistic  properties  such  as  flow  velocity,  temperature,  total 
emittance  and  rms  thermal  cmittance,  equation  of  state,  and 
Debye  length  were  studied.  Good  agreement  was  found  be¬ 
tween  the  adiabatic  thermal  equilibrium  theory  and  recent 
high-intensily  beam  experimental  measurements  at  the  Uni¬ 
versity  of  Mary  land  Electron  Ring. 
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A  small-signal  theory  of  a  nonrelativistic  magnetron  is  developed  using  a  planar  model  with  a  thin 
electron  cloud.  The  theory  includes  both  inertial  effects  and  electromagnetic  effects  in  a  Floquet 
expansion.  An  analytical  dispersion  relation  of  such  a  planar  magnetron  is  derived,  and  Ihe  growth 
rate  is  calculated  analytically.  The  found  instability  involves  the  resonance  between  the  electron 
cloud  and  the  magnetron  cavity  vacuum  slow  waves.  A  good  agreement  is  found  between  the  theory 
and  the  self-consistent  particle-in-cell  MAGIC  simulations.  ©  2008  American  Institute  of  Physics. 
[DO  I:  10.1063/1.2996577] 


Crossed-lield  devices,  such  as  magnetrons  and  erossed- 
field  amplifiers  (CFA),  are  of  great  interest  because  of  their 
low  cost,  high  efficiency,  compactness,  and  robustness  They 
have  many  civilian,  industrial,  and  military  applications.1 
Such  devices  generate  high  power  rf  radiation  via  stimulated 
interactions,  as  high  density  electron  clouds  drift  across  the 
crossed  electric  and  magnetic  fields  in  periodic  cavities.  On 
the  qualitative  level,  the  generation  and  operation  of  such 
high  density  electron  clouds  are  relatively  straightforward 
following  electrostatics  and  magnetostatics.  The  mechanisms 
of  the  rf  wave  generation,  however,  have  yet  to  be  fully 
described  due  to  the  complicated  processes  of  the  rf  interac¬ 
tion. 

The  onset  of  (he  unstable  oscillations  in  magnetrons 
has  not  been  analytically  described  although  extensive 
particle-in-cell  (PIC)  simulations  can  make  good  pre¬ 
dictions  for  the  instability  characteristics.  Previous  studies 
included  various  models  utilizing  linear  theories.  Earlier 
work  focused  on  the  diocotron  instability  in  the  guiding- 
center  approximation,  which  ignores  inertial  effects  in  the 
electron  cloud.  The  recent  work  by  Riyopoulos,  using  a 
guiding-center  model,  provided  insight  into  the  magnetron 
instability  in  the  low-space-charge  limit.*  A  linear  theory  tak¬ 
ing  a  single  rf  mode  in  the  Floquet  expansion  was  developed 
to  include  electromagnetic  effects.9  Despite  these  theoretical 
and  PIC  simulation  efforts,  a  quantitative  agreement  between 
theory  and  PIC  simulations*  has  not  been  reported  until  this 
letter. 

In  this*  letter,  we  develop  a  small-signal  (linear)  theory 
that  includes  both  inertial  and  electromagnetic  effects  in  a 
Floquet  expansion.  We  consider  a  planar  magnetron  or  CFA 
with  the  geometry  as  shown  in  Fig.  I.  The  A-K  gap  of  the 
cylindrical  magnetron  is  assumed  to  be  small  compared  with 
the  cathode  radius.  As  shown  in  Fig.  1 ,  the  cathode  is  located 
at  x-0  and  the  anode  is  a  slow-wave  structure  consisting  of 
a  periodic  array  of  vanes  and  slots.  The  applied  magnetic 
field  points  out  of  the  page,  which  forces  the  electron  drift 
velocity  in  the  y-direction.  In  a  conventional  magnetron  or 
CFA  there  is  an  electron  sheath  that  extends  out  for  some 
distance  from  the  cathode.  However,  such  a  finite  electron 
sheath  introduces  a  problem  of  multiple  poles  (singularities) 
in  the  theoretical  computation  of  the  small-signal  gain.  To 
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avoid  the  singularity  problem,  we  approximate  the  electron 
sheath  as  an  infinite  thin  electron  cloud  stream,  which  is* 
located  at  x-h  Because  of  the  resonant  interaction,  our  thin 
electron  cloud  stream  model  captures  the  basic  gain  mecha¬ 
nism  in  the  realistic  shear  (low  in  a  magnetron. 

In  the  small-signal  theory,  we  assume  that  (he  system 
starts  from  an  equilibrium  (dc)  state  and  a  small  oscillation 
(rf)  builds  up  as  the  electron  cloud  drifts  in  the  A-K  gap.  The 
general  form  of  any  field  variable  <//  is 

(//=  i Mx.y)  +  Sif^x.y.t).  ( I ) 

where  t/'o U,y)  denotes  the  equilibrium  field  variable.  Be¬ 
cause  the  magnetron  is  periodic  in  the  y-direction  with  peri¬ 
odicity  L,  all  of  the  perturbations  in  Eq.  (I)  must  be  of  the 
Floquet  form, 

X 

<5<A (x,y,t)=  S  (2) 

<7=-x  q 

where  k(=kz+2irql L.  The  small-signal  theory  is  valid  for 
the  initial  growth  stage  of  the  magnetron  interaction.  Near  or 
after  saturation,  a  nonlinear  theory  is  needed.1 

The  corrugation  on  the  anode  induces  a  small  static  pe¬ 
riodic  perturbation  on  the  equilibrium  profile  of  the  electron 
sheath,  as  observed  in  our  two-dimensional  (2D)  magic 
simulations.  In  our  simplified  model,  we  approximate  the 
thin  electron  sheath  by  a  number  of  electron  cloud  streams 
with  slightly  different  flow  velocities  and  express  the  total 
equilibrium  charge  and  current  densities  as  p(x)  = 


FIG  I.  Schematic  of  a  planar  magnetron  with  a  thin  electron  cloud  stream. 
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-eljl,(TbjS(x-li)  and  J (x)=-e’ZjLla-bjVhjS(x-h)ex,  respec- 
tively,  where  M  is  the  number  of  electron  cloud  streams,  -e 
is  the  electron  charge,  crb=^=lrrbj  is  the  total  surface  number 
density  of  the  electron  cloud,  and  x=h  specifies  the  trans¬ 
verse  displacement  of  the  electron  sheath.  The  equilibrium 
How  velocity  Vbj  is  determined  by  adding  a  small  velocity 
A Vhj  from  the  averaged  flow  velocity, 


Yr 

+  eJIt( 

_i+*V 

_  b 

2  b)  _ 

(3) 


where  VD  is  the  applied  A-K  gap  voltage,  B0  is  the  applied 
magnetic  field  strength,  /-;()  is  the  permittivity  in  vacuum,  and 
b  is  the  A-K  gap  width,  as  shown  in  Fig.  1.  Typically,  A Vbj  is 

a  few  percentage  of  Vh,  whose  value  can  be  estimated  from 
the  2D  MAGIC  simulations. 

Using  the  linearized  nonrelativistic  cold-fluid  equations 
and  the  linearized  Maxwell  equations,  the  eigenvalue  equa¬ 
tion  for  SEvL,  is  derived  to  be 
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where  t oc-eB0fm  is  the  electron  cyclotron  frequency, 
(o^i(.x)  =  e2cThjS(x-li)l e0/n  is  the  effective  plasma  frequency 
squared  of  the  jth  electron  cloud  stream,  and 
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and  hTy ( o) , kq)  is  the  element  of  the  inverse  matrix 
H ~l(o),kq).  Here,  m  is  the  electron  mass  and  c  is  the  speed  of 
light  in  vacuum. 

The  boundary  conditions  for  the  electric  field  SEyk  are 
that  the  electric  held  vanishes  at  the  cathode  jc=0,  and  the 
anode  wall  x=cl  and  is  continuous  in  the  y-direction  at  the 
electron  cloud  layer  .v=6;  the  admittance1  is  continuous  at 
the  vane  tip  x-b  and 
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Equation  (6)  can  be  derived  by  integrating  Eq.  (4)  for  x=li 
-0  to  x  =/;  +  ()  and  making  use  of  fT\{o),kq) 

Solving  Eq.  (4)  in  the  regions  O^x^h,  h^x^b,  and 
b^x^cl  and  matching  the  solutions  at  the  boundaries,  we 
arrive  at  the  small-signal  gain  equation  or  the  loaded  disper¬ 
sion  relation 


FIG.  2.  Plot  of  the  maximum  temporal  growth  rate  Imfco)  as  a  function  of 
the  surface  density  crb  for  a  single  ihin  electron  cloud  stream  with  the  struc¬ 
ture  parameters  of  C=0.478  cm,  a =0.382  cm.  6=0.478  cm.  and  <7 
=4.25  cm  in  a  uniform  applied  magnetic  field  of  B  =180  G  and  a  A-K  gap 
voltage  of  Vp=2.52  kV 
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and  pq- \(o2/c2-k~r  Setting  rk  =0,  Eq.  (7)  is  the  vacuum 
dispersion  relation  for  the  corrugated  structure.  A  numerical 
code,  named  crossed-lield  amplilier  system  simulator 
(cfass),  is  developed  to  solve  Eq.  (7). 

To  show  the  dependence  of  the  growth  rate  of  an  un¬ 
stable  mode  on  the  electron  density,  we  plot  the  maximum 
temporal  growth  rate  as  a  function  of  the  surface  density  ah 
in  Fig.  2.  The  parameters  of  the  structure  are  chosen  to  be 
L=0.478  cm,  <7  =  0.382  cm,  6=0.478  cm,  and  d- 4.25  cm  in 
a  uniform  applied  magnetic  field  of  Bz=  180  G  and  an  A-K 
gap  voltage  of  VD=2.52  kV.  The  electron  cloud  is  assumed 
to  be  an  infinitely  long  sheet  with  the  surface  density  ah> 
which  propagates  at  the  velocity  given  by  Eq.  (3),  i.e..  the 
ExB  drift  velocity,  at  the  location  of  6=0.382  cm.  As  the 
surface  density  is  varied,  the  phase  shift  of  the  unstable  mode 
is  tuned  to  achieve  the  maximum  temporal  growth  rate.  As 
shown  in  Fig.  2,  the  maximum  growth  rate  increases  as  the 
surface  density  increases. 

To  compare  our  small-signal  gain  theory  with  2D  MAGIC 
simulations,  we  restrict  our  discussion  to  the  low-current  re¬ 
gime,  where  the  thin-beam  equilibrium  model  is  a  good  ap¬ 
proximation  to  the  thin  electron  cloud  in  the  2D  MAGIC  simu¬ 
lations.  At  high  currents,  the  electron  cloud  diffuses  due  to 
the  strong  image  effects  on  the  corrugated  anode,  and  an 
improved  model  of  electron  cloud  equilibrium  is  needed, 
which  is  beyond  the  scope  of  this  letter. 

As  an  example,  we  consider  the  same  structure  as  in  Fig. 
2.  The  electron  cloud  is  assumed  to  be  an  infinitely  long 
sheet  with  crh= 2.1 2  X  108  cm-2,  which  propagates  at  the  av¬ 
eraged  velocity  V/,=0.098r  at  6=0.382  cm.  The  lowest  un- 
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I  IG.  3.  Plot  of  the  amplitude  of  oscillating  voltage  as  a  function  of  time 
as  obtained  from  the  magic  simulation.  Here,  the  parameters  arc  L 
=0.478  cm,  «=0.382  cm,  ft=0.478  cm,  d- 4.25  cm,  h =0.382  cm.  B 
=  180  G,  £>-5.27  kV/cm.  cr*=2.l2X  10K  cm'2,  and  V*«0.098c. 

stable  transverse  mode  (TM)  is  calculated  to  achieve  the 
maximum  growth  rate  at  a  phase  of  2 tt/ 3  and  a  frequency  of 
1 .979  GHz. 

The  self-consistent  PIC  code,  2D  magic,  is  used  to 
simulate  the  planar  magnetron  system.  Because  2D  magic 
can  handle  only  a  few  vanes,  a  three-vane  slow-wave  corru¬ 
gated  structure  with  the  same  parameters  as  in  Fig.  2  is  used 
in  the  simulation.  Periodic  boundary  conditions  are  used 
such  that  the  27t/3  mode  is  supported  by  the  three- vane 
structure.  The  uniform  crossed  electric  and  magnetic  fields 
arc  applied  with  B:=  180  G  and  £>-5.27  kV/cm.  The  elec¬ 
tron  beam  is  initialized  as  an  infinitely  long  slab  in  the  z 
direction  with  a  width  of  0.05  cm  in  the  ,r  direction.  The 
electron  beam  propagates  with  an  initial  velocity  of  Vh 
= 0.098c. 

As  the  beam  propagates,  the  instability  starts  to  build  up, 
which  is  illustrated  by  the  voltage  crossing  the  vane  tip  of 
the  slow-wave  structure.  In  Fig.  3,  the  amplitude  of  the  os¬ 
cillating  voltage  filtered  by  a  filter  that  selects  the  2ir/?> 
inode  is  plotted.  In  the  early  stage  of  the  instability  (e.g.,  for 
/  <  1 5  ns),  the  oscillation  has  a  very  small  amplitude  and 
exhibits  a  relatively  broad  frequency  spectrum.  Starting  from 
f  =  10  ns,  the  27r/3  mode  grows  exponentially.  It  saturates  at 
about  /  =  26  ns.  For  this  MAGIC  simulation,  the  2tt/3  mode  is 
determined  to  have  a  frequency  of  1.93  GHz  and  an  ampli¬ 
tude  growth  rate  of  3.12  dB/cm,  shown  as  a  circle  and  a 
cross  in  Fig.  4,  respectively. 

To  compare  the  MAGIC  simulation  results  with  our 
theory,  we  solve  the  loaded  dispersion  relation  in  Eq.  (7)  to 
calculate  the  real  frequency  and  the  instability  growth  rate. 
Due  to  the  influence  of  the  anode  corrugation,  the  electron 
cloud  velocity  prior  to  the  linear  growth  is  observed  to  vary 
sinusoidally  in  the  v-direction  with  a  small  amplitude  around 
the  averaged  flow  velocity  in  the  MAGIC  simulation.  To 
model  the  velocity  variation,  we  use  three  electron  cloud 
streams  at  the  same  location  .v=0.0384  cm.  each  with  one- 
third  of  the  total  surface  density  ab ( =  crh2 = = °7>/ 3  an^ 
slightly  different  velocities,  i.e.,  ^>0.98^,  Vb2=Vb,  and 


EIG.  4.  Plot  of  the  temporal  growth  rate  Im(oj)  (solid  curve)  and  the  real 
frequency  Rc(tu)  (dashed  curve)  as  a  function  of  wave  number  k.  for  the 
lowest  resonant  TM  mode  The  solid  and  dashed  curves  are  obtained  from 
Eq.  (7)  using  the  three  electron  cloud  streams  with  surface  densities  ahi 
=  (Th2=crbi=  (j/,/3  =  7.07  X  107  cm"2  and  flow  velocities  V^j^O.OK  >  Vh2 
=  Vh,  =0.098c,  and  V*>  1.02V*,  respectively.  The  circle  and  cross  arc  the 
real  frequency  and  instability  grow  rate  from  the  magic  simulation,  respec¬ 
tively.  Here,  the  other  parameters  are  the  same  as  those  in  Fig.  3. 

1^3=1.021^.  The  theoretical  growth  rate  is  in  good  jgree- 
ment  with  the  magic  simulation,  as  shown  in  Fig.  4. 

In  conclusion,  a  small-signal  theory  of  a  nonrelativislic 
magnetron  was  developed  using  a  planar  model  with  a  thin 
electron  cloud.  The  theory  includes  inertial  effects  and  elec¬ 
tromagnetic  effects  in  a  Floquet  expansion  The  present  pla¬ 
nar  model  is  valid  provided  that  the  A-K  gap  is  small  com¬ 
pared  with  the  cathode  radius.  The  thin  electron  cloud  avoids 
the  problem  of  multiple  poles  (singularities)  in  the  dispersion 
relation.  An  analytical  dispersion  relation  of  such  a  planar 
magnetron  was  derived,  and  the  growth  rate  was  calculated 
analytically.  2D  MAGIC  simulations  were  performed  to  verify 
the  theory.  A  good  agreement  was  found  between  the  theory 
and  magic  simulations. 
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Abstract 

A  self-consistent  solution  is  developed  for  the  focusing  of  coasting,  sheet-like,  space- 
charge-dominated  elliptic  beams  using  a  hybrid  of  non-axisymmetric  periodic  permanent 
magnets  (PPM)  and  quadrupole  magnets.  The  beam  envelope  equations  and  equilibrium 
flow  profiles  are  obtained  using  a  paraxial  cold-fluid  model.  Further,  a  kinetic  theory  is 
developed  to  determine  the  self-consistent  beam  distribution  and  its  evolution  and  to 
study  the  temperature  effects.  Both  cold-fluid  and  kinetic  equilibrium  theories  are 
validated  by  two-dimensional  (2D)  particle-in-cell  (PIC)  simulations  using  the  2D 
Periodically  Focused  Beam  (PFB2D)  code,  and  by  three-dimensional  (3D)  self-consistent 
trajectory  simulations  using  OMNITRAK.  The  theories  and  PFB2D  and  OMNITRAK 
simulations  are  applied  to  design  such  high-space-charge  elliptic  beams  for  applications 
in  klystrons  and  vacuum  electron  devices.  Numerical  results  show  that  the  beam  edges  in 
both  transverse  directions  are  well  confined  without  twisting  and  the  elliptic  beam  density 
profile  is  well  preserved.  For  space-charge-dominated  elliptic  beams  used  in  ribbon-beam 


klystron  applications,  the  temperature  effects  are  studied  using  both  the  kinetic 
equilibrium  theory  and  PFB2D  simulations. 


PACS:  29.27. Bd,  52.59.Sa,  47.75.+f,  52.25.Dg 
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I.  Introduction 


There  are  vigorous  activities  in  the  research  and  development  of  elliptic-beam  sources, 
traveling  wave  amplifiers,  klystrons,  and  accelerator  focusing  systems.  Over  600  high- 
power,  high-efficiency  klystrons,  for  example,  may  be  needed  to  provide  rf  power  for  the 
acceleration  cavities  of  the  proposed  TeV  International  Linear  Collider  (ILC).  The 
Stanford  Linear  Accelerator  Center  (SLAC)  has  proposed  a  10  MW  sheet-beam  klystron 
to  meet  this  need  [lj.  Other  groups,  such  as  Los  Alamos  National  Laboratory  (LANL) 
and  Scientific  Applications  International  Corporation  (SAIC),  are  also  interested  in  sheet- 
beam  technology  for  microwave  amplifier  applications  [2,  3].  The  LANL  team  has 
employed  a  solenoid/quadrupole  magnet  combination  to  transform  an  incident  circular 
beam  into  an  emergent  elliptic  beam  [4].  The  Massachusetts  Institute  of  Technology 
(MIT)  has  also  initiated  an  elliptic  beam  project  for  communications  and  accelerator 
applications  [5].  Beam  Power  Technology,  Inc.  is  developing  a  ribbon-beam  amplifier 
based  on  the  elliptic  beam  technology  developed  at  MIT. 

Elliptic  beams  have  attracted  broad  interests  because  they  have  the  following 
remarkable  properties.  First,  they  can  transport  large  amounts  of  beam  currents  at  reduced 
intrinsic  space-charge  forces  and  energies.  Second,  they  couple  efficiently  to  rectangular 
or  elliptic  rf  structures.  The  combination  of  the  space-charge  reduction  and  efficient 
coupling  allows  efficient  rf  generation  in  vacuum  electronic  devices,  and  efficient 
acceleration  in  particle  accelerators.  Third,  elliptic  beams  provide  an  additional  adjustable 
parameter  (e.g.,  the  aspect  ratio)  which  may  be  useful  for  better  matching  a  beam  into  a 
non-axisymmetric  periodic  focusing  channel. 
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Transport  and  focusing  of  sheet  (ribbon)  beams  has  been  discussed  in  the  literature  for 
four  decades.  Periodic  transverse  (wiggler)  magnetic  focusing  [6,  7]  has  been  used  for 
free-electron  laser  applications,  but  it  can  lead  to  excessive  centroid  motion  for  space- 
charge-dominated  beams  [8].  Promising  results  have  been  obtained  through  recent  studies 
of  period-averaged  focusing  in  the  periodic  cusp  magnetic  (PCM)  field  [8,  9]  or  the 
hybrid  of  PCM  and  the  periodic  quadrupole  magnet  (PQM)  field  [10,  11]  for  space- 
charge-dominated  and  emittance-dominated  sheet  beams,  however,  significant  envelope 
oscillations  and  emittance  growth  are  sometimes  seen.  They  may  be  rectified  by  the 
present  thorough  treatment  which  does  not  employ  period-averaging  but  self-consistently 
includes  the  effects  of  beam  flow  profiles  and  evolving  self-fields  neglected  in  the 
period-averaging  approximation. 

Recent  efforts  in  this  vein  [5,  12]  have  led  to  a  deeper  understanding  of  high-space- 
charge  elliptic  beam  propagation  in  a  non-axisymmetric  periodic  permanent  magnet 
(PPM)  field.  Pure  non-axisymmetric  PPM  focusing,  however,  is  unsuited  for  sheet-like 
elliptic  beams  with  very  large  aspect-ratio  because  the  twisting  of  the  beam  introduces 
instabilities,  and  because  the  magnetic  field  nonlinearities  in  the  wide  tails  of  the  beam 
become  appreciable. 

In  this  paper,  we  develop  a  self-consistent  solution  for  the  focusing  of  coasting,  sheet¬ 
like,  space-charge-dominated  elliptic  beams  using  the  most  general  formulation  of 
centroid  preserving  linear  fields  -  a  hybrid  of  non-axisymmetric  PPM  and  quadrupole 
magnets  [13].  A  paraxial  cold-fluid  model  is  employed  to  derive  generalized  envelope 
equations  that  determine  the  equilibrium  flowr  properties  of  ellipse-shaped  beams  with 
negligibly  small  emittance.  Furthermore,  a  kinetic  analysis  is  developed  to  determine  the 
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beam  distribution  consistent  with  the  equilibrium  flow  properties  of  a  uniform  ellipse¬ 
shaped  beam.  In  the  kinetic  model,  the  emittances  are  taken  into  account.  A  matched 
envelope  solution  is  obtained  numerically  from  the  generalized  envelope  equations,  and 
the  results  show  that  the  beam  edges  in  both  transverse  directions  are  well  confined 
without  twisting.  Two  dimensional  (2D)  particle-m-cell  (PIC)  simulations  with  our 
Periodic  Focused  Beam  2D  (PFB2D)  code  and  3D  Omnitrak  simulations  show  good 
agreement  with  the  predictions  of  equilibrium  theory  as  well  as  beam  stability.  For  space- 
charge-dominated  elliptic  beams  used  in  ribbon-beam  klystron  applications,  the 
temperature  effects  are  studied  using  both  the  kinetic  equilibrium  theory  and  PFB2D 
simulations. 

The  organization  of  the  present  paper  is  as  follows.  In  Sec.  II,  the  cold-fluid 
equilibrium  theory  is  used  to  derive  the  generalized  envelope  equations  for  a  high-space- 
charge  elliptic  beam.  In  Sec.  Ill,  the  kinetic  equilibrium  theory  of  a  large-aspect-ratio 
ellipse-shaped  charged-particle  beam  is  presented  following  Sacherer’s  matrix  method 
[14],  and  the  generalized  envelope  equations  which  include  the  emittances  are  derived.  In 
Sec.  IV,  PIC  PFB2D  simulations  are  used  to  verify  the  theoretical  results  and  examples  of 
large-aspect-ratio  elliptic  beams  are  discussed.  Conclusions  are  presented  in  Sec.  V. 
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II.  Cold-fluid  equilibrium  theory 

We  consider  a  high-intensity,  space-charge-dominated  beam,  in  which  kinetic 


(thermal  emittance)  effects  are  negligibly  small,  The  beam  can  be  adequately  described 
by  cold-fluid  equations.  In  the  paraxial  approximation,  the  steady-state  cold-fluid 
equations  for  time-stationary  flow  ( d/dt  =  0)  in  cgs  units  are  [15,  16] 


(1) 


(2) 


where  s-z  ,  \  =  xex +yev ,  V±  =  d/d\x ,  q  and  rn  are  the  particle  charge  and  rest 
mass,  respectively,  nb  is  the  particle  density,  Vx  is  the  transverse  flow 


velocity,  yh  —  (l  -  Pi )  ”  is  the  relativistic  mass  factor,  use  has  been  made  oi 


p.  =  V.  / c  =  Ph  -  const ,  c  is  the  speed  of  light  in  vacuum,  and  the  self-electric  field 
E'  and  self-magnetic  field  B5  are  determined  from  the  scalar  potential  <j)s  and  vector 
potential  Aik. ,  i.e.,  Ev  =  -V  <f>s  and  Bs  =V  ±x  Ale„. 

Wc  use  the  combination  of  a  periodic  non-axisymmetric  magnetic  field  and  a 
quadrupole  magnetic  field  to  focus  a  nearly  straight  large-aspect-ratio  elliptic  beam, 
whose  twist  angle  is  minimized.  For  the  beam  dimensions  small  relative  to  the 

characteristic  scale  of  magnetic  variations,  i.e.,  (k0xx)2  /6  «\  and  (vd2/6  «  1  ,  the 

combined  magnetic  field  can  be  described  to  the  lowest  order  in  the  transverse 
dimensions  as 
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(4) 


B"'  =  B,(s]e.  - 


dB,(s ) 
ds 


xcx  +  — ycy 


+ 


K(s’i}'ex+xe 


o  o  'y 

where  k0  =  2ft/S,  £qv+A'ov=Aq  ,  S  is  the  axial  periodicity  length,  and 


We  seek  solutions  to  Eqs.  ( 1  )-(3)  of  the  form 


ha(x±,s)  = 


Nl 


-© 


y 


m{s)b(s)  a2(.?)  />2(.s) 

vi  (x±  .‘0  =  \Mx (*)*  ~  a* (4v]/^e?  +  [//,  (s)y  +  a  (.v)xj/?4ref . 


(5) 

(6) 


In  Eqs.  (5)  and  (6),  \±  =  xev  +  yc^  is  a  transverse  displacement  in  the  twisted  coordinate 
system  illustrated  in  Fig.  1;  <9(5)  is  the  twist  angle  of  the  ellipse;  ©(.x)  =  1  if  jc  >  0  and 
©(x)  =  0  if  .Y<0;  and  the  functions  a(s),  6(5),  /^(s),  a^),  tfv.(s)  and  0(v) 

are  to  be  determined  self-consistently  [see  Eqs.  (8)-(12)].  We  carry  the  twist  angle  #(.v) 
in  our  calculation  for  the  purpose  of  generalization,  and  we  will  take  #(s)  =  0  later  to 
obtain  the  solution  for  a  nearly  straight  elliptic  beam. 

The  self-electric  and  self-magnetic  fields  are  well  known  for  an  elliptical  beam  with 
density  distribution  specified  in  Eq.  (5)  [17],  i.e., 

/-2  t?2\ 


bs  =  pb  'At  = 

ci  +  b 


x-+y- 

V  a  h  J 


(7) 


Using  the  expressions  in  Eqs.  (4)-(7),  it  can  be  shown  that  both  the  equilibrium 
continuity  equation  (1)  and  force  equation  (3)  are  satisfied  if  the  dynamical  variables 
a(s)  ,  b(s)  ,  jux(.^)  =  a  'da/ds  ,  /jv(s)  =  b  'dbl ds  ,  art(s) ,  ay  (s)  and  B(s)  obey  the 
generalized  beam  envelope  equations  (see  Appendix  A  for  the  derivation) 
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(8) 


d2a  b2[a]-2axa v)+a2a2v  2K  .  , — 7-c 

— r - — r~r 2 - a  -  mv " 2  v^- 

ds  a  -b  (a  +  b) 

kox  ~^0v  d  -J  (•?)  .  /- -a  /  \  /t/i\ 

- = — - — - - sin(2(9)+  /qJsjcos^tf) 


kl 

ds 

2  axav 

)+b~a; 

ds’ 


^(a?-2^.w)+^  b_JW__2J^a  h 

^ 1  a2-b 2  (a  +  b)  V  ‘V 

Ox  oy  sm(2^)+  k-  (s)cos(2fl)  />  =  0, 

/Tg  ds 

a^{ax  ~av)  d  ( a' 


(9) 


d  I  2  \ 

—  a  a.,)-  .  , 

ds '  a^-b"  ds^b 


+  2a’ayficj(s) 


+a‘ 


klx  cos2  6+ kl  sin 2  6  djic.  (5)  ,  ,  .  ,  * 

kl  ds  qX  V 


(10) 


-0, 


—  (&2ar)- 
ds '  x) 


+b- 


a~-b -  ds\a) 

A'()  v  sin 2  <9  +  A'q  v.  cos“  0  d ■Jtci (5) 


kl 


ds 


+  /q.(.s)sin(2$) 


(ID 


=  0, 


dO  a2av-b2ax 


ds 


2  id 
a  -  b 


(12) 


where 


2ybphmc  YbPhmc 


and  A'  .  -2?2'V‘ 


yIpW  ' 


(13) 


Equations  (8)-(  1 2)  have  “time”  reversal  symmetry  under  the  transformation 
{s,a,b,a',b',ax,av,0)  — >  (- s,a,b,-a',- b'-ax—  a ,#).  This  implies  that  the  dynamical 

system  described  by  Eqs.  (8)-(12)  has  the  hyper  symmetry  plane  (a'.b',ax,av)- 
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We  look  for  a  solution  that  has  negligible  twisting,  i.e.,  #(.s)  =  0.  In  such  a  case, 
equation  (12)  leads  to 


arr(s)  a2(s) 
ay{s)  b2{s) 


(14) 


Correspondingly,  the  generalized  envelope  equations  (8)-(l  1)  become 

——5-^  +  xq (5)«(5)-  2yJK~:(shv (5)0(5)- ax(s)ay (5)0(5)-  ,  ,  t  =  0 , 


ds 


0(5)+  b(s) 

dl.^  ~  Kq  (J  2a/*"TK  (#)"  Ux  (s)ay  Ws)  2K 


ds 


0(^9)  +  b(s) 


=  0 , 


dax(s)  |  «'{.?)  |  bjs ) 

ds  o(j)  b(s) 


day(s )  f  a '(5)  A'fc) 
ds  0(5)  b(s ) 


«x(‘O+2V^0ttt  +  27 2-^£0  =  °, 

b{s )  ki  ds 


a(s )  ki  ds 


(15) 


(16) 


(17) 


(18) 


The  cold-fluid  equilibrium  provides  a  set  of  envelope  equations,  i.e.,  Eqs.  ( 1 5)-(  18) 
which  determines  the  semi-axes  of  the  elliptic  beam  and  the  flow  velocity  profile.  The  set 
of  envelope  equations  is  convenient  to  use  for  design  of  an  elliptic  focusing  system. 
Normally,  for  given  periodic  lattice  parameters,  we  solve  Eq.  (15)-(18)  to  obtain  the 
matched  solution  of  beam  envelopes  using  periodic  boundary  conditions,  which  is  a 
forward  problem.  On  the  other  hand,  for  given  beam  envelope  dimensions,  we  can  solve 
Eq.  (15)-(18)  to  determine  the  required  lattice  parameters,  which  is  an  inverse  problem. 
As  will  be  discussed  in  Sec.  IV,  we  use  the  set  of  envelope  equations  to  design  two 
periodic  lattice  systems  for  two  examples  of  elliptic  beams  given  the  beam  parameters  in 
Table  1.  However,  the  cold-fluid  theory  does  not  provide  detailed  information  of  the 
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beam  distribution  and  it  dose  not  take  the  emittance  and  temperature  effects  into  account, 
which  might  be  important  in  practical  devices.  The  kinetic  treatment  required  to  address 
the  emittance  and  temperature  effects  is  discussed  in  Sec.  III. 


10 


III.  Kinetic  Equilibrium  Theory 


In  this  section,  we  develop  the  kinetic  equilibrium  of  the  periodic  focused  large-aspcct- 
ratio  elliptic  beams  when  the  emittance  is  no  longer  small  and  the  cold-fluid  theory 
presented  in  Sec.  11  is  no  longer  valid.  Following  Sacherer’s  theory  [14,  18]  that  there 
exits  a  self-consistent  uniformly  charged  beam  in  situations  where  both  external  forces 
and  self  forces  acting  on  the  charged  particles  are  linear,  we  can  construct  the  beam 
distribution  and  find  the  self-consistent  envelope  equations  for  the  applied  magnetic  field 
configuration  in  Eq.  (4).  For  simplicity,  we  consider  the  twist  angle  0  =  0  in  the 
following  calculations.  The  results  recover  the  cold-fluid  theory'  in  the  cold-fluid  limit. 

We  express  the  transverse  single-particle  equations  of  motion  in  the  paraxial 
approximation, 


(19) 


p  =  my — —  s  m 
dt 


and 


(20) 


where  p  is  the  relativistic  particle  momentum.  Here,  E  =  -V^T  with  <f>s  given  in  Eq.  (7) 

and  B  =  Bs  +  Bex  with  B  =  Vx  x  As:c:  and  B^v/  given  in  Eq.  (4). 

Since  all  the  terms  in  Eqs.  (19)  and  (20)  are  linear,  we  reorganize  the  equations  of 
motion  (19)  and  (20)  into  a  matrix  format, 


(21) 


where 


x  = 


and 


(0 


F  = 


y 

\y  j 


FXf>x  0  0 


(22) 


F  x  0  F  n  v  Fn  „ 

PxX  Pxy  PxPy 

0  0  0 


ypy 


F  F  F  0 

,  Py *  PyPx  pyy 


(23) 


The  prime  denotes  the  derivative  with  respect  to  s  .  For  6  -  0  >  the  elements  of  F  are 


F.xpx  FyPy  1 


2  K 


PxX  a(s)[a(s)  +  b(s)\  ^ 


_  0  kly  dFz(s) 


F  =2- 
P'y  kt  ds 


FPxPy  =2FF)> 


^klx  d^Kz{s) 
PyX  ~  4  ds 


FPvPx  =~2FF), 


2  K 


j  +  *V(T 


(24) 

(25) 

(26) 

(27) 

(28) 

(29) 

(30) 


For  a  single  charged  particle,  there  exists  a  transfer  matrix  T(s,s0)  that  transfonns  an 
initial  location  Xo^o)  'nto  a  corresponding  location  at  a  later  time  x(s),  i.e., 

x(^)=T(-s>5o)xo(5o)  (31) 
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and  note  that  the  elements  of  T(s,s0)  satisfy 

^fed  =  F(s).T(„0).  (32) 

as 

Since  we  know'  the  equations  of  motion  for  individual  charged  particles,  vve  can 
determine  the  evolution  of  any  distribution  of  charged  particles  in  phase  space  following 
Sacherer’s  theory.  In  particular,  the  distribution  that  is  consistent  with  the  linear  self  and 
focusing  forces  is  a  hyperelhpsoid  in  the  four  dimensional  phase  space  that  gives  a 
uniform  density  profile  in  the  transverse  plane  (jc,  y)  .  The  hyperelhpsoid  can  be 
described  by  a  symmetric  matrix  M  through  the  equation 

//>(*>0>/>s)  =  /*(k7T)  M'  (s)-M  '(O-x(s)-l).  (33) 

where  the  superscript  ”  denotes  the  transpose  operation  of  a  matrix  and  £(v)  is  the 
Direc-Delta  function.  At  any  position  s  ,  the  distribution  has  the  elliptical  boundary 
Xr(4M“'(5>x(.5)=l  .  Making  use  of  Eqs.  (3 1 )  and  (33),  we  have 

M{s)=T{s,s0)-M{s0)Tt(s,s0).  (34) 

The  evolution  of  the  beam  distribution  is  fully  characterized  by  the  distribution  matrix 
M(j),  which  evolves  according  to 

=  p(g)-M(f)+M(j)-  Fr(.v),  (35) 

ds 

where  we  have  used  Eqs.  (32)  and  (34). 

We  relate  the  ten  independent  elements  of  the  matrix  M(.9)  to  the  physical  parameters 
in  Appendix  B,  yielding 

A/,  ,(5)  =  a2  (5),  (36. a) 
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^22 (»)  -  - 6S"™^  +  a2 {s)t£ (s) , 
a-(s) 

(36. b) 

«»(*)- »2M. 

(36. e) 

b  (5) 

(36. d) 

M,2(.s)  =  tf2(sK(.s), 

(36. e) 

M13(^)  =  0, 

(36.0 

M 14 (.S')  O  (s'jcXy . 

(36. g) 

A/23(.v)=  -/)2(^)av(.v). 

(36.h) 

M24  (5)  =  a 2  {s)mx  C0«v  (s)~b2  (s)t'y  (-0, 

(36. i) 

M34{s)  =  b2(s)py{s), 

(36.j) 

where  and  b(s)  are  the  semi  axes  of  the  ellipse  whieh  bounds  the  beam  distribution 
in  the  configuration  spaee,  £xrms  (0  and  £yrms (0  are  the  mis  emittanees  in  the  x—  and 
y-  directions,  respectively.  Furthermore,  /jx (.9),  av(0’  and  ^v(.s’)  are  the  flow 

velocities  as  defined  in  Eq.  (6)  with  6  =  0  . 

We  use  the  evolution  equation  of  the  matrix  M(s),  i.e,  Eq.  (35),  to  determine  the 

evolution  equations  for  £(^),  A* CO’  /*y(0>  <2*  (O’  a\(s)^  sxrmsis)  anc*  €vrms (0 • 

Substituting  Eq.  (36)  into  Eq.  (35),  the  first  and  third  diagonal  elements  give  the 
equations  for  the  flow  velocities, 


t*xis)  = 


1  da(s ) 
a(s)  ds 


(37) 
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(38) 


My{s)  = 


1  db(s ) 
^(5)  ds 


The  third  element  on  the  first  row  of  Eq.  (35)  gives  a  relation  for  the  rotational  flow 
velocity,  i.e.. 


«vfe)  =  £lt) 

«y(5) 


(39) 


We  use  the  second  element  on  the  first  row  and  the  fourth  element  on  the  third  row  of  Eq. 
(35)  to  determine  the  envelope  equations  for  the  beam  semi-axes,  yielding 


“T"  +  K-q{s)a{s)-  2^K-:(s)av(s]a(s)- 
ds~ 


2K 


1 6ei 


a(s)+b(s)  o3(j) 


-  Kq  {s)b{s)~  2y[tcjs)ax  {s)b{s) — 
ds*  at?  +1 


2 

l  6Syrms 


<?(<>)+ />(s)  b^(s) 


(40) 


(41) 


where  use  has  been  made  of  Eqs.  (37)  and  (38).  Further,  the  second  element  of  the  third 
row  and  the  fourth  element  on  the  first  row  of  Eq.  (35)  determine  the  evolution  of  the 
rotational  flow  velocities  tfv(s)  and  ay(s) 


dax(s ) 


ds 


+ 


f4g) . 

_a(s)  b(s)_ 


b[s )  k(j  ds 


0, 


(42) 


da  ..(j) 


ds 


-4- 


a(s) 


k0 


ds 


(43) 


Finally,  we  use  the  second  and  fourth  diagonal  elements  of  Eq.  (35)  to  determine  the 
evolution  of  the  mis  emittances.  After  considerable  simplification  with  aid  of  Eqs.  (39)- 
(43),  we  find  that  the  evolution  equations  of  the  mis  emittances  are  two  total  differential 
equations  which  is  expressed  as,  after  integration. 
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^/-mi(-s)=T7°4('yK('yH(5)+cl  > 


(44) 


£yrms  CO  (-V  )('^v  (.V  V^y  (.V ) +  C2 , 


(45) 


where  C\  and  c2  are  integration  constants.  The  constants  c\  and  c2  are  related  to  the  rms 
thermal  emittances  which  are  defined  as 


£?h,x  =  (x2)((x'  ~vx! PbcY)  = 


a2{s) 
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£7  (5) 


e,  ,  (46) 


£ky  s  (>’2)((>’'  -  Vy/fcf  )  =  ^ 


1  6  £  yrmS  (-*'  ) 

^T 


-a2(s)tz2(s) 


—  C~t 


(47) 


where  Vx  -  (s)x -  ax (,v)y]  and  Vv  =  ft bc\juy(s)y  + a  y(s)x\  are  the  tlow  velocities 

in  the  x-  and  y  —  directions,  respectively,  and  the  average  denotes 
{/)  =  J  Xf 'bdxdydx'dy' .  It  should  be  stressed  that  the  thermal  emittances  are  kept  constant, 
although  the  rms  emittances  are  not  constant  and  evolve  as  the  beam  propagates.  In  the 
cold-fluid  limit,  ethx  =  sthy  =  0  ,  s\rms  =  aA(s)ax(s)aY(s)l\b  , 

£vnm  =  ft4(j)aJf(5)a>;(5)/l6 ,  and  Eqs.  (40)-(43)  recover  the  cold-fluid  results  in  Eqs. 
( 1 5)-(  1 8)  for  0  =  0. 
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IV.  Examples  of  intense  elliptic  beam  equilibria 

In  this  section  we  present  two  examples  of  periodically  focused  large-aspect-ratio 
elliptic  beam  equilibria  in  a  hybrid  magnetic  field  configuration  consisting  of  a  periodic 
non-axisymmetric  magnetic  focusing  field  and  a  quadrupole  magnetic  field.  One  example 
is  a  high-space-charge  elliptic  beam  for  a  high-efficiency  200  W  ribbon-beam  amplifier 
(RBA)  under  development  at  Massachusetts  Institute  of  Technology  (MIT)  and  Beam 
Power  Technology  for  wireless  communication.  The  other  is  a  relativistic  elliptic  beam 
which  can  be  used  in  a  10  MW  L-Band  ribbon-beam  klystron  (RBK)  for  the  International 
Linear  Collider  (ILC). 

To  aid  the  design  of  high-space-charge  elliptic  beam  transport,  a  numerical  module  in 
the  PFB2D  code  [5,  12]  has  been  developed  to  solve  the  generalized  envelope  equations 
(40)-(43),  which  determines  the  required  magnetic  field  for  the  desired  major-axis  a(v) 
and  minor-axis  b(s)  of  the  beam  ellipse.  Those  information  are  used  in  2D  self-consistcnt 
PIC  PFB2D  simulations  and  3D  OMNITRAK  self-consistent  particle  trajectory 
simulations  to  verify  the  theory. 

As  the  first  example,  we  consider  a  space-charge-dominatcd  6:1  elliptic  electron  beam 
with  desired  beam  semi-axes  a  =  0.373  cm  and  b  =  0.062  cm  propagating  with  current 
Ih  =  0. 1 1 A  along  a  beam  tunnel  with  a  constant  axial  potential  of  =  2290  V.  For  this 

beam,  we  have  ~  0.094  and  yh  =  1 .0045 .  Based  on  the  discussion  in  Ref.  [  1 9],  we  can 
take  the  intrinsic  thermal  emittance  to  be  zero  in  this  case,  i.e.,  £th  -  £thy  =  0.  We  use 

the  PFB2D  code  to  solve  the  set  of  the  generalized  envelope  equations  ( 1 5)-(  1 8)  [or  Eqs. 
(40)-(43)]  to  determine  the  required  magnetic  fields  for  focusing  such  a  beam.  The 
parameters  of  the  hybrid  periodic  permanent  magnet  focusing  fields  are  found  to  be 
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5r(s)  =  -263sin(Ar0s)  G  ,  £^=30.0G/cm,  5  =  1.912  cm,  k0y/k0x  =  6  (see  Table  1, 

Column  2).  As  seen  in  Fig.  2,  the  matched  beam  semi-axes  and  A(s)  (solid  curves) 
solve  the  generalized  envelope  equations  (15)-(18).  The  semi-major  axis  is  almost  a 
constant,  whereas  the  semi-minor  axis  oscillates  slightly  about  a  constant  value. 

Shown  in  Fig.  2,  the  dotted  curves  are  the  envelopes  of  the  beam  ellipse  obtained  from 
the  PFB2D  self-consistent  PIC  simulation  for  the  example  of  a  nonrelativistic  elliptic 
beam  with  voltage  1^=2290  eV  ,  current  Ih  =  0.11  A,  aspect  ratio  a/b-  6  and 

sth  =  £th  =  0  in  the  hybrid  magnetic  field  with  Bz(s)=  -263sin(A'0.v)  G  , 

B'q  =  30.0G/cm,  5  =  1.912  cm,  and  k0}/k0x  =6.  In  the  simulation,  5x10  particles  are 

used.  The  thermal  emittance  is  negligibly  small,  and  the  cold-fluid  approximation  is  valid. 
As  shown  in  Fig.  2,  there  is  excellent  agreement  between  the  theoretical  envelope 
solution  (solid  curves)  and  the  self-consistent  PIC  simulation  results  (dotted  curves).  We 
measured  the  angle  of  the  beam  ellipse  with  respect  to  the  laboratory  frame.  The  angle 
oscillations  have  a  small  amplitude  of  0.1  degrees,  which  are  due  to  fluctuations  (noise) 
in  the  simulation. 

The  PFB2D  simulation  also  shows  that  the  transverse  beam  distribution  preserves  the 
equilibrium  profile  as  it  propagates.  In  Fig.  3,  5,000  particles  (a  sample  of  the 
5  x  105  particles  in  the  PFB2D  simulation)  are  plotted  in  the  (*,>>)  plane  and  (\\dy/ds) 
plane  for  five  snapshots  within  one  period:  s/S  =  9.0, 9.25, 9.5, 9.75  and  10.0  for  the  same 
elliptic  beam  shown  in  Fig.  2.  The  results  in  Fig.  3  also  suggest  that  the  beam  equilibrium 
is  stable. 
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As  a  separate  verification  of  the  theories  and  PFB2D  code,  a  3D  OMNITRAK  [13,  19] 
simulation  is  performed  for  the  6:1  nonrelativistic  elliptic  beam.  Since  3D  trajectory 
simulations  are  time-consuming,  only  a  2-period  interval  is  used  for  this  test,  as  shown  in 
Fig.  4.  The  beam  is  sent  through  a  conducting  cylindrical  beam  tunnel  (not  shown)  of 
radius  10.0  mm.  In  this  simulation,  10,000  macroparticle  arrays  are  used.  However,  only 
a  fraction  of  them  (i.e.,  64  macropaticle  arrays)  are  shown  in  Fig.  4(a).  Substantially 
parallel,  non-twisting  transport  is  achieved,  and  the  simulated  beam  envelopes  agree  with 
the  theoretical  predictions  as  shown  in  Fig.  4(b). 

As  the  second  example,  we  consider  a  relativistic  elliptic  beam  that  can  be  used  in  a  10 
MW  L-Band  ribbon-beam  klystron  (RBK)  for  the  International  Linear  Collider  (1LC). 
The  beam  has  a  current  of  Ib  -  1 1 1.1  A,  a  voltage  of  Vb  =  120kV  and  an  aspect  ratio  of 

20:1,  which  corresponds  to  fib  =  0.094  and  yb  =  1 .0045 .  The  other  beam  parameters  are 
listed  in  the  third  column  of  Table  1.  Solving  the  generalized  envelope  equations  (40)- 
(43)  with  £tbx  -  £{hy  -  0,  the  hybrid  magnetic  fields  are  determined  to  be  the  form  of 

Eq.  (4)  with  5z(s)  =  -2000sin(A'0x)  G  ,  =  80.8  G/cm ,  S  =  2.2  cm ,  and  A'0  v jk0x  =  20 . 

In  Fig.  5,  the  solid  curves  are  the  beam  semi-axes  tf(s)  and  />(s)  calculated  from  the 
generalized  envelope  equations  (40)-(43)  with  £thx  =  £tb  v  =0,  whereas  dotted  curves 

are  from  the  self-consistent  PIC  PFB2D  simulation. 

For  practical  devices,  elliptic  beams  usually  have  nonzero  temperature  due  to  the 
heating  of  sources  or  other  effects.  To  study  the  temperature  effects,  we  solve  the 
generalized  envelope  equations  (40)-(43)  with  nonzero  initial  thermal  emittances,  i.e., 

£jjl  x  =  kBTa2(s  =  0 )/niyb  flbc 2  and  £2h  y  =  kBTb2{s  =  0 )/myb  pic2  .  As  shown  in  Fig.  6, 
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the  elliptic  beam  envelopes  are  calculated  for  three  different  temperature  choices:  0  eV, 
50  eV  and  100  eV.  Compared  with  the  cold  beam  envelopes  shown  as  solid  curves  in  Fig. 
6,  the  warm  beam  envelopes  are  found  to  increase  slightly  as  the  beam  temperature 
increases,  while  the  aspect  ratio  of  the  beam  decreases  from  20:1  to  16.4:1  as  the  beam 
temperature  increases  from  0  eV  to  100  eV. 

For  the  second  example  of  elliptic  beams,  we  assume  that  the  elliptic  beam  is 
generated  from  an  electron  gun  with  an  intrinsic  temperature  of  0.1  eV  and  a  current 
density  of  1.5  A/cm  [19].  The  elliptic  beam  has  to  be  compressed  by  a  factor  of  471.5  in 
area  to  achieve  a  current  density  of  707.3  A/cm*'  in  the  focusing  channel.  During  the 
compression,  the  temperature  increases  by  a  factor  of  47 1 .5  to  47.2  eV.  Therefore,  in  our 
calculations,  a  temperature  of  50  eV  is  a  reasonable  assumption.  As  shown  in  Fig.  7,  the 
envelopes  of  the  elliptic  beam  with  a  temperature  of  50  eV  are  obtained  by  solving  the 
envelope  equations  (40)-(43)  (solid  curves)  and  by  the  PFB2D  PIC  simulations  (dotted 
curves).  Both  results  showed  a  slight  increase  in  the  two  envelope  dimensions  and  the 
aspect  ratio  of  the  elliptic  beam  decreases  to  17.8,  compared  with  the  cold  elliptic  beam. 
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V.  Conclusions 


A  self-consistent  solution  for  the  focusing  of  coasting,  sheet-like,  space-charge- 
dominated  elliptic  beams  using  a  hybrid  of  non-axisymmetric  periodic  permanent 
magnets  (PPM)  and  quadrupole  magnets  was  developed.  The  beam  envelope  equation 
and  equilibrium  flow  profiles  were  obtained  using  a  paraxial  cold-fluid  model.  Further,  a 
kinetic  theory  was  developed  to  determine  the  self-consistent  beam  distribution  and  its 
evolution.  Both  the  cold-fluid  and  kinetic  equilibrium  theories  were  validated  by  the  two- 
dimensional  (2D1  particle-in-cell  (PIC)  simulations  using  the  2D  Periodically  Focused 
Beam  (PFB2D)  code,  and  the  three-dimensional  (3D)  self-consistent  trajectory 
simulations  using  OMNITRAK.  The  theories  and  PfB2D  and  OMNITRAK  simulations 
were  applied  to  design  such  high-space-charge  elliptic  beams  for  applications  in  klystrons 
and  vacuum  electron  devices.  Numerical  results  showed  that  the  beam  edges  in  both 
transverse  directions  are  well  confined  without  twisting  and  the  beam  density  profile  is 
well  preserved.  For  space-charge-dominated  elliptic  beams  used  in  ribbon-beam  klystron 
applications,  the  temperature  effects  are  studied  using  both  the  kinetic  equilibrium  theory 
and  PFB2D  simulations. 
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Figure  Caption 


Figure  1  Twisted  coordinates. 

Figure  2  Plots  of  the  beam  envelopes  a(s)  and  b(s )  versus  the  axial  distance  s  for  the 
nonrelativistic  elliptic  beam  in  Table  1.  The  solid  curves  are  the  generalized 
envelope  solution,  whereas  the  dotted  curves  are  from  the  PFB2D  simulation. 

Figure  3  Plots  of  5,000  particles  (a  sample  of  the  1  x  105  particles  in  the  PFB2D 
simulation)  in  the  (jc,y)  plane  and  ( x^dy/ds )  plane  for  five  snapshots  within 
one  period:  s/S  =  9.0, 9.25,9.5, 9.75  and  10.0  for  the  same  beam  shown  in  Fig. 
2. 

Figure  4  3D  OMN1TRAK  simulation  results:  (a)  image  of  the  simulated  beam  and  (b) 
plots  of  the  beam  envelopes  from  the  3D  OMN1TRAK  simulation  (dotted 
curves)  and  the  theory  (solid  curves).  Here,  the  parameters  of  the  nonrelativistic 
ellipse-shaped  beam  are  listed  in  Table  1. 

Figure  5  Plots  of  the  beam  envelopes  a(s)  and  />(.v)  versus  the  axial  distance  s  for  the 
relativistic  elliptic  beam  in  Table  1  with  zero  temperature.  The  solid  curves  are 
the  generalized  envelope  solution,  whereas  the  dotted  curves  are  from  the 
PFB2D  simulation. 

Figure  6  Plots  of  the  beam  envelopes  (a)  a(s)  and  (b)  b(s)  versus  the  axial  distance  s  for 
the  relativistic  elliptic  beam  in  Table  1  for  three  different  temperature  choices: 
OeV,  50  eV  and  100  eV. 

Figure  7  Plots  of  the  beam  envelopes  a(.s')  and  b(s)  versus  the  axial  distance  s  for  the 
relativistic  elliptic  beam  in  Table  1  for  a  temeprature  of  50  eV.  The  solid  curves 
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are  the  generalized  envelope  solution,  whereas  the  dotted  curves  are  from  the 
PFB2D  simulation. 
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Table  1  System  parameters  for  elliptic-beam  examples 


Parameter 

Nonrelativistic 

Relativistic 

Wireless 

ILC 

Application 

Communication 

Frequency  (GHz) 

1.95 

1.3 

RF  Power  (kW) 

0.2  (cw) 

10  (pulsed) 

Current  (A) 

0.11 

111.1 

Voltage  (kV) 

2.29 

120 

a/b 

6.0 

20 

a  (cm) 

0.373 

1.0 

S  (cm) 

1.912 

2.2 

^0  y  Ojc 

6.0 

20 

B0  (kG) 

0.263 

2.0 

B'q  (G/cm) 

30.0 

80.8 
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28 


Envelopes  (cm) 


0.5 


.  PIC 

a(s) 

b(s) 
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5  50  -0.25  0  00  0  25  0  50 

x  (cm) 


3  50  -0  25  0.00  0.25  0.50 

x  (cm) 


&DO0 


3.50  -0.25  0  00  0.25  0.50 

x  (cm) 


Fig.  3 


30 


Fig.  4 


31 


Envelopes  (cm) 


Fig.  5 


32 


1.10 


1.05 


E 

2 

Q 


1  00 


0.95 

0 


(a) 


M ••••••• ••#«•••• •! 


iMMlItllMl 


5.5 


11.0 

s  (cm) 


16.5 


22.0 


s  (cm) 


Fig.  6 


33 


1.5 


Theory 

PIC 


E  1 .0 

CO 

<D 

Q_ 

O 

<U 

> 

c  0.5 
LU 


0.0 


a(s) 


b(s) 


o.o 


5.5 


11.0 

s  (cm) 


16.5 


22.0 


Fig.  7 


34 


Appendix  A 


The  derivation  of  the  generalized  envelope  equations  (8)-(  1 2)  is  as  follows.  Wc  first 
solve  the  continuity  equation  (1).  Substituting  the  density  profile  in  Eq.  (5)  and  the  flow 
velocity  profile  in  Eq.  (6)  into  Eq  (1),  we  obtain 

a'  b' 


f.ix  +  //, 

a  b 


© 


~2  ~2 

£ _ y_ 

2  ,2 

a  b 


+  2 


a 

- Mx 

Va  J 


■  + 


b' 

~~Vy 


\  < 


+  1  -  —  0'  +  —  O'  +—ax-—av 
1  a  b  a  4  b  y 


\  ~~  — 

r  ~2  ~2i 

|  xy 

s 

j  ab 

£J2  b2\ 

(A.l) 


=  0, 


where  the  prune  denotes  8/ ds ,  and  the  relations  x'-y6',  y'--x6',  dc  x  /ds  =  e  ^6' 
and  dex  jds  =  -evt9'  have  been  used.  Since  Eq.  (A.l)  must  be  satisfied  for  all  x  and  y , 

the  coefficients  of  the  terms  proportional  to  ©  ,  x2S  ,  y2S  and  xyS  must  vanish 
independently.  This  leads  to  the  following  equations 

da 


l‘x  ,  » 

a  ds 


1  db 
b  ds 


Vy=T  —  > 


d6  a2ay-b2a 


(A.l) 

(A. 3) 

(A.4) 


ds  a2  -b1 

Second,  we  solve  the  force  balance  equation  (3).  Substituting  the  density  profile  in  Eq. 
(5),  the  flow  velocity  profile  in  Eq.  (6),  and  the  self-field  potential  in  Eq.  (7)  into  Eq.  (3), 
we  obtain  two  equations  in  the  x  -  and  y  -  directions,  respectively.  In  the  x  -  direction, 
we  have 
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pic1  [(//;  -  ay&  +  -  axay  +  ay0'pc  +  (-  a'x  +  /.ix&  -  nyff  -  ax/ix  -  ax/iy  ).y  ] 

^7-— -  PbBy  +  (pvy  +  ay*)PbBz 


7  bm 


7b 


2  a  +  b  a 


In  the  y  —  direction,  we  have 

PhC2  -  /.i yff  +  f-ix6 '  +  ayf.ix  +  (XyfJy )x  +  (//',  +  ay6' -  ax0’  +  f.i2y-  ax ay  )y 


J_ IsHl  "-1  +  pbB~  -  (pxx  -  axy)phB. 
a  +  b  a 


Ybm 


where 


2  K. 

fx  =  Mx  -  ayQ'  +  Mx  -  axay  +axB' - 7~ — T7  _  2^/a:-  (.v)«v  h - - — —  B ^ 


a{a  +  b) 


7bPbmc 


2  2  K 


/,  =  //'v  +  ayff -  ax9'  +  /p  -  axav  -  ■  -  2^/c:(s)ax - C/ — T B ^ , 

b{a  +  b)  7bPbmc 


gx  =  a'y  -  t-iyd' +  /.ix0'  +  ayf!x  +  avfiv  +  2xJkz(s)/jx - - — =- 5^ , 

7bPbmc 

g  y  —  OCx  +  fiyO  f-lyd  C(  xf.ix  CCy^ly  (.V  )/^  H  ? 


7bPhmc 


with  the  definitions  of 


JZfr)  e  AifL  and  ft  s 


2yhphmcL  ylPW 

We  can  express  the  magnetic  field  in  Eq.  (4)  in  the  twisted  coordinate  as 


(A. 5) 


(A. 6) 


We  can  rewrite  Eqs.  (A. 5)  and  (A. 6)  as 

f,x+g,y  =  o. 

(A. 7) 

g*x+fyy  =  0. 

(A. 8) 

(A. 9) 


(A.  10) 


(A.ll) 


(A. 12) 


(A- 13) 
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(A. 14) 


B  Bz (.s' )e .  +  {bxxx  +  Bxyy  +  [Byxx  +  B y^,  y'jc~ , 

where  the  magnetic  fields  contributions  are 


dB;{s) 


.v.v 

ds 

%  = 

dB.(s) 

ds 

%  = 

dB,  (s) 
ds 

Byy  = 

dB.(s) 

ds 

&  —  *" 

,2 


cos2  #(.?)+ -^-sin2  &(s) 


+  2 B'  sin#(.s)cos0(,y), 


(A. 14) 


,2  a  \ 
*qx_  +  _2l 
'  a  a 

y  K0  *1)  ; 


,2  /  2  \ 
k0x  .  k°y 

L-2  L2 

y  A0  A0  j 


cos#(s)sin<9(s)+.6^  cos2  4s)  -sin  2  4*)].  (A- 15) 


cos#(.?)sin<9(.9)  +  B'q  cos2  4s)  -sin  2  4-4  (A- 16) 


-^-sin  #(s)  + — —  cos2  6>(s) 

k1  k2 

Ao  h0 


-2 B'  sinks’)  cos  4A 


(A.  17) 


Since  Eq.  (A. 7)  and  (A. 8)  must  be  satisfied  for  all  3c  and  y ,  the  coefficients  of  the  terms 
proportional  to  3c  and  y  must  vanish  independently,  which  lead  to 


II 

o 

(A.  18) 

o 

II 

(A.  19) 

g,  =  0  , 

(A.20) 

gy  =  0  • 

(A. 21) 

Substituting  Eqs.  (A.2)-(A.4)  into  Eqs.  (A.18)-(A.21),  we  obtain 
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fx  =  Bx  -  ax°' +  Bx  ~  axay  +  ax°'  - 


2  -  ~  '  ~  n'  2K  -2xjfc:(s)a>.  +  7 


a(a  +  b ) 


YbPbmcl 


d  l  a'  \  ( a  ^ 


ds 


a 

J 


+ 


\a  J 


2  2 

/  \a^av  —  b  a 

\ax-ay) - - 


2  K 

.  - Cl  Cl - 

y)  2  j 2  a  jv 


az-/V 


a(a  +  /)) 


-2^K.(s)i 


\av  + 


qh 


YbPbmc 


1  d2a  b2(a2-2axay)+a 


a  ds1 


2«v  2 K  ^  ,  qB 

,  i  — - —  -2JKAs)a„  + - ! 

a2-/)2  a(a  +  6)  V  “w  ' 


7  r*A»«2 


=  o. 


(A. 22) 


2.K 

fv  =  //'  +  «v^'  -  ar<9'  +  f.i2  -  axav - 2 yjicjs)ax  - 


+  b) 


YbPbmc 


d_ 

ds 


fh'\  sv\2 


+ 


b'Y  (  \a  av-b  ax 


V  v  ) 


-  {ax  ~  ay  ) 


2  K 


--aYa„  -- 


b)  VA  ""  a2-/)2  v  y  6(a  +  6) 


-2jicjsjax- 


qB~v 


YbPb,nc‘ 


d2b  a2(a)-2axay)+b2i 


ax2  2  K 


bds 2 
=  0, 


2 

a  -  b 


b(a  +  b) 


-2  fijsjax- 


qB 


XV 


YbPbmc 


(A. 23) 


qBx 


Sx  =  a'v  -  By0'  +  /A*9'  +  CtyMx  +  avBv  +  2x[^)Bx  “  ‘  2 

YbPbmc 


<  +(Bx-  By) 


2  2 

a  av-b  aY 

t~3-  +  ay\-+b 


az-y 


a[  b'^ 
a 


2^kz{s)Bx 


CJBXX 


YbPbmc~ 


(A. 24) 


1  d  (  2  )  \bi(ax-ct)d 

~  - U*  av  I - o - : - 

a2  ds^  a  a^-b2  ds{b j 


+  2—t[kJs)- 


qBx 


a  YbPbmc 


=  0. 


qBy 


gy  =  -a'x  +  Bx°'  -  By°'  -  axBx  -  axBy  -  2  JkJs)bv  +■  2 

■  YbPbmCl 


-  ~a'x  +  (/fx  ~Bv) 


a2  a  x  —b2ax 

2  72 

a  —  0 


a' 


-«v  — +  — 
a  b 


-2 -JkJs)b 


qB~y 


YbPbmc2 


(A. 25) 


=  — [T-(b2ax)+-^  -)-2^)-  +  ■ 

6  a2-/)2  ‘  7  ft  YbPbmcl 


=  0. 
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The  generalized  envelope  equations  are 


2  /  2  \  2  2 

d~a  b  \ax  -2axav  J+a  av  2 K  .  / — 7-r 

— r - 1 ^ - —  a - 2  Jk.  (5  )a..a  + 

ds2  a2~b 2  (a +  6)  V  *  ' 


YbPbmcl 


-a  =  0 , 


t/2/l '  t2.  v )+  ^  h  _  _2K_  _  _ 


ds 


2  1 2 

a  -b 


(a  +  b) 


d^xy 

nPbmc2 


-b  =  0, 


£\ 

[rayi 

UJ 

1 

a 

t: 

)  d 

ds 

2  1 2 

(7  —  b 

ds 

UJ 

d  | 

(.  2  ) 

a2b[ax-av) 

d  | 

ds 

1 

Ci 

tO 

1 

O 

to 

fife  1 

>7; 

+  la'a^l fc:(s) - -  a2  =  0  , 


YbPbmc 


Q^yy  ,2  _ 


+  2b'bMcJs)~ 


r  =0, 


d6  a2a  -b2a v 


<&  a2  -  62 


(A. 26) 

(A.27) 

(A. 28) 

(A.29) 

(A. 30) 


Substituting  Eqs.  (A.14)-(A.17)  into  Eqs.  (A.26)-(A.30),  we  obtain  the  set  of  generalized 
envelope  equations  (8)-(  1 2). 
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Appendix  B 

In  this  Appendix,  we  determine  the  elements  of  the  matrix  M(s).  In  order  to  determine 
the  elements  of  the  matrix  M(s),  we  make  the  projection  of  a  four-dimensional  (4D) 
hyperellipsoid  to  a  two  dimensional  ellipse.  Let  us  first  define  the  4D  hyperellipsoid 
through  the  equation 

Xrfc)-L(s)-X(s)=l>  (B.l) 

where  L(s)  =  M  1  (5)  and 

(r  \ 

x\ 

x2 

1  = 

*3 

We  expand  Eq.  (B.  1 )  as 

1  =  L\  ,xf  +  L21x  2  +  E33xf  +  L44X4 

+  2L|2A|X2  +  2L|3X|X3  +  2L|4.V|A'4  +  2L23x2x3  +  2Z,24x2x4  +  2L34x3x4. 

To  project  the  4D  hyperellipsoid  defined  by  Eq.  (B.3)  onto  the  subspace  (.Y|,x2,x3),  we 
find  the  extremal  points  in  these  coordinates  by  differentiating  Eq.  (B.3)  implicitly  and 
setting  clx |  =  d.x2  =  t/x3  =  0 .  In  this  case,  we  find 


(B.2) 


(B.3) 


x,  =  -- 


-(Z.|4X]  +-£-24x2  +L34x3), 


(B.4) 


-44 


which  we  substitute  into  Eq.  (B.3)  to  obtain  the  equation  for  the  3D  ellipsoid  that  bounds 
the  projection  of  the  4D  hyperellipsoid  defined  by  Eq.  (B.3)  onto  the  subspace  (x|,x2,x3). 


i.e.. 


^44  -  (A  1^44  “  L\4  k2  +  (^22^44  ~  E24f  2  +  (T33/.44  -  Z.34  )t3 

+  2(A|  1 T44  —  LuL24  )x ( x 2  +  2(  L\  3Z,44  -L|4Z,34)x1x3  +  2(L23L44  -  L24L34)x2x3. 


(B.4) 
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The  3D  ellipsoid  of  Eq.  (B.4)  can  be  further  projected  onto  the  2D  subspace  (x| ,  x2 ) . 
We  again  differentiate  Eq.  (B.4)  implicitly  and  set  dx \  =  dx2  =  0  to  obtain 


x  _  (^14-^34  ~  ^13^44  )*1  tiL 24-^34  ~  ^23^44  )x2 

L33L44  -  L 34 


(B.5) 


which  we  substitute  into  Eq.  (B.4)  to  obtain 


_  ^44^33  ~  2Z-13Z.14Z.34  +  Lx  |  Z-34  +  Tj3L44  1, 1T33T44  2 


1  = 


+  2 


^34  ^33^44 


X, 


L|4L24^33  ~  ^14-^23^34  ~  ^13^24^34  +  ^12^34  +  ^13^23^44  ~  ^12^33^ 


44 


^34  “  ^33^44 


XjX2  (B.6) 


2 

,  ^24^33  -  2^23^24^34  +  ^22^34  +  ^23^44  _  ^22^33^44  ..2 

+ - “2 - x2- 

^34  _  ^33^44 


Equation  (B.6)  is  the  2D  ellipse  boundary  of  the  projection  of  the  4D  hyperellipsoid 
defined  by  Eq.  (B.l)  onto  the  subspace  (xi,x2). 


Expressing  the  elements  in  terms  of  the  elements  A/,;  via  the  relation 


L(s)  =  M_l  (5),  we  obtain  after  lengthy  manipulation, 


1  = 


M 


22 


A/ j  j  A/22  —  A/  1 2 


-xj-2 


M 


12 


M  |  j  h/1 22  ilf  ]  2 


rx,x2  + 


Mi 


M\\M22  ~  Af|2 


2  x22' 


(B.6) 


Figure  B.l  The  boundary  ellipse  defined  in  Eq.  (B6). 


41 


Note  that  there  are  a  few  useful  properties  associated  with  the  ellipse  defined  in  Eq.  (B.6). 


The  area  bounded  by  the  ellipse  is  n^MuM22  -  Afj^  .  The  maximum  extents  of  the 

distribution  in  the  X|  -  and  x2  -  directions  are  yjMu  and  -y]M22  ,  respectively,  as 
shown  in  Fig.  B.l. 

Since  the  projection  of  the  4D  hyperelhpsoid  distribution  onto  any  2D  subspace 
produces  a  uniform  distribution  bounded  by  the  ellipse  defined  in  Eq.  (B.6),  we  can  relate 
the  M (s)  matrix  elements  to  the  statistical  averages, 
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It  follows  from  Eqs.  (B.7)-(B.9)  that, 


xl  )\x2  /  (xi*2 )  0^1 1 M22  -  M]2). 


(B.IO) 


We  use  the  result  in  Eq.  (B.6)  to  project  the  4D  hyperellipsoid  defined  in  Eq.  (33)  onto 
the  subspace  (x,_y)  =  (x,  ,x3) .  By  the  definition  of  the  density  profile,  such  a  project 


should  produce  an  elliptic  boundary  defined  by 


-  + 


y 


aHs)  b2(s) 


(Bll) 


Comparing  the  coefficients  in  Eqs.  (B.6)  and  (B.7),  we  find  that  Mu(.9)=  a2(.v)  , 

A/]3 (s)  =  0,  and  (s)  =  62(s),  which  give  Eqs.  (36. a),  (36. f),  and  (36. c),  repsectvicly. 

We  project  the  4D  hyperellipsoid  defined  in  Eq.  (33)  onto  the  subspace 
(x,x')  =  (jcj , JC2 )  to  yield  an  elliptic  boundary.  We  find  that  the  usual  nns  emittance  is 
related  to  the  area  of  the  ellipse  as 

£xrms  CO  =  Pj(^f  =  ^MUM22  ~  ,  (B.12) 
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where  £xrms.(s)  is  the  nns  emittance  in  the  x—  direction.  Using  the  relation  in  Eq.  (B.9), 


the  correlation  coefficient  of  *  and  x'  is  (xx'Y/7x2\  =  ]  ,  which  allows  us  to 


relate  the  flow  velocity  defined  in  Eq.  (6)  to  the  element  A/)2  as  [see  Eq.  (36. e)] 


Mi2(s)  =  a2(sK(s), 


(B.13) 


where  use  has  been  made  of  (jcx')  =  (xVxj =  (vxX1')  -  •  Solving  Eq.  (B.12) 


yields 


Ma{s)- -6-if &  +  a2(s)tf(s), 
a  (. s  ) 


(B.  14) 


which  gives  Eq.  (36. b). 

In  a  similar  manner,  we  determine  the  corresponding  elements  of  the  matrix  M(s)  for 
the  subspace  (y,y ')=  (x3,x4),  i.e., 

A^34(S)=  (B.15) 


and 


,  v  16f2  (s)  2/  \  2 1  \ 

M44 (s)  =  ,27\  +  b  (5W  (s) ,  (B.  1 6) 

b  (5) 

which  give  Eqs.  (36 .j)  and  (36. d),  respectively.  Here,  //v(s)  is  the  flow  velocity  in  the 
v  -  direction,  and 

Gyms (5)  =  )(a'2 )  "  O’/)2  =  ^M33M44-Afl  (B.17) 

is  the  nns  emittance  in  the  y-  direction. 
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Wc  project  the  4D  hyperellipsoid  defined  in  Eq.  (33)  onto  the  subspacc 
(x,  v')=  (jc, , jc4 )  to  yield  an  elliptic  boundary  describing  the  coupling  of  the  particle 
motion  between  the  x-  and  y  -  directions.  The  projection  involves  the  elements 
A/n(.9)  ,  M,4(.?)  and  A/44(,?)  .  We  use  the  relation  in  Eq.  (B.9),  i.e., 

(.vv')/ lx2)  =  M  )4  M  n  and  the  definition  of  the  flow  velocity  in  Eq.  (6)  to  express 


Ml4(s)=  a2  (.?)«>,(.?)  (B.  IS) 

where  or  (s)  is  the  rotational  flow  velocity  in  the  y—  direction.  Similarly,  we  have 

M23(5)=-62^)aJ((5),  (B.  19) 

where  av(s)  is  the  rotational  velocity  in  x-  direction.  Equations  (B.18)  and  (B.19)  give 
Eqs.  (36. g)  and  (36.h),  respectively. 

Finally,  we  project  the  4D  hyperellipsoid  defined  in  Eq.  (33)  onto  the  subspace 
(*',_/)  =  (a'2,a4)  .  The  correlation  of  the  particle  motion  between  the  jc'-  and 
v'- directions  is  calculated  to  be 


M24{s)=(xy)  ----  }  +  }  =  a 2 (,v )//r (s )ay (,v )-b2{s )//v (.9 )orv (.v ) . 
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(B.20) 


which  gives  Eq.  (36. i). 
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NON-AXISYMMETRIC 
CHARGED- PARTICLE  BEAM  SYSTEM 

PRIORITY  INFORMATION 

71ns  application  claims  priority  from  provisional  appli¬ 
cation  Ser,  No.  60/577 ,1 32  filed  Jun.  4,  2004  which  is 
incorporated  herein  by  reference  in  its  entirety. 

BACKGROUND  OF  THE  INVENTION 

'Hie  invention  relates  to  the  Held  of  charged-particle 
systems,  and  in  particular  to  a  non-axisymmetric  charged- 
particle  system. 

’Hie  generation,  acceleration  and  transport  of  a  high¬ 
brightness,  space-charge -dominated,  charged-particle  (elec¬ 
tron  or  ion)  beam  are  the  most  challenging  aspects  in  the 
design  and  operation  of  vacuum  electron  devices  and  par¬ 
ticle  accelerators.  A  beam  is  said  to  be  space-charge-domi- 
nated  if  its  self-electric  and  self-magnetic  field  energy  is 
greater  than  its  thermal  energy.  Because  the  beam  brightness 
is  proportional  to  the  beam  current  and  inversely  propor¬ 
tional  to  the  product  of  the  beam  cross-sectional  area  and  the 
beam  temperature,  generating  and  maintaining  a  beam  at  a 
low  temperature  is  most  critical  in  the  design  of  a  high- 
bnghtness  beam.  If  a  beam  is  designed  not  to  reside  in  an 
equilibrium  state,  a  sizable  exchange  occurs  between  the 
field  mid  mean-flow  energy  and  thermal  energy  in  the  beam. 
When  the  beam  is  space-charge-dominated,  the  energy 
exchange  results  in  mi  increase  in  the  beam  temperature  (or 
degradation  in  the  beam  brightness)  as  it  propagates. 

If  brightness  degradation  is  not  well  contained,  it  can 
cause  beam  interception  by  radio-frequency  (RF)  structures 
in  vacuum  electron  devices  and  particle  accelerators,  pre¬ 
venting  them  from  operation,  especially  from  high-duty 
operation.  It  can  also  make  the  beam  from  the  accelerator 
unusable  because  of  the  difficulty  of  focusing  the  beam  to  a 
small  spot  size,  as  often  required  in  accelerator  applications. 

Hie  design  of  high-brightness,  space-charge-dominated, 
cliarged-partiele  beams  relies  on  equilibrium  beam  theories 
and  computer  modeling.  Equilibrium  beam  theories  provide 
the  guideline  and  set  certain  design  goals,  whereas  computer 
modeling  provides  detailed  implementation  in  the  design 

While  some  equilibrium  states  are  known  to  exist,  match¬ 
ing  them  between  the  continuous  beam  generation  section 
and  the  continuous  beam  transport  section  has  been  a 
difficult  task  for  beam  designers  and  users,  because  none  of 
the  known  equilibrium  states  for  continuous  beam  genera¬ 
tion  can  be  perfectly  matched  into  any  of  the  known  equi¬ 
librium  states  for  continuous  beam  transport. 

For  example,  the  equilibrium  state  from  the  Pierce  diode 
in  round  two  dimensional  (2D)  geometry  cannot  be  matched 
into  a  periodic  quadrupole  magnetic  field  to  create  a 
Kapachinskij -Vladimirskij  (KV)  beam  equilibrium.  A  rect¬ 
angular  beam  made  by  cutting  off  the  ends  of  the  equilibrium 
state  from  the  Pierce  diode  in  infinite,  2D  slab  geometry’ 
ruins  the  equilibrium  state. 

However,  imperfection  of  beam  matching  in  the  beam 
system  design  yields  the  growth  of  beam  temperature  and 
the  degradation  of  beam  brightness  as  the  beam  propagates 
in  an  actual  device. 

SUMMARY  OF  THE  INVENTION 

According  to  one  aspect  of  the  invention,  there  is  pro¬ 
vided  a  chargcd-particlc  beam  system.  The  chargcd-particlc 
beam  system  includes  a  non-axisymmetric  diode  which 


2 

forms  a  non-axisymmetric  beam  having  an  elliptic  cross- 
section.  A  focusing  channel  utilizes  a  magnetic  field  for 
focusing  and  transporting  a  non-axisymmetne  beam. 

According  to  another  aspect  of  the  invention,  there  is 
provided  a  non-axi symmetric  diode.  The  non-axisymmetric 
diode  comprises  at  least  one  electrical  terminal  for  emitting 
charged-particles  and  at  least  one  electrical  terminal  for 
establishing  an  electric  field  and  accelerating  charged- par¬ 
ticles  to  form  a  charged-particle  beam.  These  terminals  arc 
arranged  such  that  the  charged- particle  beam  possesses  an 
elliptic  cross-section. 

According  to  another  aspect  of  the  invention,  there  is 
provided  a  method  of  forming  a  non-axisymmetric  diode 
comprising  forming  at  least  one  electrical  terminal  for 
emitting  charged-particles,  forming  at  least  one  electrical 
terminal  for  establishing  an  electric  field  and  accelerating 
charged-particles  to  form  a  charged-particle  beam,  and 
arranging  said  terminals  such  that  the  charged-particle  beam 
possesses  an  elliptic  cross-section. 

According  to  another  aspect  of  the  invention,  there  is 
provided  a  charged-particle  focusing  and  transport  channel 
wherein  a  non-axisymmetric  magnetic  field  is  used  to  focus 
and  transport  a  charged- particle  beam  of  elliptic  cross- 
section. 

According  to  another  aspect  of  the  invention,  there  is 
provided  a  method  of  designing  a  charged-particle  focusing 
and  transport  cliannel  wherein  a  non-axi  symmetric  magnetic 
field  is  used  to  focus  and  transport  a  charged- particle  beam 
of  elliptic  cross-section. 

According  to  another  aspect  of  the  invention,  there  is 
provided  a  method  of  designing  an  interlace  for  matcliing  a 
charged-particle  beam  of  elliptic-eross  section  between  a 
non-axisymmetric  diode  and  a  non-axisymmetric  magnetic 
focusing  and  transport  channel. 

According  to  another  aspect  of  the  invention,  there  is 
provided  a  method  of  forming  a  charged-particle  beam 
system.  The  method  includes  forming  a  non-axisymmetric 
diode  that  includes  a  non-axisymmetric  beam  having  an 
elliptic  cross-section.  Also,  the  method  includes  forming  a 
focusing  channel  that  utilizes  a  magnetic  field  for  focusing 
and  transporting  the  elliptic  cross-section  beam 

BRIEF  DESCRIPTION  OF  THE  DRAWINGS 

FIGS.  1A-1C  are  schematic  diagrams  demonstrating  a 
non-axisymmetric  diode; 

FIG.  2  is  a  graph  demonstrating  the  Integration  Contour 
C  for  the  potential 

FIG.  3  is  a  graph  demonstrating  the  cross-section  of  the 
<£=0  electrode  at  various  positions  along  the  beam  axis; 

FIG  4  is  a  graph  demonstrating  the  cross-section  of  the 
<I>=V  electrode  at  various  positions  along  the  beam  axis.; 

FIG.  5  is  a  schematic  diagram  demonstrating  the  electrode 
geometry  of  a  well-confined,  parallel  beam  of  elliptic  cross 
section; 

FIG.  6  is  a  schematic  diagram  of  a  non-axisymmetric 
periodic  magnetic  field; 

FIG.  7  is  a  schematic  diagram  of  the  field  distribution  of 
a  non-axisymmetric  periodic  magnetic  field; 

FIG.  8  is  a  schematic  diagram  demonstrating  the  labora¬ 
tory  and  rotating  coordinate  systems; 

FIGS.  9A-9E  are  graphs  demonstrating  matched  solutions 
of  the  generalized  envelope  equations  for  a  non-axisymmet¬ 
ric  beam  system  with  parameters  corresponding  to:  ko*  3.22 
cm-1,  k^-5.39  cm”1,  >/k^=0.805  cm”1,  K  1.53xlO”2  and 
axial  periodicity  length  S=0.956  cm; 
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FIGS.  10A-10E  are  graphs  demonstrating  the  envelopes 
and  dow  velocities  for  a  non-axisymmetric  beam  system 
with  the  choice  of  system  parameters  corresponding  to: 
k^-3.22  cm"1,  kov  5.39  cnf\  vic-0.805  cm’1,  K  1.53x 
10”2,  axial  periodicity  length  S -0.956  cm,  and  a  slight 
mismatch; 

FIG.  11  is  a  graph  demonstrating  the  focusing  parameter 
for  a  periodic  quadrupole  magnetic  field, 

FIG.  12  is  a  graph  demonstrating  the  beam  envelopes  of 
a  pulsating  elliptic  heam  equihhnum  state  in  the  periodic 
qmidrupole  magnetic  held  shown  in  FIG.  II; 

FIG.  13  is  a  graph  demonstrating  the  focusing  parameter 
for  a  non-axisymmetric  periodic  permanent  magnetic  field1 
and 

FIG  14  is  a  graph  demonstrating  the  beam  envelopes  of 
an  elliptic  beam  equilibrium  state  in  the  non-axisymmetric 
periodic  permanent  magnetic  field  shown  in  FIG.  13. 

DETAILED  DESCRIPTION  OF  THE 
INVENTION 

The  invention  comprises  a  non-axisymmetric  charged- 
particle  beam  system  having  a  novel  design  and  method  of 
design  for  non-axisymmetric  charged-particle  diodes. 

A  non-axi symmetric  diode  2  is  shown  schematically  in 
FIGS.  1A-1C.  FIG.  1A  shows  the  non-axi  symmetric  diode  2 
with  a  Child-Langmuir  electron  beam  8  with  an  elliptic 
cross-section  having  an  anode  4  and  cathode  6  FIG.  IB  is 
a  vertical  eross-seclional  view  of  the  non-axi  symmetric 
diode  2  and  FIG.  1C  is  a  horizontal  cross-sectional  view'  of 
the  non-axisymmetric  diode  2  showing  an  electron  beam  8 
and  the  cathode  6  and  anode  4  electrodes. 

lhe  electron  beam  8  has  an  elliptic  cross  section  and  the 
characteristics  of  Child-Langmuir  flow.  The  particles  are 
emitted  from  the  cathode  6  and  accelerated  by  the  electric 
field  between  the  cathode  6  and  anode  4.  For  an  ion  beam, 
the  roles  of  cathode  and  anode  arc  reversed. 

To  describe  the  method  of  designing  an  non-axisymmetric 
diode  with  an  elliptic  cross-section,  one  can  introduce  the 
elliptic  coordinate  system  (§,r|,z;  0*  defined  in  terms  of  the 
usual  Cartesian  coordinates  hy 

x  f  cos  h  (i-)cos(Tp,  y*f  sia  h(|)s»i(r|),  z-z  (1.1) 

where  §£[0,  so)  is  a  radial  coordinate,  r|e[0f2;t)  is  an  angular 
coordinate,  and  f  is  a  constant  scaling  parameter.  A  charged- 
particle  beam  flowing  in  the  ez  direction  and  taking  the 
Child-Langmuir  profile  of  parallel  flow  with  uniform  trans¬ 
verse  density  will  possess  an  internal  electrostatic  potential 
of 


n-  z)  = 


(1.2) 


where  one  can  have  defined  fh(z-O) r  0  along  a  planar 
charge-emitting  surface  and  <F(z-d)-V  along  a  planar 
charge-accepting  surface. 

If  both  planes  have  transverse  boundaries  of  elliptic 
shape,  specified  by  the  surface  <;=l;0=constant,  then  a  solu¬ 
tion  exists  for  a  parallel  flow,  uniform  transverse  density, 
Child-Langmuir  chargcd-particlc  hcam  of  elliptic  cross- 
section,  flowing  between  the  planes  at  z  0  and  z~d.  Due  to 
the  muuial  space-charge  repulsion  of  the  particles  constitut¬ 
ing  the  heam,  this  Child-Umgmuir  profile  must  be  supported 
hy  the  imposition  of  an  external  electric  field  through  the 


4 

construction  of  appropriately  shaped  electrodes.  The  design 
of  said  electrodes  requires  knowledge  of  the  electrostatic- 
potential  function  external  to  the  beam  which  satisfies 
appropriate  boundary  conditions  on  the  beam  edge 


to 


ic°- 


0.3' 


As  the  potential  and  its  normal  derivative  are  specified 
independently  on  the  surface  £=£0,  this  forms  an  elliptic 
15  Cauchy  prohlem,  for  which  standard  analytic  and  numerical 
solution  methods  fail  due  to  the  exponential  growth  of  errors 
which  is  characteristic  of  all  elliptic  Cauchy  problems.  The 
present  technique  builds  on  the  2-dimensional  technique  of 
Radley  in  order  to  formulate  a  method  of  solution  lor  the  lull 
20  3D  problem  of  determining  the  electrostatic  potential  out¬ 
side  a  Child-Langmuir  charged-particle  beam  of  elliptic 
cross-section. 

In  the  region  external  to  the  beam,  the  potential  satisfies 
Laplace's  equation,  which  is  written  in  elliptic  coordinates 
25  as 


30 


40 


0  = 


1 


^(C  tj,  z) 

/■'(£,  77,  z)[  f2[cosh  2f  -cos  2q){d^2  ’’  dt]: 
d2 


dz: 


/TC  7?,  z) 


1 


£_ 

Z(Z)T{f.t})[fHcosh  -cos  27;)^tf£2  +  dtj2 

a 2  | 

t-j  Z(Z)7\£,  77) 
dzi 


1 


T(f,ri)  /2(cos/i  2^-cos  2q)(d# 


F  a1  \ 


11.4) 


45 


—  —7 

zm  az 2  u 


where  one  can  follow  the  usual  technique  of  separation  of 
50  variables,  writing  F(^,r),z)^Z(z)T(H,r|)  and  introducing  the 
separation  constant  k2.  The  separated  equations  can  now  be 
written  as 


55 

(1.5) 

1  ( a1  a2  \  ,  /z 

11.6) 

60 

1  a2  k2/2 

=  mWm*—c<Kh74* 

65 
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where  one  can  have  performed  another  separation  on  the 
transverse  equatiou,  writing.  T(|,r|)-  R(5)©(rj)  and  introduc¬ 
ing  the  separation  constant  a  This  last  equation  thus  yields 


R^-^a-iJ-Cosh  2^j*(£). 

(1.7) 

■G(rj)  +-  -  ^-cos 

(1.8) 

6 

sets  of  eigenvalues  denoted  by  a2„,  a2w+1,  b2„,  b2„+1  lor  N 
e{(),  1,2,  .  }  which  differ  in  their  symmetry  properties. 

Only  the  set  a2„  and  the  corresponding  cosine-elliptic  solu¬ 
tions  denoted  by  0(rj)  cc2„(r|'k)  possess  the  appropriate 
5  symmetries,  and  the  integral  over  a  becomes  a  sum  of  the 
form 


n=0 


Solutions  to  the  separated  transverse  equations  are  known 
as  the  Radial  Mathieu  functions  R(£)  and  Angular  Mathieu 
Functions  0(rj),  respectively,  while  the  solutions  to  the 
separated  longitudinal  equation  are  easily  expressed  in  terms 
of  exponentials  7.{z)*Q*kz . 

The  solution  for  the  potential  is  now  represented  as  a 
superposition  of  separable  solutions  which,  jointly,  satisfy 
the  boundary  conditions  on  ^  One  can  write 

dkfAik^j  B{a)Ra(^k)Sa(i]:k)daJ  (1.9) 


Moreover,  the  set  of  solutions  ce2„  is  orthogonal  and 
complete  over  the  space  of  functions  with  the  desired 
symmetry  and  periodicity  properties.  Thus  one  can  expand 
unity  as 


cez 


ccznin I  k)dr} 


tU.V' 


where  the  amplitude  functions  A(k)  and  B(a)  are  introduced 
and  the  integration  contours  are  as  yet  unspecified.  In  order 
to  satisfy  the  boundary  condition  on  <I>  along  the  beam  edge  25 
using  the  analytic  continuation  of  the  Gamma  function,  one 
can  write 


The  boundary  condition  on  <t>  is  then  satisfied  by  choosing 


2  sJA4)Jc 


= 


Jo 


k)drj 


■HMt) 


A  7/3  cik. 


(1.10)  30  k)\2dtj 

and 

1. 


(1.16 


(1.171 


where  the  integration  contour  C  is  taken  around  the  branch 
cut  as  shown  in  FIG.  2 

One  can  then  write  the  boundary  condition  as 


The  condition  that  the  normal  derivative  of  the  potential 
vanishes  along  the  beam  surface  implies 


{Z\4 

«)*  v[^j 


7Z*)fc 


’kzk  113  dk 


■hm?  )- 

=  j^dk^A (k)ekz  ^B{a)^«(fo;*)@fi(^;A)rfa]. 


(1.11)  40 


=  0, 


■<-<o 


(1.18' 


45  which,  along  with  the  boundary  value  of  Rdi  and  the  eigen¬ 
value  a^,,  fully  specify  the  second-order  ftadial  Mathieu 
Equation.  It  can  then  be  integrated  by  standard  methods  in 
order  to  determine  the  radial  solutions 


The  boundary'  condition  is  satisfied  by  choosing  C  as  the 
integration  contour  for  the  representation  of  and  making 
the  correspondences 


Mk)  = 


and 


/ 


R(a)Ra(£o,  k)Ga(r);  k)rfa  =  1. 


(1.12) 

55 


(1.13) 

60 


llius,  one  may  rewrite  the  expansion  for  as 


v<r1/5  i 

<t>(£  '/«  2)  =  — 4T 


(1.19* 


•HM?) 


r 

I  re2/i(rj:k)drj 

r 

k  7/3 e*:< 

)  - —  ) 

^  X  (ln\ 

The  physical  system  requires  a  solution  periodic  in  r)  and 
symmetric  about  r|=0  and  r|  Jc/2.  In  general,  the  Angular 
Mathieu  Functions  0a(r|)  are  not  periodic.  Indeed,  a  periodic  65 
solution  arises  only  for  certain  characteristic  eigenvalues  of 
the  separation  constant  a.  There  are  4  infinite  and  discrete 


A  number  of  methods  may  be  used  to  evaluate  the  charac¬ 
teristic  values  a2w  and  the  corresponding  Angular  Mathieu 
Functions  ce2„.  These  can  be  integrated  by  standard  meth¬ 
ods.  In  practice,  only  the  first  few  terms  of  the  infinite  series 
need  be  retained  in  order  to  reduce  fractional  errors  to  below 
10"' .  The  integral  along  the  contour  C  can  be  transformed 
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into  definite  integrals  of  complex-valued  functions  along  the 
real  line,  and  thus  it,  too.  can  he  evaluated  using  standard 
methods. 

Once  the  potential  profile  is  known,  one  can  employ  a 
root-linding  technique  in  order  to  determine  surfaces  along 
which  one  may  place  constant-potential  electrodes.  A 
numerical  module  has  been  developed  which  determines 
these  electrode  shapes  based  on  the  theory  described  and 
solution  methods  described  above.  Sample  electrode  designs 
are  shown  in  FIGS.  3  and  4  for  the  case  of  a  10:1  elliptical 
heam  of  semi-major  radius  6  mm  and  semi-minor  radius  0.6 
mm.  These  electrodes  serve  to  enforce  the  analytically- 
derived  potential  profile  along  the  beam  edge,  which  in  turn 
serves  to  confine  the  beam  and  maintain  it  in  the  Child- 
Langmuir  form. 

The  3 -dimensional  charged-particle  optics  tool  Omni- 
Trak  has  been  used  to  simulate  the  emission  and  transport  of 
charge  particles  in  the  geometry'  of  FIGS.  3  and  4  The 
resulting  particle  trajectories,  shown  in  FIG.  5,  are  indeed 
parallel,  as  predicted  by  the  theory.  The  results  of  the 
Omni-Trak  simulation  also  provide  a  validation  of  the 
analytical  method  presented  above. 

One  will  ol\en  wish  to  extract  this  heam  and  inject  it  into 
another  device  by  excising  a  portion  of  the  charge-collecting 
plate  Doing  so  will  modify  the  boundary  conditions  of  the 
problem  such  that  the  above  solution  can  no  longer  be 
considered  exact,  however,  the  errors  introduced  by  rela¬ 
tively  small  excisions  will  be  negligible,  and  the  appropriate 
electrode  shapes  w  ill  be  substantially  unchanged  from  those 
produced  by  the  method  outlined  above. 

It  should  also  be  noted  that  additional  electrodes,  inter¬ 
mediate  in  potential  between  the  cathode  and  anode,  may  be 
added  in  order  to  aid  the  enforcement  of  the  Child-Langmuir 
How  condition.  The  above  prescription  allows  for  their 
design.  As  with  the  charge-collecting  plate,  neither  the 
cathode  electrode  nor  the  intermediate  electrodes  need  be 
extended  arbitrarily  close  to  the  beam  edge  in  order  to 
enforce  the  Child-L  angmuir  flow  condition.  The  portion  of 
these  electrodes  nearest  the  beam  may  be  excised  without 
substantially  affecting  the  beam  solution. 

Along  similar  lines,  in  a  physical  device,  one  cannot 
extend  the  electrodes  infinitely  far  in  the  transverse  direc¬ 
tions.  The  analytically-prescribed  electrodes  correspond  to 
the  surfaces  of  conductors  separated  hy  vacuum  and^r  other 
insulating  materials  and  (in  some  region  distant  from  the 
heam)  deviating  from  the  analytical ly-prescrihed  profiles 
Nevertheless,  as  the  influence  of  distant  portions  of  the 
electrodes  dimmish  exponentially  with  distance  from  the 
beam  edge,  these  deviations  will  have  a  negligible  effect  on 
the  beam  profile,  provided  that  they  occur  at  a  sufficient 
distance  from  the  beam  edge. 

FIG.  5  depicts  an  Omni-Trak  simulation  in  which  the 
finiteness  of  the  electrodes  is  evident  without  affecting  the 
parallel-flow  of  the  charged  particle  beam.  Note  FIG.  5 
illustrates  the  charge  collection  surface  10,  charge  emittmg 
surface  14.  parallel  particle  trajectories  12,  and  analytically 
designed  electrodes  16.  By  equating  the  electrode  geometry 
w  ith  equipotential  surfaces,  the  analytic  method  of  electrode 
design  detailed  herein  specifies  the  precise  geometry  of  the 
charge-emitting  14  and  charge-collecting  10  surfaces  as  well 
as  the  precise  geometry  of  external  conductors  16.  These 
external  conductors  may  be  held  at  any  potential,  however, 
generally,  two  external  conductors  are  used — one  held  at  the 
emitter  potential  and  the  other  at  the  collector  potential.  A 
charged -particle  system  designed  conformally  to  this  geom¬ 


etry  will  generate  a  high-quality,  laminar,  parallel-flow. 
Child-I angmuir  beam  of  elliptic  cross-section  as  show  n  in 
FIG.  5 

As  an  illustrated  example,  a  non-axisymmctric  periodic 
5  magnetic  field  for  focusing  and  transporting  a  non-ax lsym- 
mctric  beam  is  shown  FIG.  6.  FIG.  6  shows  the  iron  pole 
pieces  18  and  magnets  19  used  to  form  the  periodic  magnetic 
field.  The  iron  pole  pieces  are  optional  and  may  be  omitted 
in  other  embodiments.  The  period  of  the  magnetic  field  is 
10  defined  by  the  line  20  The  field  distribution  is  illustrated 
FIG.  7.  Note  FIG.  7  illustrates  the  field  lines  form  by  the  iron 
pole  pieces  18  and  magnets  19  of  FIG.  6 

For  a  high-brightness,  space-charge-dominated  beam, 
kinetic  (emittance)  effects  are  negligibly  small,  and  the  beam 
15  can  be  adequately  described  by  cold-fluid  equations.  In  the 
paraxial  approximation,  the  steady-state  cold-lliiid  equations 
for  time-stationary  flow  (d/dt=0)  in  cgs  units  are: 


20  d 

Pt>fTnb  +  V  .  •  infrVj  =  0, 
os 

P  =  /£*  V‘  A]  =  -4 nqnb, 

(  d  d  \ 

25  4*si +  V‘- Mir- 

— -  -  -i  V  +  fat,  X  B?  +  —  x 

W  n  c 


30  where  s-z,  q  and  m  are  the  particle  charge  and  rest  mass, 
respectively. 


l 

v  i-$ 


is  the  relativistic  mass  factor,  use  has  been  made  of 
const,  and  the  self-electric  field  E*  and  self-magnetic 
40  field  W  are  determined  from  the  scalar  potential  <ps  and 
vector  potential  A/e?,  i.e.,  EJ— and  LT  VxA/e,. 

One  seeks  solutions  to  Eqs.  (2.1)-(2.3)  of  the  form 


nb{x^,  s)  =  ■ 


-0 


x 2  y2  1 

<72(5)  b~(S)  J* 


(2.4) 


na{s  )£(5) 

s)  =  1/1,  (5)x  -  +  LfM-sly  +  ay{s)x]fari}r  t2.5) 


50 

In  Eqs.  (2,4)  and  (2.5),  Xy^xe^+yCp  is  a  transverse  displace¬ 
ment  in  a  rotating  frame  illustrated  in  FIG.  8,  0(s)  is  the 
angle  of  rotation  of  the  ellipse  with  respect  to  the  laboratory 
frame;  0(x)=l  if  x>0  and  0(x)^O  if  x<0;  and  the  functions 
55  a(s).  b(s),  m(s),  pv(s),  at(s),  av.(s )  and  0(s)  are  to  be 
determined  self-consistently  [see  Eqs.  (2.1 1  )-(2.1 5)]. 

For  the  self-elcctric  and  self-magnetic  fields,  Eqs.  (2.2) 
and  (2  4)  arc  solved  to  obtain  the  scalar  and  vector  potentials 


60 


a  +  cj  b  ) 


(2.6) 


65  For  a  3D  no n -a xi symmetric  periodic  magnetic  field  with 
an  axial  periodicity  length  of  S,  one  can  describe  it  as  the 
fundamental  mode, 
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-continued 


Bo  U—  sinh(koxx)cosh{koy  y)cos{kos)ex  +  -7-^co$/i(Ao.**HinMAo>.y) 
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2V*7cos(A<>y) 


aLcos2#  ,  , —  db 

- - b‘  -  2'J*z  b—sin(kos)  =  0. 

*0  as 


cos (kos]ey  -  cos/»(Ao,*)cos/i(Aoy>,)$in(A(>s)2i 


and  further  expand  it  to  the  lowest  order  in  the  transverse 
dimension  to  obtain 


dG  a2ay  -  tfa x 


a2-*?2 


,  where 


* 


<?Bq 
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In  Eqs.  (2.7)  and  (2.8), 


2 n 


A0  -  ■ 


An,  +  A0k,  —  Ay. 


(2.9)  35 


(2.10) 


Equations  (2.1 1)-(2.1 5)  have  the  time  reversal  symmetry 
under  the  transformations  (s,a,b,a\b\av,av,0)— * ►(-s,a,b,-a\- 
b1, -0^,-0^ 0).  This  implies  that  the  dynamical  system 
described  by  Eqs.  (2  1 1  )-(2  1 5)  has  the  hyper  symmetry 
plane  (a*  b’^c^,). 

A  numerical  module  was  developed  to  solve  the  gener¬ 
alized  envelope  equations  (2.11)-(2.15).  There  are.  in  total, 
seven  functions  a(s).  b(s),  a’(s),  b’(s),  a,(s).  at(s)  and  0(s)  to 
be  solved.  The  time  inverse  symmetry  of  the  dynamical 
system  requires  the  quantities  (a\b\(xx,av)  vanish  at  s~0  for 
matched  solutions,  therefore,  only  the  three  initial  values 
a(0),  b(0)  and  0(0)  corresponding  to  a  matched  solution  need 
to  be  determined  by  using  Newton’s  method.  Hie  matched 
solutions  of  the  generalized  envelope  equations  are  shown  in 
FIGS.  9A-9E  for  a  non-axi symmetric  beam  system  with  the 
choice  of  system  parameters  corresponding  to:  1 3.22 
cm-1,  k^,-5.39  cm"1,  vd<-0.805  cm”1,  K  1.53x10" 2  and 
axial  periodicity  length  S=0.956  cm. 

In  particular,  FIG.  9 A  demonstrates  the  envelopes  asso¬ 
ciated  with  the  functions  a(s)  and  b(s).  FIG.  9B  is  graphical 
representation  of  rotating  angle  0(s).  FIG.  9C  is  a  graph 
illustrating  velocity 


M s)  =  : 


Idle  31)  magnetic  field  is  specified  by  the  three  parameters 
B0,  S  and  k^/kp^. 

Using  the  expressions  in  Eqs.  (2.5),  (2.6)  and  (2.8),  it  can 
be  shown  that  both  the  equilibrium  continuity  and  force 
equations  (2.1)  and  (2  3)  are  satisfied  if  the  dynamical 
variables  a(s),  b(s),  gx(s)»a~Ida/ds,  pv(s)»b-1db/ds,  a*(s), 
ar(s)  and  0(s)  obey  the  generalized  beam  envelope  equa¬ 
tions: 


FIG.  9D  is  a  graph  demonstrating  velocity 
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IK 

b - -  =0. 

a  +  b 


d  ,  ab\(*x  -av)  d  ,  a  \ 
ds  a^-b2  ds\b' 

1 2  r.n*2a  . 


, —  AqjCOs2#  +  Afl  sin‘0  —  da 

2vk,  cos(Aoi) - 7 - a  -  2 v*v  a  —  ^in{kos)a  =  0, 

A0  ds 


FIG.  9E  is  a  graph  demonstrating  velocities  oq.(s)  and  av(s) 
50  versus  the  axial  distance  s  for  a  fiat,  ellipse-shaped,  uniform- 
density  charged -particle  beam  in  a  3D  non-axisymmelrie 
magnetic  field. 

The  matching  from  the  charged-particle  diode  to  the 
focusing  channel  might  not  be  perfect  in  experiments.  If  a 
=  0,  55  mismatch  is  unstable,  it  might  ruin  the  beam.  However. 

investigations  of  small-mismatch  beams  show  that  the  enve- 
(2.12)  lopes  are  stable  against  small  mismatch. 

For  example,  the  envelopes  and  flow  velocities  are  plotted 
in  FIGS.  10A-10E  for  a  non-axisymmetric  beam  system 
60  with  the  choice  of  system  parameters  corresponding  to: 
^=3.22  cm"1,  ^=5.39  cm"1,  ^=0.805  cm"1,  K 
1 .53  xlO"2  and  axial  periodicity  length  S  0.956  cm  w  ith  an 
initial  5%  mismatch  of  0,  i.e.  0(s=O)=Owa,cW(s=T))x(l  .05). 
In  particular,  FIG.  10A  demonstrates  the  envelopes  asso- 
65  ciated  with  the  functions  a(s)  and  b(s).  FIG.  10B  is  graphical 
representation  of  rotating  angle  0(s).  FIG.  IOC  is  a  graph 
illustrating  velocity 


(2.13) 
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FIG.  101)  is  a  graph  demonstrating  velocity 


t  db 
b  ds ' 
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FIG  1  OB  is  a  graph  demonstrating  velocities  ax(s)  and  ay(s) 
versus  the  axial  distance  s  for  a  Hat,  ellipse-shaped,  uniform-  15 
density  charged-particle  beam  in  a  3D  non-ax isymmetric 
magnetic  field. 

By  the  technique  described  herein,  one  can  design  a 
non-axisymmetric  magnetic  focusing  channel  which  pre-  20 
serves  a  uniform-density,  laminar  charged-particle  beam  of 
elliptic  cross-section. 

One  can  illustrate  how  to  match  an  elliptic  charged- 
particle  beam  from  the  non-axisymmctric  diode,  described 
herein,  into  a  periodic  quadrupole  magnetic  Held  In  the  25 
paraxial  approximation,  the  periodic  quadrupole  magnetic 
field  is  described  by 


ir" 


( v£t  +  xiy ). 


(3.1)  30 


The  concept  of  matching  is  illustrated  in  FIGS.  11  and  12. 

1  1G.  11  shows  an  example  of  the  magnetic  focusing 
parameter 


q  (dRi) 

k,*(S)  YbPb*nc2{  dy  J0 


(3.2)  40 


associated  with  the  periodic  quadrupole  magnetic  field  for  a 
beam  of  charged  particles  with  charge  q,  rest  mass  m,  and  45 
axial  momentum  y^mc. 

FIG.  12  shows  the  envelopes  for  pulsating  elliptic  beam 
equilibrium  in  the  periodic  quadrupole  magnetic  field,  as 
described  previously.  50 

Hie  matching  of  the  equilibrium  state  from  the  diode  to 
the  equilibrium  state  for  the  periodic  quadrupole  magnetic 
field  at  s  0  is  feasible,  because  the  transverse  density  profile 
and  How  velocity  of  the  two  equilibrium  states  are  identical 
at  s^O.  In  particular,  the  transverse  particle  density  is  uni-  55 
form  within  the  beam  ellipse  and  the  transverse  How  veloc¬ 
ity  vanishes  at  s~0. 

Also,  one  can  illustrate  how  to  match  an  elliptic  charged- 
particlc  beam  from  the  non-axisymmetric  diode,  as 
described  herein,  into  a  non-axisymmetric  periodic  penna-  6 
nent  magnetic  field  In  the  paraxial  approximation,  the 
non-axisymmctric  permanent  magnetic  field  is  described  by 
Eq.  (2.8).  The  concept  of  matching  is  illustrated  in  FIGS.  13 
and  14 

FIG.  13  shows  an  example  of  the  magnetic  focusing 
parameter 


2  Jkfornc1 
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associated  with  the  non-axisymmetric  periodic  permanent 
magnetic  field  (presented  for  a  beam  of  charged  particles 
with  charge  q,  rest  mass  m,  and  axial  momentum  y^mc. 

FIG.  14  shows  the  envelopes  for  a  fiat,  elliptic  beam 
equilibrium  state  in  the  non-axisymmetric  periodic  perma¬ 
nent  magnetic  field.  The  angle  of  the  ellipse  exhibits  slight 
oscillations.  However,  these  oscillations  can  be  corrected  by 
utilizing  higher  longitudinal  harmonics  of  the  magnetic  field 
profile. 

The  matching  of  the  equilibrium  state  from  the  diode  to 
the  equilibrium  state  for  the  non-axisymmetric  periodic 
permanent  magnetic  field  at  s  -0  is  feasible,  because  the 
transverse  density  profile  and  flow  velocity  of  the  two 
equilibrium  states  are  identical.  In  particular,  the  transverse 
particle  density  is  uniform  within  the  heam  ellipse  and  the 
transverse  flow  velocity  vanishes  at  s~0. 

The  matching  procedure  discussed  herein  illustrates  a 
high  quality  interface  between  a  non-axisymmetric  diode 
and  a  non-axisymmetric  magnetic  focusing  channel  lor 
charged-particle  beam. 

fhis  beam  system  will  find  application  in  vacuum  elec¬ 
tron  devices  and  particle  accelerators  where  high  brightness, 
low  emittance,  low  temperature  beams  are  desired. 

Although  the  present  invention  has  been  shown  and 
described  with  respect  to  several  preferred  embodiments 
thereof,  various  changes,  omissions  and  additions  to  the 
form  and  detail  thereof,  may  be  made  therein,  without 
departing  from  the  spirit  and  scope  of  the  invention. 

What  is  claimed  is: 

1 .  A  eharged-particle  beam  system  comprising 

a  non-axisymmctric  diode  that  forms  a  non-axisymmetric 
beam  having  an  elliptic  cross-section:  and 

a  focusing  channel  that  utilizes  a  magnetic  field  for 
focusing  and  transporting  said  elliptic  cross-section 
beam. 

2  The  charged-particle  beam  system  of  claim  of  1, 
wherein  said  charged-particle  beam  possesses  a  uniform 
transverse  density. 

3.  The  charged-particle  beam  system  of  claim  of  1, 
wherein  said  charged-particle  beam  possesses  a  laminar 
flow. 

4.  The  chargcd-particle  beam  system  of  claim  of  1, 
wherein  said  charged-particle  beam  possesses  a  parallel 
longitudinal  flow. 

5.  rhe  charged-particle  beam  system  of  claim  1,  wherein 
said  focusing  channel  comprises  a  non-axisymmctric  mag¬ 
netic  field  for  focusing  and  transporting  said  chargcd-par¬ 
ticle  beam. 

6.  The  charged-particle  beam  system  of  claim  5,  wherein 
said  non-axisymmetric  magnetic  field  includes  a  non-axi¬ 
symmctric  periodic  magnetic  field. 

7.  The  charged-particle  beam  system  of  claim  5  wherein 
said  non-axisymmetric  magnetic  field  includes  a  non-axi¬ 
symmetric  permanent  magnetic  field. 

8.  The  chargcd-particlc  beam  system  of  claim  5,  wherein 
said  non-axisymmetric  magnetic  field  includes  a  non-axi¬ 
symmetric  periodic  permanent  magnetic  field. 

9  The  charged-particle  beam  system  of  claim  of  5, 
wherein  said  non-axisymmctric  magnetic  field  includes  at 
least  one  quadrupole  magnetic  field. 
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10.  The  charged- particle  beam  system  of  claim  of  5. 
wherein  said  non-ax i symmetric  magnetic  field  includes  a 
periodic  quadrupole  magnetic  field. 

11.  The  charged-particle  beam  system  of  claim  2,  w  herein 
said  focusing  channel  comprises  a  non-axisymmetric  mag¬ 
netic  field  for  focusing  and  transporting  said  charged-par¬ 
ticle  beam. 

12.  The  charged -particle  beam  system  of  claim  11, 
wherein  said  non-axisymmetric  magnetic  field  includes  a 
non-axisymmetric  periodic  magnetic  field. 

13.  The  charged-particle  beam  system  of  claim  11 
wherein  said  non-axisymmetric  magnetic  field  includes  a 
non-axisymmetric  permanent  magnetic  field. 

14.  The  charged -particle  beam  system  of  claim  11, 
wherein  said  non-axisymmetric  magnetic  field  includes  a 
non-axisymmetric  periodic  permanent  magnetic  field. 

15.  The  charged- particle  beam  system  of  claim  of  11 
wherein  said  non-axisymmetric  magnetic  field  includes  at 
least  one  quadrupole  magnetic  field. 

16.  The  chargcd-particle  beam  system  of  claim  of  11 
wherein  said  non-axisymmetric  magnetic  field  includes  a 
periodic  quadrupole  magnetic  field. 

17.  A  method  of  forming  a  charged-particle  beam  system 
comprising 

forming  a  non-axisymmetric  diode  that  includes  a  non- 
axisymmetric  beam  having  an  elliptic  cross-section 
and 

forming  a  focusing  channel  that  utilizes  a  magnetic  field 
for  focusing  and  transporting  said  elliptic  cross-section 
beam. 

18.  The  method  of  claim  17,  wherein  said  charged-particle 
beam  possesses  a  uniform  transverse  density. 

19.  The  method  of  claim  1 7,  wherein  said  charged-partiele 
beam  possesses  a  laminar  flow. 

20  The  method  of  claim  1 7,  wherein  said  charged -particle 
beam  possesses  a  parallel  longitudinal  flow. 

21  The  method  of  claim  17,  wherein  said  focusing 
channel  comprises  a  non-axisymmetric  magnetic  field  for 
focusing  and  transporting  said  charged-particle  beam. 
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22.  'lhe  method  of  claim  21,  wherein  said  non-axisyin- 
metnc  magnetic  field  includes  a  non-axisymmetric  periodic 
magnetic  field. 

23.  The  method  of  claim  21,  wherein  said  non-axisym- 
metrie  magnetic  field  includes  a  non-axisymmetric  perma¬ 
nent  magnetic  field. 

24.  lhe  method  of  claim  21,  wherein  said  non-axisym¬ 
metric  magnetic  field  includes  a  non-axisymmetne  periodic 
permanent  magnetic  field. 

25.  The  method  of  claim  21,  wherein  said  non-axisym¬ 
metric  magnetic  field  includes  at  least  one  quadrupole 
magnetic  field 

26.  The  method  of  claim  21,  wherein  said  non-axisym¬ 
metric  magnetic  field  includes  a  periodic  quadrupole  mag¬ 
netic  field. 

27.  The  method  of  claim  18,  wherein  said  focusing 
channel  comprises  a  non-axisymmetric  magnetic  field  for 
focusing  and  transporting  said  charged-particlc  beam 

28.  The  charged-particle  beam  system  of  claim  27, 
wherein  said  non-axisymmetrie  magnetic  field  includes  a 
non-axisymmetric  periodic  magnetic  field. 

29  lhe  charged-particle  beam  system  of  claim  27, 
wherein  said  non-axisymmetric  magnetic  field  includes  a 
non-axisymmetric  permanent  magnetic  field 

30.  The  charged-particle  beam  system  of  claim  27, 
wherein  said  non-axisymmetric  magnetic  field  includes  a 
non-axisymmetne  periodic  permanent  magnetic  field. 

31  The  chargcd-particle  beam  system  of  claim  of  27. 
wherein  said  non-axisymmetric  magnetic  field  includes  at 
least  one  quadrupole  magnetic  field. 

32.  The  chargcd-particle  beam  system  of  claim  of  27. 
wherein  said  non-axisymmetric  magnetic  field  includes  a 
periodic  quadrupole  magnetic  field. 
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